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Einleitung 



Die Kleintierbildgebung ist ein wichtiger Bereich der Grundlagen- und präklinischen 
Forschung, in dein verschiedene Wissenschaftsdisziplinen wie Medizin. Biologie. Chemie. 
Physik und Medizintechnik eng zusammenarbeiten. Die voranschreitende Entwicklung 
transgener Mausmodelle bietet vielfältige Möglichkeiten biomedizinische Fragestellungen 
präzise im lebenden Organismus zu untersuchen. Stoffwechselfunktionen lassen sich 
in diesen Tiermodellen mittels gentechnischer Manipulation gezielt beeinflussen oder 
ausschalten. Von großem Interesse sind dabei Studien, bei denen die lebenden Tiere 
über einen langen Zeitraum immer wieder untersucht werden können, um so funk- 
tionelle Abläufe und deren Veränderungen in einem komplexen Zusammenhang zu 
studieren. Die Visualisierung dieser Prozesse im lebenden Tier ist die Aufgabe der 
bildgebenden Verfahren. Viele der Verfahren, die ursprünglich für den klinischen Einsatz 
entwickelt wurden, sind mittlerweile mit miniaturisierten Geräten an die Anforderungen 
der Kleintierbildgebung angepasst. Das Bindeglied zwischen der Visualisierung und 
den tatsächlichen funktionellen Vorgängen in einem Mausmodell ist die Kopplung von 
Signalgebern an die zu beobachtenden Abläufe. Die herausfordernde Aufgabe der Bio- 
chemie besteht darin, geeignete Biomoleküle zu finden und diese für die Darstellung 
mit der entsprechenden Modalität zu präparieren. Ausführliche Ubersichtsartikel mit 
detaillierter Gegenüberstellung der bildgebenden Verfahren und deren Anwendungen in 
der Kleintierbildgebung linden sich in [1. 2. 3. 4j. 

Neben der fortschreitenden Entwicklung in den einzelnen Modalitäten verstärkt sich 
innner mehr der Trend, unterschiedliche Verfahren miteinander zu kombinieren und 
von den Vorteilen der jeweiligen Modalität während einer simultanen Bildgebung zu 
profitieren. Für die Untersuchung am Mausmodell verspricht z. B. die überlagerte 
Darstellung von Anatomie und funktionellen Abläufen oder die Kombination schneller 
und langsanier Bildgebung einen erheblichen Informationsgewinn. 
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1. Einleitung 



1.1 Motivation und Schwerpunkt der Arbeit 

Die Integration unterschiedlicher bildgebender Modalitäten erfordert speziell angepasste 
und optimierte Mosssvstemo. mit denen die jeweiligen Signalintensitäten möglichst un- 
gestört detektiert werden können. Das Prinzip eines Hachen, optischen Detektors wurde 
bereits in experimentellen Studien und mit kombinierter PET-Bildgobung 1 gezeigt [.">. 
Für die optische Bilderzeugung wurde zu diesem Zeitpunkt jedoch nur ein kleiner Teil 
der verfügbaren Detektorinformationen ausgewertet (Fokuspunktbilder). Diese Arbeit 
knüpft nun unmittelbar an die zuvor gewonnenen Erkenntnisse über dieses Detektor- 
prinzip an und beinhaltet die weiterführende Forschung mit Blick auf die optischen 
Bildgebung. Das Ziel dieser Arbeit ist. optische Mikrolinsen-Detektoren für die kon- 
taktlose überflächensignalmessung anstelle einer herkömmlichen Kamera einzusetzen, 
um mit minimalem Platzaufwand zukünftig optische Tomographie für kleine Objekte 
zu ermöglichen. 

Die Schwerpunkte der Arbeit lassen sich hierfür in zwei Etappen unterteilen: 

• Hildaehuna mit einem optischen Detektor: Im ersten Schritt soll das Defekt orprinzip 
in einen einsatzbereiten, optischen Detektor umgesetzt werden. Im Mittelpunkt 
steht dabei die Entwicklung eines Abbildungsmodells, das die bisher verworfenen 
Detektorinformationen einbindet und somit die Bildgebung gegenüber der Fokus- 
punktmethode verbessert. Ausgehend von den vorhandenen Komponenten soll ein 
robuster, kalibrierter Detektor zusammengesetzt werden, dessen Abbildungseigen- 
sehaften anschließend zu charakterisieren sind. 

• Aufbau eines optischen Kleintiers canners mit mehreren Detektoren: Im zweiten 
Schritt ist vorgesehen, mehrere optische Detektoren auf einer rotierbaren Plattform 
anzuordnen und damit einen optischen Kleintierscanner zu realisieren. Das Ziel 
ist hier, über eine platzeffiziente Detektoranordnung eine Nacktmaus großflächig 
abzudecken. In experimentellen St udien soll dann die Oberflächensignahnessung 
durchgeführt und bewertet werden. 

Die Arbeit gliedert sich in eine Einleitung, einen Methodenteil, ein Ergebniskapitel sowie 
eine absehli« ßende Diskussion, eine Zusammenfassung und ein- -:i Ausblick, Im Mötho- 
denteil (Kap. 2) wird der Aufbau und das Rahmenkonzept eines Detektormoduls, das 
Detekt.o-Abl ildutigsmodell. die Methoden ffu di< Ch.ir.ikvüsi.-mug der Abbildungsei- 
genschaften und der experimentelle Seaimeraufbau beschrieben. Die Ergebnisse sind 
in Kap. 3 getrennt nach der Bildgebung mit einem Detektor und nach den Messungen 

1 Posit runen-Emissions- Tomographie 
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1.2. Optische Tomographie 



a 



mit dein Seaimeraufhau darfst eilt . In der Diskussion in Kap. 4 werden die Ergebnisse 
schließlich bewertet. Die Arbeit ist in Kap. 5 zusammengefasst und ein Ausblick ist in 
Kap. 6 gegeben. 

1.2 Optische Tomographie 

Optische Methoden sind in der Kleintierbildgebung weit verbreitet und routinemäßig im 
Einsatz. Gründe dafür sind die fehlende ionisierende Strahlung, die hohe Kosteneffizienz 
und die gute Verfügbarkeit der optischen Farbstoffe. Die Anwendungen der optischen 
Bildgebung reichen von schnellen, planaren Aufnahmen bis hin zur tomographischen 
Rekonstruktion. 

Die optische Tomographie ist ein vielversprechendes Werkzeug für die Kleintierbild- 
gebung und ein aktives Forschungsfeld, sowohl in Bezug auf die Rekonstruktionsalgo- 
rithmen als auch auf die erforderliche Signalmessung. Einen umfassenden Überblick zur 
optischen Tomographie bieten die Übersichtsartikel in [7. 8. 9. 10]. Ausführliche Beschrei- 
bungen zur mathematischen Behandlung der Rekonstruktion finden sich z. B. in [11, 12]. 
Grundlage der optischen Tomographie ist die Lichtausbreitung im Untersuchungsobjekt, 
durch die auf die optischen Gewebeeigenschaften (endogener Kontrast) oder die Konzen- 
trationsverteilung eines optisch aktiven Signalgebers (exogener Kontrast) geschlossen 
werden soll. Im Gegensatz zu den hochenergetischen Röntgenphotonen (~ IkeV bis 
150 keV), die Gewebe nahezu ungestreut passieren, können für die niederenergetische 
elektromagnetische Strahlung (~ leV bis 2,5 eV) des sichtbaren Lichtes die Streuei- 
genschaften des Gewebes nicht mehr vernachlässigt werden. Die Herausforderung der 
tomographischen Rekonstruktion liegt nun darin, die Lichtausbreitung im Gewebe ma- 
thematisch adäquat zu beschreiben, um damit eine Lösung zu finden, die die tatsächliche 
Signahnessung erklärt. Eine allgemein akzeptierte, wenngleich auch nicht ganz unkri- 
tische Methode für die Kleintierbildgebung, basiert auf der Diffusionsapproximation der 
Strahlungstransport-Gleichung. Die Lösung dieses Ansatzes erfordert die Signahnessung 
an der Objektoberfiäche und somit auch die Kenntnis der Objektoberfläche selbst. 
Erschwerend kommt hinzu, dass das Problem mathematisch schlecht gestellt ist (engl.. 
ill-posed problem). d.h. es gibt eine Vielzahl von Lösungen, die sich für ein solches 
Gleichungssystem linden lassen und fraglich ist immer, ob die gefundene Lösung auch tat- 
sächlich der reellen Situation entspricht. Die Ausgangssituation dieser Problemstellung 
durch zusätzliche a-pre'on'-Information (z.B. aus der multimodalen Bildgebung) oder 
durch verschiedenen Messmodi (Zeitauflösung, Frequenzdiskriminierung) zu verbessern 
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1. Einleitung 



-•11:] ( ;>-;-i-i;i>l .Im aktuellen Fnrs.-l,in». 

Die diffuse optische Tomographie (DOT) bezeichnet allgemein die Rekonstruktion der 
optischen Streu- und Absorptionseigenschaften von Gewebe. Der endogene optische 
Kontrast des Gewebes wird z. B. durch die Absorptionseigenschaften von Zellfarbstoffen 
wie Hämoglobin hervorgerufen. Für die diffuse optische Tomographie mit exogenen Si- 
gnalgebern unterscheidet man zwischen FMT ( viuß. fluorescence mediated tomography) 
mit Fluoreszenzfarbstoffen und BLT (Biolumineszenztomographie) mit Biolumineszenz- 
markern. Für die FMT werden Fluoreszenzmarker durch eine externe Quelle angeregt, 
die dann selbst Licht mit einer längeren Wellenlänge emittieren. Bei der Signalmessung 
müssen Anregungs- und Emissionsücht getrennt werden, da deren Intensitäten zu- 
sammen in den Rekonstruktionsprozes.- .•'•n.-eh-ii. l'ber die bekannte Position der 
Anregungsquelle auf der Objektoberfiäche wird versucht, auf die unbekannte Position 
und Konzentration des Fluoreszonzmarkers im Objekt zu schließen. Biolumineszenz- 
marker dagegen emittieren Licht unter bestimmten biochemischen Voraussetzungen, 
sodass eine externe Anregung nicht erforderlich ist. 

Die Signalmessung für die optische Tomographie wird unmittelbar an der Objektober- 
fläche durchgeführt, indem die Anregungsquelle (im Falle der FMT) und die Detektoren 
durch Glasfasern mit der Oberfläche in Kontakt gebracht werden. Aufgrund der limi- 
tierten Anzahl von Glasfasern ist hier jedoch nur eine räumlich grobe Signalaufnahnie 
möglich. In jüngster Zeit hat sich immer mehr die kontaktlose Oberllächensignahnessung 
mit CCD-Kameras 2 durchgesetzt, mit der sowohl größere Datenmengen orfasst und somit 
die Oberflächenabtastung erhöht werden kann, als auch der experimentelle Aufwand 
vereinfacht wird. Die kontaktlose Signalmessung mit einer CCD-Kamera für die FMT 
ist in [12] demonstriert. 

1.3 Bilderzeugung mit Mikrolinsenarrays 

Die Funktion des refraktiven Mikrolinsenarrays (MLA) im optischen Detektor in dieser 
Arbeit ist es. analog zu einer einzelnen Objektivlinse einer klassischen Kamera, das 
Gesichtsfeld für die optische Abbildung zu definieren. Der optische Detektor besteht aus 
einem MLA für die optische Abbildung, einem CM OS-Sensor 3 in der MLA-Fokalebene 
und einer Septenmaske. die das Überlappen der einzelnen Mikrolinsenprojektionen auf 
der Sensoroberfläche verhindert . Der Detekoraufbau wird noch ausführlich in Abschn. 2.1 

2 engl, chargo coupied device 

' ! cn/y/. eompleiuentary metal oxide semiconduetor 
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1.3. Bilderzeugung mit Mikrolinsenarrays 



beschrieben. 




W (b) (c) 



Abb. 1.1. Bilderzeugiing mit einem MLA: Fokuspiinktbilder sind aus den abgetasteten 
Intensitäten der zentralen Fokuspunkte zusammengesetzt, (b) Fokuspiinktbilder, die sieh aus 
dezentralen Fokuspunkten zusammensetzten, entstehen aus einem um den entsprechenden 
Winkel geneigten Gesichtsfeld, (c) Das in dieser Arbeit entwickelte Abbildungsmodell verwendet 
gleichzeitig alle Fixelintensitäten für die Bilderzeugiing. Die Abbildung ist entnommen aus [13|. 

Eine praktikable Methode, um mit dem optischen Detektor eine Bild zu erzeugen, ist 
die Betrachtung der Fokuspunktinleusiläten wie in [5] gezeigt und in Abb. 1.1 darge- 
stellt. Parallel einfallende Lichtstrahlen, die senkrecht auf die plankonvexen Mikrolinsen 
treffen (a), werden in die zentralen Fokuspunkte gebrochen. Tastet man dann die In- 
tensitäten in den zentralen Fokuspunkten ab, so entsteht ein Bild aus einem parallel 
kollimierten Gesichtsfeld. 4 Diese Bilder werden im weiteren Verlauf als Fokuspiinktbilder 
bezeichnet, wobei immer die Abtastung der zentralen Fokuspunkte der Mikrolinsen 
auf der Sensoroberfläche gemeint ist. falls nicht anders erwähnt. Grundsätzlich muss 
ein Fokuspunktbild nicht zwingend aus den zentralen Fokuspunkten bestehen, sondern 
kann sich aus benachbarten dezentralen Fokuspunkten zusammensetzen. Dadurch wird 
das Gesichtsfeld um den entsprechenden Winkel geneigt, wie in Abb. 1.1 (b) skizziert. 
Das simultane Zusammenspiel aller Fokuspunktintensitäten (c) für die Bilderzeugung 
mit dem optischen Detektor ist schließlich die Grundlage des Abbildungsmodells. das 
in Abschn. 2.2 vorgestellt wird. Abb. 1.2 zeigt verschiedene Fokuspiinktbilder. die aus 
einem Sensordatensatz (s. Abschn. 3.1.1) extrahiert wurden. Jedes Fokuspunktbild 
besteht hier aus 51 x 101) px n . entsprechend der Anzahl der Mikrolinsen. und hat eine 



* Eine Analogie zu diesem Prinzip findet sich in den Parallelloch-Kollimatoren in der SPECT 
(engl, single photon Omission computed tomography) [14]. 

8 Ein Pixel [px] leitet sich ab von (eng! ) picture element und bezeichnet die kleinste Flä- 
cheneinheit eines digitalen Bildes, die durch eine bestimmte Position in der Bildmatrix und 
durch einen Grauwert gekennzeichnet ist. 
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1. Einleitung 



leicht unterschiedliche Perspektive auf das Objekt. Die Anordnung der Fokuspunkt- 
bilder entspricht dabei dem Otl'se! der extrahierten Fokuspunkte zur Optischen Achse, 
d.h. das zentrale Fokuspunktbild hat ein senkrecht kolliiuiertes Gesichtsfeld wie in 
Abb. l.l(a) dargestellt. Alle anderen Fokuspunkt bilder haben ein entsprechend geneigtes 
Gesichtsfeld, das sich dann in d> i l> änderten Perspektive äußert. Die Fokuspunktbilder 
in den Ecken fehlen, da hier die Sensorpixel durch die Septenmaske (s. Abb. 2.3) verdeckt 
werden. Es sei an dieser Stelle auch erwähnt, dass die parallele Abbildungsgeomtrie 
der l ; okuspunktbilder nur näherungsweise gilt, da Sensorpixel nicht punktförmig sind, 
sondern ein Ausdehnung besitzen. Jedes Fokuspunktbild wird grundsätzlich normiert. 
An den gleichen Pixelpositionen werden dafür die Intensitäten aus dein homogen ausge- 
leuchteten Detektor (Nonnierungsbild. s. Abschn. 3.1.1) extrahiert, durch die dann das 
Fokuspunkt bikl pixelweise dividiert wird. 




Abb. 1.2. Fokuspunkt bilder mit unterschiedlichem Offset zur optischen Achse. Jedes Bild 
unterscheidet sich dadurch in der Perspektive. Aus dem Sensordatensatz könnten 10 x 10 
Fokuspiinktbilder erzeugt werden, von denen jedoch einige durch die Septenmaske verdeckt 
sind. Der Übersichtlichkeit wegen ist hier nur jedes vierte Fokuspuuktbild abgebildet. Der 
schwarze Streifen durch die drei rechts liegenden Fokuspunktbilder ist durch fehlerhafte 
Sensorpixel verursacht. 

In der Literatur linden sich ebenfalls verschiedene Prinzipen für die MLA-Bilderzeugung. 
Das Prinzip des künstlichen Insektenauges (artificial apposition Compound eye) er- 
möglicht ein sehr dünnes optisches System [15, IG], das zehnmal kleiner ist als der 
optische Detektor. Jeder der optischen Kanäle besteht hier aus einer Mikrolinse und 
einer einzigen photosensitiven Fläche, deren Anordnung zueinander das Gesichtsfeld 
des Systems definiert. Durch lichtdichte Abtrennungen werden die Kanäle optisch 
voneinander isoliert. Die Bilderzeugung ist vergleichbar mit der Fokuspunktmethode. 
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Ein weiteres Konzept, aus der Photographie ist gezeigt in [17, 18|. Hier wird das M LA 

a\\ is> heu die ( ibjektiviiiis 1 und den Bildsens' r oh trag »areii Kamer« eingesetzl , an 

damit nachträglich scharfe, synthetische Bilder in unterschiedlichen Tiefen zu berechnen, 
was als digitale Refokussieruug bezeichnet wird. Das Überlappen der einzelnen Mikro- 
linsenprojektioncn auf dem Sensor wird vermieden, indem die relative Größe der Haupt- 
und Mikrolinsenaperturen aufeinander abgestimmt werden. Ebenfalls in Verbindung 
mit. einer CCD-Kamera wird im Bereich des Integral Im agnig ein MLA diesmal vor der 
Objektivlinse der Kamera verwendet [19, 20, 21]. Die aufgenommene Szene lässt sich 
dann entweder auf einem LCD-Display darstellen und durch ein MLA betrachten oder 
durch ein idgorithinisches Verfahren (computatioual integral imagiug) berechnen. Diese 
Methoden sind vergleichbar mit dem Abbildungsmodell, da hier eine Mikrolinse mit 
mehreren Sensorpixeln assoziiert ist. 
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Material und Methoden 



In diesem Kapitel werden die Expcrimontalaufbauten und Methoden vorgestellt, mit 
denen sich die Oberllächensignahnessung mit Mikrolinsen-Detektoren realisieren lässt. 
Zu Beginn wird in Absehn. 2.1 die experimentelle Anordnung und das Rahmenkonzept 
des optischen Detektors vorgestellt. Anschließend wird in Abschn. 2.2 das Detektor- 
Abbildungsmodell erklärt, mit dem sieh die Bilderzeugung gegenüber einem Fokus- 
punkt bild deutlich verbessern lässt. Um die Abbildungseigenschaften wie die räumliche 
Auflösung und das Signal- Rausch- Verhältnis zu charakterisieren, werden in Abschn. 2.3 
hierzu geeignete Methoden vorgestellt. Abschließend wird der experimentelle Scanner- 
aufbau in Abschn. 2.4 beschrieben, der aus vier Einzeldetektoren besteht. In diesem 
Zusammenhang wird auch die Scannerkalibrierung, die Oberflächenberechnung sowie 
ein Messprotokoll für die Signalmessung gezeigt. 

Jegliche Formalismen, die in diesem Kapitel vorgestellt werden, sind für die experi- 
mentelle Umsetzung in der plattformunabhängigen Programmiersprache Java™ [22] 
implementiert. Die Algorithmen lassen sich damit als Plugin in die Open-Source-Software 
bnageJ [23] einbetten, die ein sehr hilfreiches Werkzeug für die Bilddarstellung ist und 
für die Bildbearbeilung und -analyse eine umfangreiche Funktionalität, bietet. 

2.1 Experimenteller Aufbau des optischen Detektors 

Grundlage dieser Arbeit ist der optische Detektor, mit dem die Oberflächensignal- 
messung durchgeführt werden soll. Der optische Detektor besteht aus einem MLA, 
einer Septenmaske und einem CMOS-Sensor. Diese Komponenten ermöglichen einen 
besonders flachen Aufbau im Vergleich zur aktiven Fläche. Der Detektor soll zuerst in 
Einzehnessungen charakterisiert und anschließend mit mehreren baugleichen Detektoren 
auf einer Rotationsplatt form kombiniert werden. Zu diesem Zweck wurde ein spezielles 
Rahmenkonzept entwickelt und optimiert, das jeden Detektor als unabhängiges Modul 
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2. 



Material und Methoden 



kapselt . Ein einsatzbereiter Einzeldetektor i.sl in \bb. 2.1 dnrgesieili . 

Die Detektorkoinponenten und «leren Funktion werden ausführlich in Unterabschn. 2.1.1 
erörtert. Der Detektorrahmen, der die drei Komponenten genau gegeneinander ausge- 
richtet und gleichzeitig stabil zusammenhält wird in Unterabschn. 2.1.2 vorgestellt. In 
Unterabschn. 2.1.-5 wird die effektive Fokuslänge unter Berücksichtigung der MLA-Dicke 
festgelegt, die für die Berechnung tiefenkalibrierter Objektbilder mit dem Abbildungs- 
modell benötigt wird. 



110 mm 




Abb. 2.1. Einsatzbereiter optischer Detektor hl Auf- und Seitenansicht. Die aktive Fläche 
umfaast etwa 5 x 2,5 cm 2 bei einer Detektordicke von 8 mm. Der abbildbare Objektraum 
beginnt in unmittelbarer Nähe (wenige Millimeter) vor dem Detektor und erstreckt sich bis 
zu mehreren Zentimeter. Der Detektor ist im Hinblick auf einen kompakten Scanncraufbau 
für die Kleintierbildgebung konzipiert. Über den Montagewinkel wird der Detektor späier 
anf der Rotationsplattform (s. Absclm. 2.4) befestigt und über das Flachbandkabel mit der 
Ausleseelektronik verbunden. In der Aufsicht ist das MLA zu sehen, das die anderen beiden 
Komponenten verdeckt. Die Innenansicht des geöffneten Detektors ist In Abb. 2.5 zu sehen. 



2.1.1 Komponenten 

Die Abmessungen der Detektorkomponenten sind auf die Kleintierbildgebung abge- 
stimmt, d. h. der abzubildende Ob jektraum entspricht der Größe einer Labomiaus. Die 
aktive Detektorfläche beträgt etwa 5 X 2, 5 cm 2 , wohingegen die Dicke des zusammenge- 
setzten Detektors lediglich bei 8 mm liegt. Die für den Aufbau relevanten Spezifikationen 
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der drei Komponenten sind in Tal). 2.1 zusammengefa.sst und gegenübergestellt. 



Tab. 2.1. Die Bauteilgrößen Sind für jeden Detektor individuell vermessen worden (Iiier 
Detektor 4 aus dem Scanneraufbau, vgl. Ahl). 2.18). Der Rahmen wird anschließend passgenau 
gefertigt. Die aktive Fläche bezieht sich auf den Bereich der entsprechenden Elemente abzüglich 
der Bmiteilränder. Die Herstellung des ML As und der Septenmaske wurde auf den Sensor 
abgestimmt. Als Abstand zwischen den Mikroliusen bzw. zwischen den Bohrlöchern ist der 
zehnfache Sensorpixelahstand gewühlt worden. 



Bauteil 


MLA 


Septenmaske 


Sensor 


Abmessung (BxLxT) [mm 3 ] 


25.01 x 50.01 x 1 


28.46 x 53.51 x 2.15 


24.7 x 70.1 x 1 8 


aktive Flache (BxL) [mm 2 ] 


24.48 x 49 44 


25,44 x 50,4 


24 58 x 49 15 


Element 


Mikrolinse 


Bohrloch 


Pixel 


Anzahl (BxH) 


[px] 


51x103 


53x105 


512x1024 


Größe 


[mm] 


0 0,48 


0 0,4 


0,048 


Abstand 


[mm] 


0,48 


0,48 


0,048 



Das MLA in Abb. 2.2(a) besteht aus einer Vielzahl von plankonvexen Mikroliusen. 
die quadratisch nebeneinander in einer Ebene angeordnet sind, sodass alle optischen 
Achsen parallel verlaufen. Refraktive MLAs lassen sich mit herkömmlichen Halbleiter- 
Fertigungsprozessen mit unterschiedlichen Parametern wie Abmessungen. Linsendurch- 
messer, -abstand, -form und Fokuslänge herstellen [24]. Die Geometrie der Zentral- 
strahlen der parallel einfallenden Strahlenbündel, wie in Abb. 1.1 dargestellt, lässt 
sich in guter Näherung mit dem Lochkamorarnodell [25] für jede Mikrolinse be- 
schreiben. Das heißt, jeder Zentralstrahl durch das optische Zentrum repräsentiert 
rlie aus dieser Richtung parallel einfallenden Lichtstrahlen. Deren Intensitäl bündell 
sich im Schnittpunkt, des Zentralstrahls mit der Fokalebene (sprich im Fokuspunkt.) 
und wird dort vom Sensor abgetastet. Eine komplette Detektorabbildung kann man 
sich daher als gleichzeitige Aufnahme vieler kleiner Lochkameras vorstellen, die in einer 
Ebene nebeneinander positioniert sind. Entsprechend besteht das Ergebnis aus vielen 
kleinen Einzelprojektionen des Objektraumes, die sich zusammen einen gemeinsamen 
Bildsensor teilen. Aufgrund der Vielzahl der Lochkanuaas steht jeder Einzelnen nur noch 
ein kleiner Sensorbereich (10 x 10 px) zur Verfügung, in den auch nur ein entsprechend 
kleiner Teil des Objektraumes projiziert wird. Aus den vielen kleinen Einzelprojektionen 
lässt. sich mit dem Abbildungsmodell aus Absehn. 2.2 ein gemeinsames Objektbild 
berechnen, das mit dein einer klassischen Kamera vergleichbar ist. Auf diese Weise lässt 
sich durch das MLA ein komplettes Kameraobjektiv ersetzen und es ermöglicht zugleich 
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einen extrem Hachen Detektoraufbau. Ein weiterer Unterschied zu einer herkömmlichen 
Kamera ist der variable Objekt f'ukiis, der in der Sensorinformation selbst enthalten ist. 
Der Objektfokus einer klassischen Kamera wird über die Optik mechanisch vor der 
Bildaufnahme eingestellt, während er beim Detektor nach der Sensoraufnahme für die 
Objektbildberechnung beliebig gewählt werden kann (s. Abschn. 2.2.")). 




(a) (b) (c) 



Abb. 2.2. Komponenten des optischen Detektors unter einem Auflichtmikroskop (MZluF. 
Leica Microsystems AG. Heerbrugg, Schweiz): (a) MLA (26). (b) Septenmaske und (c) CMOS* 
Sensor [27] . Zu sehen sind jeweils Eckbereiche der drei Bauteile mit den Mikrolinsen, den 
Bohrlöchern und den aktiven Sensorpixeln. Die [landbreiten können mit entsprechender 
Vergrößerung exakt ermittelt und in der individuellen Rahmenkonstruktion berücksichtigt 
werden. 

Die Septenmaske in Abb. 2.2(b) hat die Aufgabe die einzelnen Mikrolinsen und deren 
Abbildungen auf dem Sensor optisch voneinander zu isolieren. Im Fall von Fokuspunkt- 
bildern kann man eine Mikrolinse auch als optischen Kanal betrachten, für den das 
Ubersprechen (engl, cross-talk) auf einen anderen Kanal verhindert werden soll [15]. 
Allgemein vermeidet die Septenmaske, dass die Einzelprojektionen der Mikrolinsen 
sich überlappen, d.h. ein Sensorpixel isi immer eindeutig einer einzigen Mikrolinse 
ZUOrdenbar. Mit einer Dicke, die der Äquivalent brennweite des MLAs entspricht, dient 
die Septenmaske zusätzlich als Abstandhalter zwischen MLA und Sensor. Die ideale 
Septenmaske sollte aus lichtundurchlässigem und minimal reflektierendem Material 
gefertigt sein. Die 1 hier verwendete Septenmaske besteht aus Aluminium. Die Reflexionen 
an der Aluminiumoberfläche in den Bohrlöchern sind in kurzen Abbildungstiefen in 
den Fokuspunktbildfirn bemerkbar und werden durch das Abbildungsmodell etwas 
abgeschwächt. Für einen zukünftigen Detektoraufbau sollton diese Reflexionen gerade 
mit Blick auf eine absolute Signalquantifizierung jedoch vermieden werden. Für eine 
Septenmaske aus Aluminium bietet sich z.B. an, die Oberfläche durch ein gängiges 
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EL ocal- Verfahren abzumatten. 

Bin großflächiger CMOS-Bilrlsensor, dargestellt in Abb. 2.2(e), tastet die durch die 
Septenmaske getrennten Einzelprojektionen der Mikrolinsen in deren Fokalebene ab. 
Der Sensor besteht aus einer Matrix von 512 x 1024 Photodioden und wird eigentlich 
für die Lichtdetektion hinter einer Szintillationsfolie in der digitalen Radiographie 
eingesetzt [28] . Die Abmessungen des Sensors waren Vorgabe für die Herstellung der 
anderen beiden Bauteile. So ist der Mikrolinsendurchmesser als ganzzahliges Viel- 
faches des Photodiodenabstandes gewählt, damit sieh die zentralen Fokuspunkte in 
regelmäßigen Pixelabständen ausrichten lassen. Die externe Ausleseelektronik [29], an 
die mehrere Sensoren über entsprechende Kabelverbindungen angeschlossen werden 
können, liefert Sensordatensätze mit . r >12 x 1000 px und 12-Bit-Grauworten. Die für 
die Signalmessung relevanten Spezifikationen des Sensors sind in Tab. A.2 im Anhang 
zusannnengefasst. 

Die Septenmaske deckt einen Teil der Sensorpixel ab, die folglieh nicht mehr für die 
Signalmessung zur Verfügung stehen. Eine mögliche Ausrichtung eines Bohrlochs und 
der entsprechenden Mikrolinse auf dem zugeordneten Sensorbereich ist in Abb. 2.3 
gezeigt. Die Bauteile werden in dieser Art über eine passgenaue Rahmenkonstruktion 
zusammengefügt, die im nächsten Abschnitt vorgestellt wird. 




(a) (b) (c) 



Abb. 2.3. Vergrößerte Ansicht des MI. As (a), der Bohrlöcher (l>) und der Sensoroberfläche 
(<•)• Für jedes Bohrloch bzw. jede Mikrolinse ist ein Sensorbcreich von 10 x lOpx vorgesehen, 
der durch die Septenmaske teilweise bedeckt ist. Die Lage des Bohrlochs bzw. der Mikrolinse 
innerhall) dieses Bereichs ist hier beispielhaft angedeutet. Es nniss nicht zwingend ein Pixel 
zentral im Bohrloch liegen, da für die zentralen Fokuspunkte im AbbildungsmodeU rationale 
OfFsets erlaubt sind. Die Abbildung ist entnommen aus [13]. 



Abkürzung für elektroly tische Oxidation von Aluminium: Eine Oberflächentechnik zur 
Erzeugung einer oxidischen Schutzschicht auf Aluminium durch anodische Oxidation. 
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Die exakte Ausrichtung der drei Bauteile gegeneinander ist entscheidend für die op- 
timale Bilderzeugung. Da die Mikrolinsen seihst mir schlecht zu erkennen sind, hilft eine 
Parallellichtquelle- senkrecht zum MLA. um die Fbkuspunkte auf der Sensoroberfläche 
sichtbar zu machen [5]. Die Position der Septenmaske lässt sich gleichzeitig über eine 
schwache, diffuse Hintergrundbeleuchtung finden. So ist die Lage des MLAs gegenüber 
der Septenmaske und des Sensors erkennbar. Einen Ausschnitt einer Sensoraufnahme 
eines kalibrierten Detektors unter diesen Bedingungen ist in Abb. 2.4 gezeigt. Nach 
korrekter Justierung müssen die zentralen Fokuspunkte nun in der Mitte der Bohrlöcher 
liegen und auf Sensorpixel treffen, die sich durch einen einheitlichen Offset in x- und 
y-Richtung beschreiben lassen. Dabei dürfen die Komponenten aber nicht gegeneinander 
verdreht sein. Für die zentralen Fokuspunkte bzw. die Schnittpunkte der optischen 
V : :v ■•in, -r Mik-Mhi;-,,. n:i> der Seiisorolwrlläcli ■ sind im Vlil.iMi:nr,M:->,d 'II üinndsätzlicli 
rationale üffsetkoordinaten erlaubt. Das heißt, die parallele Ausrichtung der Kompo- 
nenten und mittig in den Bohrlöchern auftreffende zentrale MLA-Fokuspunkte sind 
ausreichend für einen kalibrierten Detektor. 



Abb. 2.4. Fokuspunkte eines kalibrierten Detektors (Ausschnitt mit 12 x 9 Mikrolinsen): 
Parallel einfallendes Lieht, das senkrecht auf das MLA trifft, lässt die zentralen Fbkuspunkte 
auf der Sensoroberfläche erkennen, während eine gleichzeitige diffuse Ilintergrundbeleuchtung 
die Umrisse der Septenmaske sichtbar macht. Die zentralen Fokuspunkte müssen für einen 
kalibrierten Detektor in der Mitte der Bohrlöcher liegen. Die Helligkeitsunterschiede der 
einz« Inen B< ■hrlö< her sind auf Reflex* men in der Septenmaske zurüekzufühn H. 

2.1.2 Rahmenkonzept 

Der Detektorrahmen hat die Aufgabe, alle drei Komponenten justiert zusammenzu- 
halten ohne die Gesamtabmessung des Detektors unnötig zu vergrößern. Ein erstes 
Rahmenkonzept sah hier vor. die einzelnen Bauteile über Halteschrauben zu fixieren, 
damit, eine nachträglich Justage möglich bleibt (vgl. Abb. 3.1). Dieses Variante zeigte 

2 Paralleles Licht kann z.B. erzeugt werden, indem eine Glasfaser im Brennpunkt einer 
.Sammellinse positioniert wird. 
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sich jedoch bereits gegenüber geringen Erschütterungen wenig robust und war somit 
nicht für den späteren Scanneraufbau geeignet. Mit Verzicht auf die Justierschrauben 
ist das Konzept dann weiterentwickelt worden, indem die einzelnen Bauteile exakt 
vermes» n und für jeden DetektOl ein passgenauer Rahmen angefertigt wurde, sodass 
keine nachträglich Ausrichtung der Komponenten mehr erforderlich ist. 

Jeder der vier Detektoren im Scanneraufbau besteht aus zwei Rahmenkonstruktionen: 
einem oberen Teil, in dem das MLA und die Septenmaske gehalten werden, und einem 
unteren Teil, in dem der Sensor sitzt. Abb. "2.5 zeigt einen geöffneten Detektor, für 
den Auszüge der Konstruktionsskizzen in Abb. A.l und A.2 im Anhang zu finden 
sind. Beide Teile lassen sich über eine Zapfenverbindung zusammenfügen, an der die 
einzelnen Komponenten ausgerichtet sind. Dadurch können die beiden Rahmenteile 
unabhängig voneinander gefertigt werden. So lässt sich der untere Teil im Falle eines 
defekten Sensors einfach tauschen, ohne dass der obere Rahmen demontiert werden muss. 

Die Komponenten sind einzeln unter einem Auflichtmikroskop mit einem halben Sen- 
sorpixel Toleranz (±0.01 mm) vermessen worden. Die Randboreiche der Komponenten 
unter dein Mikroskop sind in Abb. 2.2 gezeigt. Die Sensoren unterscheiden sich nicht 
wesentlich in ihren Abmessungen, sodass die unteren Rahmenteile identisch konstruiert 
wurden. Die in dieser Arbeit eingesetzten MLAs und Septenmasken hingegen haben 
leicht variierende Randbreiten, die individuell in den Konstruktionsskizzen der oberen 
Rahmenteile berücksichtigt sind. Unter der erforderlichen Vergrößerung wurde dafür der 
Abstand jeder Bauteilkant» las zum Mittelpunkt der ersten Mi kr- linse bzw. di - erster] 
Bohrlochs gemessen. Ein zusammengesetzter Detektor benötigt keine weitere Justierung 
und ist weitaus stabiler in der Handhabung als das Vorgängermodell. Die aktive Fläche 
des Detektors ist die Schnittmenge der drei justierten Bauteile und entspricht in etwa 
der aktiven Fläche des Sensors. 

2.1.3 Effektive Fokuslänge 

Jede Mikrolinse wird später im Detektor-Abbildungsmodell auf ihr optisches Zentrum 
reduziert, durch das die aus unterschiedlichen Richtungen einfallenden Zentralstrahlen 
verlaufen. Die Fokuslänge / beschreibt im Abbilduugsmodell den Abstand dieser op- 
tischen Zentren zur Sensoroberfläche und hat direkten EinHuss auf die Abbildungstiefe, 
die einem Objektbild zugeordnet wird. Die Dicke des MLAs entspricht etwa dem zwei- 
fachen Linsendurchmesser und hat. damit eine nicht unerhebliche Ausdehnung. Die Frage 
ist nun. wo die optischen Zentren innerhalb des MLAs näherungsweise zu lokalisieren 
sind. Die für die optische Abbildung relevanten Parameter des MLAs sind in Tal). A.l 
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Abb. 2.5. Oberer und unterer Detektorrahinen in der Innenansicht. Im oberen Teil ist das 
MLA durch die Septeninnske verdeckt. Der CMOS-Sensor sitzt im unteren Teil und ist über ein 
herausgeführtes Flachbandkabel mit der Ausleseelektronik verbunden. Heide Teile lassen sich 
über die Zapfenverbindung zusammenfügen und sind dann exakt gegeneinander ausgerichtet. 
Der zusammengesetzte Detektor ist in Aldi. 2.1 zu sehen. 



im Anhang angegeben. Die konvexe Mikrolinsenseite zeigt im Experiment silaufbau in 
Richtung des Sensors. Die Orientierung des hier verwendeten MLAs ist für die Geometrie 
des Strahlenvorlaufs nicht relevant und führt lediglich zu Unterschieden im Abstand 
zwisc hen MLA-OberHäche und der Fokalebene. 



Abb. 2.6 zeigt, einen maßstabsget reuen, schematischen Querschnitt durch drei Mikro- 
linsen des optischen Detektors. In der oberen Mikrolinse ist der vereinfachte Verlauf der 
Zentralst iah hm durch ein optisches Zentrum angedeutet. Für dünne Linsen liegt dieser 
Funkt in der Linsenmitte. Für das ausgedehnte MLA werden die optischen Zentren 
näherungsweise zwischen die beiden Haupt ebenen gelegt. Diese Annahme lässt sich 
mit den Ergebnissen der Tiefenkalibrierung rechtfertigen, auf die im nächsten Absatz 
eingegangen wird. Die Lagen hi % ho der Hauptebenen berechnen sich nach [30] über 

fo(n-l)b 

hi ~ — R^r~ (21) 

und 



(2.2) 
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MLA Septenmaske Sensor 



Abb. 2.6. Maßstnbsgetreue Qiiersehnittskizze durch drei Mikrolinsen und die Septeninaske. 
Dir konvexe Seite clor Mikrolinsen zeigt zum Sensor. .Jede Mikrolinse wird auf ihr optisches 
Zentrum O reduziert, da.s näherungsweiso zwischen die beiden Hnuptcbenen {H\.: H>) gelegt 
ist. It\. /?2 bezeichnen die Krümmungsradien, V j . Vi die Scheitelpunkte, b die Dicke des 
MLA. /h die Äquivalentbrennweile und h\ den Abstand der ersten Ilauptebene zu deren 
Scheitelpunkt. Der Abstand h> ist aufgrund der Linseiigeoiuetrie Null, sodass H, und Vo 
aufeinander liegen. Durch die obere Mikrolinse ist der im Ahhildungsmodell vereinfachte 
Verlauf der Zentralstrahlen durch da.s optische Zentruni dargestellt. In der unteren Mikrolinse 
haben die drei optischen Zentren einen maßstäblichen Abstand von {2.7: 2.4: 2.1} min zum 
Sensor, die als Modellparamerer in Abb. 14 getestet sind. 

Mit den Krümmungsradien R\ = oc. R> = 1.81 mm, Linsendicke b = 1 mm. Bre- 
cliungsindox n = 1, SU und Äquivalentbronnweite /,> = 2.15 mm ergibt sich dann der 
Abstand h x 0.54 nun /wischen der ersten Hauptebene und deren Scheitelpunkt, 
wogegen die zweite Hauptebene mit ihrem Scheitelpunkt zusammenfällt. Der Mit- 
telpunkt zwischen beiden Haiiptebenen hat schließlich einen Abstand von 2,4 ±0.1-4 nun 
zur Sensoroberf lache, der nun als effektive Fokuslänge / (weiterhin aber einfach als 
Fokuslänge bezeichnet) des MLA verwendet wird. 

Die ermittelte Fokuslänge ist mit ihrer Fehlertoleranz in Abb. 2.0 auf allen drei optischen 
Achsen skizziert. Mit einer Linienquelle in bekannter Abbildungstiofe lässt sich über- 
prüfen, ob mit dieser Fokuslänge das fokussierte Linienbild auch tatsächlich wieder in 
der richtigen Abbildungstiefe berechnet wird. Die Ergebnisse der Tiefenkalibrierung sind 
in Absehn. 3.1. iL 2 gezeigt. Mit einer abweichenden Fokuslänge werden zwar dieselben 
Objektbilder berechnet, jedoch entstehen diese nicht in der korrekten Abbildungstiefe. 
Auf der optischen Achse der linieren Mikrolinse sind zwei zusätzliche Fokuslängen / 
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({2.1; 2,7} mm) außerhalb des Toleranzbereichs eingezeichnet. Werden <liese als Modell- 
parameter verwendet, dann weicht auch die Lage des fokussierten Linienbildes von der 
ursprünglichen Abbildungstiefe ab, wie in Abb. 3.14 zu sehen ist. Für möglichst genau 
tiefenkalibrierte Objektbilder, wie sie für die Oberflächensignalmessung erforderlich 
sind, ist somit die entsprechend festgelegte Lage der optischen Zentren wichtig. 

2.2 Detektor- Abbildungsmodell 

Das Detektor-Abbildungsmodell ist in dieser Arbeit entwickelt worden, um die Bild- 
gebung mit dem optischen Detektor zu verbessern. Ausgangslage war die Fokuspunkt- 
methode, die lediglich 1 % der verfügbaren Sensorinformation einer Detektorabbildung 
gleichzeitig zur Bilderzeugung nutzt. Abzüglich der durch die Septenmaske abgedeckten 
Bereiche beinhalten jedoch rund G4 V A der Sensorpixel Intensitäten uns dem Objekiraum. 
Wie in Abb. 1.2 dargestellt, lassen sich daraus z.B. mehrere Fokuspunktbilder erzeugen, 
die alle eine leicht unterschiedliche Perspektive auf das Objekt aufweisen. Über das Ab- 
bildungsmodeU hingegen werden all diese Informationen gleichzeitig in die Bilderzeugung 
eingebunden. Als Folge können Objektbilder berechnet werden, die sich in beliebigen 
Abbildungstiefen fokussieren lassen und dabei eine bessere räumliche Auflösung und 
ein höheres Signal-Rausch-Verhältnis aufweisen als ein Fokuspunktbild (gezeigt im 
Ergebnisteil Kap. 3). Das Detektor-Abbildungsmodell ist veröffentlicht in [13]. 

\'a< hfolgend behandeln die Unterabschn. 2.2.1 ff. den Mo lellramii so« ie lie V- n\ ärts 
und Rückwärtsabbildung. In den Unterabschn. 2.2.4 f. werden Eigenschaften der Sen- 
soroberfläche bzw. des Objektraumes erörtert. Weiter impliziert die tiefabhängige 
Objekt bildberechnung die Frage, ob die unbekannten Abbildungstiefen einer Intensitäts- 
wrteilung im Objektraum aus einer einzelnen Sensoraumahme bestimmt werden können 
(Unterabschn. 2.2.G). Mittels der Vorwärtsabbildung ist eine iterative Anwendung des 
Abbildungsmodi Iis möglich, die in Unterabschn. 2.2.7 vorgestellt wird. 

2.2.1 Systembeschreibung 

Die Modellelemente für das Abbildungsmodell sind in Abb. 2.7 skizziert. Das Detek- 
torkoordiiiatensystem (DKS) hat seinen Ursprung auf der Sensoroberfläche. Parallel 
dazu liegt in z-Rkhtung im Abstand der Fokuslänge / die MLA-Ebene und in einem 
Abstand z > 2/ beginnt das Objektvolumen, dessen Voxelgröße 3 an die Abstände a der 
Sensorpixel augepasst ist. Innerhalb der MLA-Ebene liegen im Abstand d die optischen 

3 Ein Voxel fvx] ist das dreidimensionale Äquivalent eines Pixels [px] 
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Abb. 2.7. Modellheschreibung: Die Positionen der Objektvoxel r, der Sensorpixel p und der 
optischen Zentren o sind im gleic hen Koordinatensystem definiert. / beschreibt die Fokuslänge 
einer Mikrulinse. a den Pixelabstand, d den Mikrolinsenabstand, D den Durchmesser eines 
Bohrlochs der Septenmaske und N, M, L sind die Anzahl der diskreten Elemente in der 
jeweiligen x-, y- und z-Richtung. Die Abbildung ist entnommen aus [13]. 

Zentren der Mikrolinsen. Der Bohrlochdurchmesser der Seplemnaske ist als Kreis in 
der Sensorebene angedeutet. 

Für die Projektion einer Mikrulinse wird die paraxiale Näherung [30. 31] in der geo- 
metrischen Optik verwendet, da die einfallenden Zentralstrahlen bedingt durch die 
Septenmaske lediglich in einem kleinen Winkel (< 5,4°) zur optischen Achse auftreffen. 
Jede Mikrolinse lässt sich dadurch als Lochkamera auf ihr optisches Zentrum reduzieren, 
durch das alle Zentralstrahlen aus den verschiedene Richtungen verlaufen. Diffraktion, 
Reflektion und Abhängigkeiten von der Wellenlänge des Lichtes werden hier vernach- 
lässigt. Die Dicke des MLAs wird über die Fokuslänge / (s. Absehn. 2.1.3) berücksichtigt. 

Die Nomenklatur der Pixel- und Voxelpositionen sowie der Intensitätszuordnung richtet 
sich nach [25], Eine diskrete Voxelposition im Objektraum ist beschrieben durch 
i*; f/.m,ii mit i = 0, . . . , (LAIN — 1) und der Voxelwert an dieser Stelle ist de- 
finiert als 9j := g(ti). Die Intensitätsverteilung im Objektvolumen ist somit durch 
den Vektor g r (jr('o)t • • • > 9(t{LMN-l))) gegeben. Analog wird eine diskrete Sensor- 
pixelpostion durch p, p m ,„ mit i 0 (M N — 1) bezeichnet und der Pixelwert 

als #f := g(pi). Die Intensitätswerte einer Sensoroberfläche sind dann im Vektor 
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g p - (ä(Po) y(P(A/.v-i))) zusammengelässt. 

Mit dem nachfolgenden Formalismus werden die Pixel-Voxel-Beziehungen über den 
Strahlenverlauf durch die optischen Zentren der Mikrolinsen bestimmt. Diese Bezie- 
hungen werden als Elemente einer dünn besetzten Abbildungsmatrix für eine feste 
AbbiIdi;i:-is:M'.- i : ihki l für die V< rwäits- m- : einmal füi di RtVkwnits-. ,hiMum! 
betrachtet. Dadurch könn n die berechneten Intens: ,n :i . i • il' . *i :i> • einheitlich ( >bj< kt- 
raumebene normiert, werden. Die Abbildungstiefe t bezeichnet den senkrechten Abstand 
von der Sensoroberfläche zu einem Voxel im Objektraum. 

2.2.2 Vorwärt sabbildung 




Abb. 2.8. Vorwärtsabbildiing: Intensitäten aus dem Objektranin werden auf die Sensorober- 
fläehe projiziert, Über die Pixelgrundfläche und ein optisches Zentrum wird ein spezieller 
Kau m hereich definiert, in dein jr nach Ahbildungstiefe / mehrere Objektvoxel liegen, die dem 
Sensorpixel dann zugeordnet werden. Die Abbildung ist entnommen ans [13]- 

Mit der Vorwärtsabbildung (Abb. 2.8) kann für einen vorgegebenen Intensitätsvektor 
aus dem Objektraum die Projektion auf der Sensoroberfläche berechnet werden. Jedes 
Sensorpixel integriert Intensitäten aus einem speziellen Rauinwinkel. der durch die 
Mikrolinsc und die Position des Pixels vorgegeben ist. Im abstrahierten Fall des Modells 
kann eine Beziehung zwischen einem Sensorpixel und allen Voxel u hergestellt werden, 
die innerhalb dieses Bereichs liegen, der hier die Form einer Pyramide mit quadratischer 
Grundfläche hat. Mit gegebener Abbildungstiefe t lässt sich eine Objektebene auswählen 
und der Schnittpunkt mit der Mittelachse der Pyramide berechnen: 
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Der Hochindex p bezeichnet hierbei den Bereich der Sensoroberfläche und der Hochindex o 
die MLA-Ebene. Weiterhin ist dieser Schnittpunkt nur dann gültig, wenn die Mitte- 
lachse nicht außerhalb der kreisförmigen Apertur des entsprechenden Bohrlochs der 
Septenmaake liegt. Das heißt, die Bedingung 



(x - n°df | (y - m*d)* < • | ) (2.4) 

muss erfüllt sein. In der Objektebene deckt der Querschnitt der Pyramide nun eine 
bestimmte Anzahl an Voxeln ab. Hieraus wird dann die Vöxelbeziehung € [0. 1| 
abgeleitet, die das Deckungsverhältnis eines Voxels durch die von Sensorpixel p< 
erzeugte Pyramide angibt. Alle Vöxelbeziehung eines bestimmten Sensorpixels können 
so als Zeile einer dünn besetzten Vorwärtsabbildungsmatrix V mit den Dimensionen 
(NM) X (NML) verstanden werden. Die Summe aller Voxelbeziehungen. die einem 
Sensorpixel zugeordnet sind, entspricht schließlich dem Vergrößerungsfaktor der Fläche 
des Sensorpixels: 

Die Abbildung wird weiterhin über einen einheitlichen Intensitätsvektor g r0 mit 
01° = IV/ normiert, der zu einer bestimmten Intensitätsverteilung g*' auf der Sen- 
soroberfläche führt. Die normierte Vorwärtsabbildung lautet 

g<" = U''Vg Hl . (2.6) 

Die Elemente der Diagonalmatrix U p = diag(uf) ergeben sich durch Umstellung dieser 
Gleichung zu 

"^'(^t)"- (27) 

Mit dem experimentellen Detektoraufbau wird später ein Normierungsbild g'*' mit einer 
planaren. homogenen Leuchtfolie erzeugt. Durch die Normierung ist die Gesanitsensi- 
tivität eines Sensorpixels berücksicht igt. die sich aus der eigentlichen Pixelsensitivität 
sowie der Lichttransmission durch die Mikrolinse und das Bohrloch der Septenmaske 
zusammensetzt. 
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Dil' normierte Intensit ät lässt sich für jedes Sensorpixel unabhängig von anderen Sen- 
sorpixeln berechnen und mit einer individuellen Abbildungstiefe t ausdrücken nls 



(2.8) 



Dieser Intensitätswert, entspricht dem Mittelwert über alle assoziierten Voxelintensitäten, 
gewichtet mit der Sensorsensitivität. 

2.2.3 Rückwärtsabbildung 




Abb. 2.9. Rückwärtsabbildung: Ein Objektvoxel ist gleichzeitig über verschiedene optische 
Zentren mit mehreren Sensorpixeln assoziiert. Basierend auf diesen Relationen werden die 
Pixelintensitäten in den Objektraum rückprojiziert. Die Abbildung ist entnommen aus [13]. 



Die Rückwärtsabbildung (Abb. 2.9) behandelt nun analog die Rückprojektion vor- 
gegebener Intensitäten auf der Sensoroberfläche in den Objektraum. Ein Voxel im 
Objektraum ist über mehrere Mikrolinsen mit verschiedenen Sensorpixeln verbunden. 
Die Voxelintensität summiert sich über diese Pixelintensitäten. Wie bei der Vorwärts- 
abbildung die Sensorpixel, so erzeugt jetzt ein Voxel mit einem Mikrolinsenzentrum 
eine quadratische Pyramide. Der Schnittpunkt, der Mittelachse mit der Sensoroberfläche 
lässt sich durch Invertieren von Gl. 2.3 berechnen: 
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Der Hochindex r bezeichnet hier das Objektvolumen. Wiederum ist der Schnittpunkt 
nur gültig, falls sich die Mittelachse innerhalb des Bohrlochs der betrachteten Mikrolin.se 
befindet: 

(x-n°d) 2 I {y-m°d) 3 < (j\ . (2.10) 

Die Projektion des Voxels r< auf die Sensoroberfläche ist um den Faktor [//(' — /)j 2 
verkleinert, verglichen mit der Grundfläche des Voxels. Das Dockungsverhältnis dieser 
Voxclprojektion mit einem beliebigen Sensorpixel pj wird ausgedrückt über die Pixel- 
beziehung urg 6 [0. 1). Alle Pixelbeziehungen, die einem bestimmten Voxel zugeordnet 
sind, lassen sich als Zeile einer dünn besetzen Rückwärtsabbildungsmatrix W mit den 
Dimensionen (NML) x (NM) schreiben. 

Der Gesamtbereich auf der Sensoroberfläche, der einem Voxel zugeordnet ist. lässt sich 
mit w,j angeben. Da diese Fläche mit der Position der Voxel variiert und direkt von 
der Form der Bohrlöcher abhängt, kann anders als bei Gl. 2.5 kein einfacher analytische 
Ausdruck für diese Summe angegeben werden. Neben den Sonsorsonsitivitäton u p - der 
verschiedenen assoziierten Pixel, wir:] nun noch Iii:' das Voxel selbsl ein Sensitivitäts 
faktor ii r , eingeführt. Hiermit wird erreicht, da.ss das Xormierungsbild g'" wieder auf 
die ursprüngliche Einheitsfläche g'" abgebildet wird: 

g r " = u r wuY' (2.11) 

Die Elemente der Diagonalmatrix LP' = diag(«. ) ergeben sich nach der Umformung zu 

h[ = — W . 2 . (2.12) 

Ist eine Intensitätsverteilung g'' auf der Sensoroberfläche gegeben und die Abbil- 
dungstiefe t bekannt, so lässt jede Voxelintensität unabhängig berechnen mit 

Ei mj [gy) j 

Obwohl die Vorwärts- und Rückwärtsabbildung in Matrixschreibweise formuliert sind, 
ist es nicht notwendig, für die Intensitätsabbildung die kompletten Matrizen aufzustellen. 
Vielmehr können Pixelwerte mit GL 2.8 und Voxelwerte mit Gl. 2.13 individuell berechnet 
werden. Mit diesem Vorgehen lassen sich auch Intensitätsabbildungen für eine große 
Anzahl an Sensorpixeln (hier 512 x 1000 px) und Objektraumvoxeln (z.B. CO x 512 x 
1000 vx) ohne großen Speicheraufwand realisieren. 
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2.2.4 Sensoroberfläche 

Ein Sensordatensatz enthält die zweidimensionale Intetisitätsverteiluug auf der Seusor- 
oborflächo und ist das Messorgehnis des optischen Detektors sowie die Grundlage für 
die Bildberechnimg. Für den Parametersatz eines optischen Detektors kann innner auch 
ein künstlicher Sensordatensatz über die Vorwärtsabbildung erzeugt werden. Damit 
ermöglichen Modellstudien einen direkten Vergleich mit der tatsächlichen Detektorab- 
bildung. Die Modellparameter des optischen Detektors sind in Tal). A.3 im Anhang 
zusammengef'asst. Abweichungen von den Standardwerten werden bei den jeweiligen 
Studien angegelten. Beispiele für Detektormessungen sind im Ergebnisteil in Abb. 3.2 
dargestellt. Zu sehen sind dort die einzelnen Mikrolinsenpro jektionen des Objektraumes 
auf der Sensoroberfläche, die durch die Septenmaske getrennt sind. 




(a) (b) (c) 

Abb. 2.10. Künstlich erzeugte Scnsordatcn (Ausschnitte mit je 150 x )50px): (a) Darstellung 
der Septenmaske und (b) Lage der zent ralen Fokuspunkte auf der Seitsoroberrläche. (c) Senso- 
rintensitäten nach der Vorwärtsabbildung dreier Piuiktquellen. die sieb in unterschiedlichen 
Tiefen im Objektraum befinden (gelb, magenta. eyan je {10: 20: .iojinm). 

Drei Beispiele für künst liche Sensordatensätze (200 x 200 px) sind in Abb. 2.10 gezeigt . 
Um die Lage der Septenmaske auf der Sensoroberlläche (a) sichtbar zu machen, sind 
alle durchlässigen Bereiche in der Vorwärtsabbildung auf einen Einheitswert gesetzt. 
Zugleich dient, solch ein Datensatz als Normierungsbild für die Rückwärtsabbildung 
eines künstlichen Sensordatensatzes, entsprechend analog zum Xormierungsbild des 
Detektors aus Abb. 3.2(a). In Abb. 2.10(b) sind die zentralen Fokuspunkte einheitlich 
markiert, die grundsätzlich auch in der Bohrlochuiitte der Septenmaske liegen müssen. 
Für ein Fokuspunktbild werden die Pixelintensitäten an fliesen Positionen aus einem 
kompletten Sensordatensalz extrahiert. Sichtbar werden die Fokuspunkte im optischen 
Detektor mit einer senkrecht zum Sensor ausgerichteten Parallelstrahlquelle. Die Dar- 
stellung der Fokuspunkte ( vgl. Abb. 2.4) dienen in erster Linie zur Kalibrierung der 
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Detektorkomponenten. Der Sensordatensatz in Abb. 2.10(c) ist die Vorwärtsabbildung 
(/ = 2mm) von drei Punktquellen im Objektraum. Jede Quelle besteht aus einem Voxel 
mit Einheitswert und sitzt in derselben y-Schicht (100 px) an äquidistanten x-Positionen 
({75: 100; 125} px) aber in unterschiedlichen Tiefen / ({10; 20; 30} mm). Dieses Beispiel 
stellt einen typischen Messdatensatz nach. Je größer die Abbildungstiefe ist. desto mehr 
Mikrolinsen und entsprechend mehr Sensorpixel sind an der MLA-Abbildung beteiligt. 
Umgekehrt steigt die Intensität in den verbleibenden Sensorpixeln, falls sich die Ab- 
bildungstiefe verkleinert. Sind im Objektraum nur wonige Punktquellen vorhanden, so 
lassen sich deren Signalmuster auf der Sensoroberfläche durchaus unterscheiden, wie in 
diesem Beispiel zu sehen ist. Zwar greifen die drei Muster ineinander, belegen jedoch größ- 
tenteils getrennte Sensorpixel. Auf diese Weise erzeugt jede Punktquelle im Objektraum 
ein eindeutiges Signalmuster auf einer beliebig fein abtastbaren Sensoroberfläche, 

2.2.5 Objektraum 

Der Objektraum besteht aus einem diskreten Volumen (Objektdatensatz), das im 
Weiteren entsprechend schichtweise dargestellt wird und für das folgende Konvention 
gelten soll: 

• Objoktbild/z-Schicht: x-y-F.bene. adressiert über Abbildungstiefe f in [mm] 

• x-Schicht: y-z-Ebene, adressiert über x-Position in [px| 

• y-Schicht: x-z-Ebene, adressiert über y-Position in [px| 

Im Objektraum sollen nun die Voxelinteiisitäten nach einer Detektorniessung wieder 
an ihrem Ursprungsort rekonstruiert werden. Sind die Abbildungstiefen unbekannt, 
so lassen sich über schrittweise variierende Abbildungstiefen sukzessive Objektbilder 
berechnen und zu einem Objektdatensatz zusammenfassen. Für jedes Voxel in diesen 
Abbildungstiefen wird also eine Intensität bezüglich seiner assoziierten Sensorpixel be- 
rechnet unabhängig davon, ob bei der Sensoraufnahme an dieser Stelle im Objektraum 
auch tatsäi hl ich eine I J< litquelle v< rhanden war. In einem Objekt bild zeigt sich dann das 
Objekt scharf bzw. wird fokussiert abgebildet, falls die Abbildungstiefe auch tatsächlich 
dem Abstand des Objektes zur SensoroberHäche entspricht . Ob jekt e oder Ob jekt bereiche, 
die sich ursprünglich nicht in dieser Abbildungstiefe befanden, erscheinen unscharf bzw. 
werden delb kassiert abgebildet. Aus einem einzigen Sensordatensatz lässt sich somit ein 
Objektdatensatz mit beliebigen Abbildungstiefen erzeugen, der anschließend „durchfo- 
kussiert" werden kann. Dieses Phänomen wird in der Literatur als digitales Refokussieren 
(engl, digital refocusing) bezeichnet [17. 18]. Als Beispiel hierfür ist der Objektdaten- 
satzes eines dreidimensionalen Mausphantoms in Abschn. 3.1.1.1 im Ergebnisteil gezeigt. 
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Voraussetzung für die fokussierte Rüekwärtsabbildung eines Sensordatonsatzes in den 
Objekt räum ist das Wissen um die genauen Abbildungstiefen. Für künstliche Sensorda- 
tensätze. berechnet aus der Vorwärtsabbildung, sind die.se als Modellparameter festgelegt 
und somit auch für die Rückwärtsabbildung eindeutig bekannt. Anders verhält es sich 
bei der Detektormessung, für die in der Regel die Abbildungstiefen unbekannt sind. 
Grundsätzlich müssen die Abbildungstiefen nicht zwingend konstant in einer z-Schicht 
liegen. Vielmehr ist ein entscheidender Vorteil des Abbildungsmodells, dass Voxelin- 
tensitäten, die z.B. auf einer dreidimensionalen Oberfläche lokalisiert sind, mit jeweils 
unterschiedlichen Abbildungstiefen aus den Sensordaten berechnet werden können. 
Somit wird der Objektfokus nicht in eine Ebene, sondern auf diese Oberfläche gelegt . 
Auf diese Weise wird in Abschn. 3.2 im Ergebnisteil die Signalintensität direkt auf der 
zuvor bestimmten Zylinderoberflfiche berechnet. 

Die Rückwärtsabbildung einzelner Sensorpixel lässt sich anhand des künstlichen Sen- 
sordatensatzes mit den Fokuspunkten aus Abb. 2.1()(b) zeigen. Der daraus berechnete 
Objektraum (/ = 2mm) ist in Abb. 2.11 dargestellt. Die Abbildungstiefe / ist mit 
einem Inkrement von At = 0,15 mm je Objektbild zwischen 5 mm und 35 mm variiert. 
In einer y-Schicht sind die sich aufweitenden Bereiche im Objektraum zu sehen, in 
lie jedes Fokuspunktpixel rückabgebildel wird. Mit ansteigende! Abbildungstiefe ver 
größert sich der entsprechende Objektraumbereich eines Sensorpixels innerhalb einer 
z-Sehieht, auf den sich dann die gesamte Pixelintensität verteilt. Die daraus folgende 
Intensitätsabnahme ist im Profildiagramm (durchgezogene Linie) durch eines der Fo- 
kuspunktpixel dargestellt . Der stufenförmige Intensitätsverlauf hat seine Ursache in der 
Diskretisierung des Ob jektraumes. Im gleichen Diagramm ist zum Vergleich der Ideal- 
verlauf der Intensitätsabnahme einer Punktquelle nach dem Abstandsquadrat-Clesetz 
dargestellt (gestrichelte Linie). An den Randbereichen des Objektraumes kommt es 
zu einer uneinheitlichen Voxelgewichtung, da liier teilweise Sensorpixel abgeschnitten 
werden, die dann für die Mittelwertbildung fehlen. Die Randbereiche des Objektrauines 
sollten deshalb für die Signalmessung möglichst ausgespart werden. Zwar werden die 
Randarlefakte mit ansteigender Abbildungstiefe schwächer, wie in den z-Schichten in 
drei verschiedenen Abbildungsticfen t ({10; 20: 30} mm) und deren Profildiagrammen zu 
sehen ist. jedoch vergröbert sich gleichzeitig auch die Anzahl der betroffenen Randvoxol. 
Ab einer bestimmten Abbildungstiefe beginnen sich die Objektraumbereiche der Fokus- 
punktpixel zu überlappen. 

Der Sensordatensatz der drei Punktquellen von Abb. 2.1()(c) ist nun vergleichbar mit 
einer realen Punktquellenmessung des Detektors. Der Objektraum ist wie zuvor unter 
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Abb. 2.11. Objekt räum nach der Rückwärtsabbildung der zentralen Fokuspunkte aus 
Abb. 2.10(b). Anhand dieses Beispiels soll die Verknüpfung einzelner unabhängiger Sen- 
SOfpixel (liier Fukuspunktpixel) mit dem Objektraum dargestellt werden. Es ist zu erkennen, 
dass mit ansteigende Abbildungstiefe der EinHuss einzelner Seusorpixel auf die Objektvoxel 
immer schwächer wird. Zudem komml es an den Randbereichen des Objektraumes zu einer 
ungleichen Voxelgewichtung. Die Voxelintensitäten sind in den einzelnen Schichten farbkodiert 
und entsprechen der Skala der y-Schicht. 

denselben Bedingungen berechnet und in Abb. 2.12 gezeigt. In der y-Schicht und den ne- 
benstehenden Pronldiagraminen isl die rückabgebildete Inteiisitätsverteilung im Objekt- 
datensatz zu sehen. Zu beobachten ist. dass sich an den tatsächlichen Quellenpositionen 
die Voxelwerte zu lokalen Maxima aufsummieren. Mit ansteigender Abbildungstiefe 
verringert sich dieser Extremwert jedoch, da sich die ursprüngliche Quellenintensität, 
die bei allen drei Quellen gleich war. auf einen immer größeren Objektraumbereich 
verteilt. Gleiches ist in den drei Objektbildern im Abstund der Quellenpositionen zu 
sehen, die nun jeweils die tiefenabhängige Punktbildfunktion des Abbildungsmodells 
zeigen. Weiterhin sind in jedem dieser Objekt büder ttuch Intensitäten der beiden anderen 
(defokussierten) Quellen zu linden. Durch die Modellnonnierung isl die Summe aller 
Voxelintensitäten in jedem Objektbild gleich, doch abhängig von der Abbildungstiefe 
räumlich unterschiedlich verteilt. In den drei x-Schichten und den Prolildiagranimen 



28 



2. Material und Methoden 




Abb. 2.12. Rückw&rtsabbildung dos Sensordatensatzes aus Abb. 2.10(c). Die drei Punkt- 
quellen liegen auf derselben y- Koordinate (100 px) in unterschiedlichen Tiefen ({10: 20: 30} mm) 
und x- Abständen ({7Tv, 100: 125} px). In der y-Schicht sind die ursprünglichen Tiefen der 
Punktauellea zu sehen, in denen sich die nickabgehildeten Intensitäten treffen. Die Intensit.äts- 
verteilung entlang der Ahbildungstiefe ist für jede Punktquelle uls Profildiagrarnm dargestellt. 
In den z-Schichten ist die tiefenabhängige Verbreiterung der Punktbilder zu erkennen und 
ebenfalls als Profildiagramm dargestellt. In den jeweiligen K-Schichten können die Quellen 
über eine schichtweise Skalierung separat dargestellt und über ihren Maximalwert einer Abbil- 
dungS tiefe zugeordnet werden. In der mittleren x-Schicht sind im unteren Bereich Anteile einer 
fremden Punktquelle vorhanden, die ebenfalls im (Mitsprechenden Linienprofi] aus der y-Schicht 
zu sehen sind. Die Voxelintensitäten sind farbkodiert und entsprechen dem Skalenbereich 
der y-Schicht. Die Darstellungen der drei x-Schichten sind innerhalb dieses Bereiches jeweils 
separat gewichtet. 

können die Quellenpositionen durch eine separate Skalierung über ihren Maximalwert 
gut identifiziert und so einer Abbildungstiefe zugeordnet werden. 
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2.2.6 Tiefenbestimmung aus einem Sensordatensatz 

Die Modellstudie mit den drei Punktquellen hat gezeigt, dass die Tiefeninformation 
in einem einzelnen Sonsordatonsatz durchaus vorhanden ist und sich in der tiefenah- 
hängigen Intensitätsvertoilung des Objektdatensatzes äußert. Deshalb sind in dieser 
Arbeit zwei Ansätze zur Analyse des Objektdatensatzes verfolgt worden, um die Abbil- 
dungstiefen aus einer einzigen Detektorabbildung bestimmen zu können. Das Prinzip 
dabei ist. die fokussierten Bereiche über Intensität smessgrößen zu identifizieren und 
den entsprechenden Abbildungstiefen zuzuordnen. Im Ergebnisteil in Absehn. 3.1.1.2 
wird dann der Frage nachgegangen, ob diese Ansätze für die Anwendung auf eine 
dreidimensionale Oberfläche praktikabel sind. 

Ein Ansatz für die Tiefenbestimmung ergibt sich direkt aus der Modellstudie in Abi». 2.12. 
Für l'unktquellen lassen sich die ursprünglichen Koordinaten über die Suche nach lokalen 
Voxelmaxhna finden. Als Limitierung dieses Ansatzes ist aber schon in der Studie zu 

■ •rkonti) ii. lass sich ein Extren w erl ii r schlecht i limln i.'isst . je größer die Abbil 

dungstiefe der Punktquelle ist. Zudem ist ein eindeutiges Voxehnaximum schwieriger 
zu isolieren, je dichter die Punktquellen angeordnet sind, wie es z. B. bei einer zu- 
sanimenhängendeii Oberfläche ler Fall ist. Für ein beleuchtetes Mausphantom ist der 
Objektdatensatz in Abb. 3.1 im Ergebnisteil hinsichtlich dieser Messgröße analysiert 
worden. 

Ein anderer Ansatz nutzt die Strukturinformation eines Objektbildes, für den als 
Messgröße die Varianz der Voxel intensitäten innerhalb eines gewissen Objektbildbereiehs 
verwendet wird. Unter Annahme, dass die Varianz maximal wird, falls die Oberflachen- 
struktur in diesem Bereich im Objektfokus liegt, lassen sich die entsprechenden Abbil- 
dungstiefen finden. Dieses Verlähren wurde ausführlich für die Oberflächenextraktion 
aus digitalen Hologrammen [32, 33, 34] untersucht. Diese Hologramme zeigen jedoch 
gegenüber dem Objektdatensatz eine weitaus höhere räumliche Auflösung. Dennoch 
lässt sich die Methode für den Objektdatensatz adaptieren und wird folgendermaßen 
durchgeführt: Um jedes Voxel in einem Objektbild wird eine quadratische Umgebung 
(Block) gelegt. Anschließend wird in dieser Umgebung die Varianz der Voxelinten- 
sitäten berechnet und dem zentralen Voxel zugeordnet. Uber diese Varianzwert* 4 wird 
dann entlang der z-Aehse über alle Objektbilder der Maximalwert ermittelt und einer 
Abbildungstiefe zugewiesen. Diese lassen sich dann in einer Tiefenkarte zusammen- 
fassen. Grobe Sprünge in einer Tiefenkarte können danach mit einem kantenerhaltenden 
Metlianfilter geglättet werden. Grundsätzlich erfordert diese Methode immer einen 
Koinpromiss zwischen der Blockgröße und der Feinheit der Tiefen karte. Das Verfahren 
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ist ebenfalls auf den Objektdatensatz eines Mausphanlums angewendet wurden und die 
resultierenden Tiefenkarten im Ergebnisteil in Abb. 3.5 gezeigt. 

Beide Ansätze sind letztlich über die räumliche Auflösung des optischen Detektors 
limitiert, die zusätzlich noch mit der Abbildungstiefe abnimmt. Hinzu kommt, dass die 
Tiefenauflösung weitaus geringer ist als die räumliche Auflösung in einem Objektbild. 
Soviel sei der Diskussion in Kap. 4 vorweggenommen, dass obwohl die Tiefeninformation 
bereits in einem ein/einen Sensordatensatz vorhanden ist. die Oberflächenextraktion 
aus mehreren Sensoraufnahmen in Absehn. 2.4.3 deutlich bessere Ergebnisse erzielt. 

2.2.7 Iterative Anwendung des Abbildungsmodells 

Die Vorwärtsabbildung lässt sich neben Modellstudien ebenso für eine iterative An- 
wendung des Abbildungsmodells einsetzen. Für eine feste Abbildungstiefe (d.h. ein 
bestimmtes Objektbild) kann so eine Bildentfaltung (engl, deconvolution) durchgeführt 
werden, die den Einfluss des abbildenden Systems auf das Objektbild verringert. Eine 
gängige Methode ist die iterative Objektschätzung nach Van Gittert [35] mit additive 
Korrokturfaktoreu, die hier kurz vorgestellt und auf das Abbildungsmodell übertragen 
wird. 

Im eindimensionalen Fall sind die Bilddaten i(x) = s(x)0o(x) eines Objektes o(x) durch 
die Faltung mit einer Verbreiterungsfunktion s(x) verändert. Werden die Bilddaten i(x) 
als erste Approximation ö (ü) (:r) zum Objekt o(.r) betrachtet, dann führt ein erneuter 
Systemdurchgang von ö (a) (x) ebenfalls zu einer guten Näherung ? (0) (x) dieses Objektes. 
Die Methode von Van Gittert beruht nun darauf, dass die einfach zu berechnende 
Differenz der Bilddaten i(x) - j (0) (j-) in Relation zum Fehler in der Objektschätzung 
o(x) - ö (0) (ar) steht. Der Fehler in den Büdaten kann demnach als Korrekturfaktor 
zur iterativen Objektschätzung verwendet werden. In allgemeiner Form für ö w (x) und 
6 ( * +t) (x) als k-te und (k+l)-te Approximation des Objektes lautet dann die Schätzung: 

ö (k+l) (x) = ö {k) (x) + \i(x) - s(x) gl ö w (x)}. (2.14) 

In Anlehnung an Van Citterts Methode ist in Abb. 2.13 die iterative Anwendung des 
Abbilduugsmodells [36] auf einen Sensordatensatz gezeigt. Vom k-ten Objektbild wird 
mittels Vorwärtsabbildung der k-te Sensordatensatz berechnet und vom ursprünglichen 
Sensordatensatz subtrahiert. Die Differenz wird über die Rückwärtsabbildung in den 
Objektraum transformiert, dem k-ten Objektbild hinzuaddiert und so das (k f l)-te 
Objektbild erzeugt. Das Objektbild g r ' 0 ' nach der einfachen Rückwärtsabbildung ent- 
spricht hierbei der Standardanwendung des Abbildungsmodells, g'' kann entweder ein 
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künstlicher oder mit dem Detektor gemessener Sensordalensatz sein. 




Abb. 2.13. Iterative Anwendung des Abbildungsiuodells auf einen Sensordatensatz g''. Die 
Differenzdaten g'' — g f>|t> werden in den Objektraiini transformiert und zum Objektbild g'"' 1 
addiert. Dadurch verbessert sich der Kontrast im neu entstandenen Objekt bikl g , (fc * 1 ). 



Die Auswirkung der iterativen Anwendung des Abbildungsmudells auf ein Objektbild 
(t = 14 mm) ist in Abb. 2.14 demonstriert. Als Eingangsbild dient ein Rechteckmuster 
(1,7 'p/nun) der Testschablone (Abschn. 2.3.1.1). das als Profildiagramm dargestellt ist 
(a). Das mittlere Diagramm (b) zeigt das Ausgangsbild nach der iterativen Anwendung 
in einer Modellstudie. Mit jedem zusätzlichen Durchlauf verbessert sich das Verhältnis 
zwischen Minima und Maxima. Ohne Ahbruchkriterium kommt es jedoch irgendwann 
zu Uberschwingern an den Kanten, die bei den künstlichen Sensordaten zu negativen 
Werten führen können. Gleicherweise ist mit einem gemessenen Sensordatensatz (c) 
verfahren worden. Auch hier lässt sich eine Kontrastverbesserung beobachten, wenngleich 
mich hi hwä< Im als in d> •! Moi:» -Ibridic. I )')•• iterative Anw* -udum; s Ahliil:iu:..d<m •II- 
wurde in Abschn. 3.1 mit der Testschablone auf die Frage untersucht, ob diese Methode 
Vorteile für den späteren Einsatz im Scanneraufbau bringt. 



2.3 Charakterisierung der Abbildungseigenschaften 

Für den späteren Einsatz im Scanneraufbau sollen nun die relevanten Abbihlungsei- 
genschaften des optischen Detektors näher betrachtet werden. Neben der eigentlichen 
Sensoraufnahnie spielt die anschließende Bilderzeugung aus den Sensordaten dabei eine 
wesentliche Rolle. Nachfolgend werden Messmcthoden vorgestellt, mit denen sich die 
räumliche Auflösung und das Signal-Rausch- Verhältnis bestimmen und für Fokuspunkt- 
bilder und das Abbildungsmodell vergleichen lassen. 

Die relevanten Abbildungstiefen im Scannerobjektraum betragen zwischen 15 mm und 
50 mm. Für die Oberflächenextrakt um kann die Ob jektausleuchtung an die Sensor- 
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Abb. 2.14. Auswirkung der iterativen Anwendung des Abbildungsinodells auf ein Objekt bild: 
(a) Profildiagramm durch ei» Rechteckmuster (1,7'p/mm) der Testschablone, die als Vorlage für 
eine Modellstudie und eine Detektonnessuiig in einer Abbildungstiofc von t - 11 nun verwendet 
wurde, (b) Resultierender Kontrast nach einfacher Rückwärtsabbildung (0) in der Modellstudie 
und Kontrastverbesseriing während der iterativen Anwendung des Abbildungsinodolls (!,..,">). 
(c) Die Kontrastverbesseriing für die Detektormessung ist erwartungsgemäß schwächer als in 
I' r \l: ilfllsi ii' Ii. . 

emphndliolikoit angopasst worden. bei der eigentlichen OberflächensignalmosRung sind 
allerdings nur sehr schwache Quellensignale zu erwarten. Die A I )l>il< l ungseigensehaf'ten 
wurden mit Blick auf diese Voraussetzungen untersucht. 



2.3.1 Räumliche Auflösung 

Der optische Detektor kann für eine feste Abbildungstiefe / als lineares, verschiebungsin- 
variantes System (LTI-Systom) betrachtet werden, dessen Antwort auf ein Eingangssignal 
entweder ein Fokuspunktbild oder ein Objektbild des Abbildungsmodells ist. Den größten 
F.influss auf die räumliche Auflösung haben die Abbildungsgoometrie des Ml. As und die 
Abtastung der Fokalebene durch den Sensor. In einem Fokuspimkt bild stellt jode Mi- 
krolinse eine Detektoreinheit dar. deren Abstand somit gemäß «lern Abtast .theorein [37] 
die (Jrenzfrequenz f y = '/im = 1.0 1 '''/um, vorgibt. Für das Abbildungsmodells hingegen, 
mit einem Sensorpixel als Detektoreinheit, gilt eine zehnfach höhere CJrenzfrcquenz 
/,, = [/■>« = 10. 1 'i'/mm. Durch die Abbildungsgoometrie dos MLAs nimmt die erreichbare 
räumliche Auflösung mit ansteigender Abbildungstiofc ab. da ein immer größer wer- 
dender Objektraumbereich in ein Sonsorpixol pro jiziert wird. 



Mit der planaren Testscliablone ans Absclin. 2.3.1.1 lässt sieh sofort die Kontrast - 
schwäehung in den Objektbildern bestimmen und so auf die räumliche Auflösung 
schließen. Die Recht eekmust er eignen sich durch die gewählte Anzahl der Linienpaare 
besonders für kleine Abbildungstiefen (8mm < i < 25 mm) und erlauben zudem einen 
direkten Vergleich mit einem Fokuspunkt bild. 
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Eingehender und auch für größere Abbildungstiefen (15mm < t < 70 mm) lässt sich 
die räumliche Auflösung mittels einer Linienquelle untersuchen (Abschn. 2.3.1.2). Ein 
LTI-Systeni ist durch die Inipulsantwort vollständig charakterisiert, die sich über die 
Abbildung einer Punktquelle oder alternativ mit einer Linienquelle ermitteln lässt. Die 
Verbreiterung des Linienbildes kann in unterschiedlichen Abbildungstiefen bestimmt 
und so ein Zusammenhang zur räumlichen Auflösung hergestellt werden. Zusätzlich 
lässt sich damit eine tiefenabhängige Modulationsübertragungsfunktion (MTF. engl. 
modulation transfer function) für das Abbildungsmodell ableiten. Die MTF beschreibt 
die Abschwächung einzelner Signalfrequenzen (hier Raumfrequenzen) durch ein abbil- 
dendes System. Außerdem hat sich gezeigt, dass die Liuienbildfimktion effizient für eine 
Tiefenkalibrierung des Abbildungsmodells verwendet werden kann. 

2.3.1.1 Testschablone 
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Abb. 2.15. Planare Testschablone für die Kontrastübertragung des optischen Detektors. 
Die Linienbreiten variieren in der digitalen Vorlage von 1 px bis 10 px und erzeugen somit 
im Folienausdnick für die Detektormessimg Linienpaare von 10. 4 'p/mm bis 1,0'p/mm. Alle 
Strukturen haben einen einheitlichen Kontrast. Die Testschabloue ist in Anlehnung an das 
USAF-1951-Testinuster [38] entworfen. 

Die Testschablone in Abb. 2.1ö besteht aus verschiedenen Linienanordnungen mit ein- 
heitlichem Kontrast, die durch die Detektorabbildung entsprechend verändert und in 
den resultierenden Objektbildern beurteilt werden können. Die Schablone ist speziell 
für den Vergleich der Detektormessung mit einer Modellstudie entworfen worden und 
lässt sich gleichzeitig als Mess- und digitale Schablone verwenden. In beiden Fällen 
lassen sich so Sensordatensätze in bekannten Abbildungstiefen erzeugen und daraus die 
fokussierten Objektbilder berechnen. Das Auflösungsvermögen kann dann für jedes Li- 
nienpaar per Augenmali oder über das resultierende Verhältnis des Ausgangskoni rast es 
zum Eingangskontrast bestimmt werden. Außerdem ist mit dieser Methode ein direkter 
Vergleich mit den Fokuspunktbildem möglich, die aus denselben Sensordaten extrahiert 
werden. 
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Als Strukturen sind zehn Linienpaare in horizontaler und vertikaler Orientierung neben- 
einander angeorduet. Die Linienbreite variiert zwischen 1 px und 10 px in der digitalen 
Schablone und entsprechend zwischen lO. l'i'/mm und l.O'i'/mm auf clor Messschablone. 
Die Testschablone ist für die Detektonnessung mit einer Pixelbreite von 0,048 mm auf 
durchsichtiger Folie gedruckt und vor einer Leuchtfolie angebracht. Die Linienbreiten 
wurden so gewählt, dass nur das gröbste Muster gerade noch von einem Fokuspunkt bild 
aufgelöst werden kann. Die schmaleren Linienpaare hingegen sind abhängig von der 
Abbildungstiefe nur noch mit dem Abbildungsmodell rekonstruierbar. Aufgrund dieser 
feinen Strukturen eignet sich das Testmuster speziell für die Auflösungsmessung in 
kleineren Abbildungstiefen {t < 25 mm). Die Resultate der einfachen Rückwärtsab- 
bildung und der Effekt auf die Kontrastübertragung nach iterativer Anwendung des 
Abbilduugsmodells sind im Ergebnisteil in Abschn. 3.1.2 gezeigt. 

2.3.1.2 Linienbildfunktion 

Aus Sicht der Systemtheorie gilt für ein LTI-System, dass die Systemcharakteristik 
vollständig durch die Impulsantwort oder Punktbildfunktion (PSF. engl, point spread 
function) beschrieben ist [37. 39]. Die tiefenabhängige Punktbildfunktion h t (x,y) des 
optischen Detektors lässt sich messen, indem eine Punktquelle in einer bekannten 
Abbildungstiefe aufgenommen und anschließend über die Rückwärtsabbildung das 
entsprechende Objektbild berechnet wird. Nach der Fouriertransformation erhält man 
daraus die komplexe Übertragungsfunktion H t (u,v). Die tiefenabhängige, normierte 
MTF ist dann definiert als 



MTF,(u. r) = 



H t (u,v) 



(2.13) 



//i(0,0) 

Für gewöhnlich erweist sich jedoch die Messung der Linienbildfunktion ( LSF. engl, line 
spread function) übei einen schmalen Schlitz als prakl is.. In r g- -gi -nübi r d> -r Mi ssung d> t 
PSF. Linienquellen werden häufig in der digitalen Radiographie für die Messung der 
MTF verwendet [40]. Es lässt sich zeigen, dass man zur komplexen 1 bertragungsfunktion 
ebenso durch eine eindimensionale Fouriertransformation der LSF gelangen kann [41]. 
Unter Annahme einer isotropen Abbildung ist die Orientierung der Linienquelle sogar 
unerheblich. 



Als Linienquellen dienen für die Detektormessung zwei verschiedene Schlitzschablonen 
(Schlitzbreite je 25 /im bzw. 100 /im), die vor einer Leuchtfolie angebracht sind. Für die 
Modellstudie wird eine digitale Linienquelle (Linienbreite lpx) verwendet. In allen drei 
Fällen ist für genügend große Abstände (f > 8 mm) gewährleistet, dass die Ausdehnung 
der jeweiligen Linienquelle klein gegenüber deren Bildfunktion nach der Rückwärtsab- 
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bildung ist, Der Vorteil dieser Methude liegt auch hier wieder darin, dass die Liiiien- 
bildfunktion gleichzeitig aus den gemessenen und künstlichen Sensordaten für dieselben 
Abbildungstiefen bestimmt und verglichen werden kann. Mit einem Fokuspunktbild lässt 
sich die Linienbildfunktion aufgrund der zu geringen Abtastung nicht mehr adäquat 
ermitteln. Abhilfe schafft hier z.B. . den Schlitz gegenüber dem Detektor leicht zu 
verdrehen und durch geschickte Kombination verschiedener Pixelwerte die LSF höher 
abzutasten, aus der sich dann die MTF berechnen lässt [40]. Allerdings sei erwähnt, 
dass es sich hierbei tatsächlich um die MTF vor der Digitalisierung (engl, presampled 
MTF) handelt, die eben nicht durch ein Fokuspunktbild aliasingfrei übertragen werden 
kann. 

Aufgrund der Abbildungsgeometrie des MLAs wird die räumliche Auflösung mit anstei- 
gender Abbildungstiefe immer schlechter. Die Frage ist, wie sieh das Auflösungsvermögen 
für den Tiefenbereich des Scannerobjektraumes verhält. Für die Charakterisierung der 
räumlichen Auflösung ist primär dieser tiefenabhängige Zusammenhang interessant. 
Als relevante Messgröße eignet sich hierfür die Verbreiterung der Linienquelle in einem 
Objektbild, die reziprok mit der räumlichen Auflösung in Verbindung steht, wie im 
Folgenden ausgeführt wird. Die Messung verläuft nach folgendem Prinzip: Von einer 
Linienquelle wird in bekannter Abbildungstiefe ein Sensordatensatz aufgenommen, 
anschließend mittels Rückwärtsabbildung ein Objektbild berechnet und dann das Lini- 
enbildprofil durch eine Gaussfunktion approximiert mit 



Die Gaussfunktion wird über eine nichtlineare Regression (Nelder-Mead-Simplex) an 
die Messdaten angepasst, Die notwendige Funktionalität, dafür wird von einer Java- 
Bibliuthek [42] bereit gestellt, die für die Datenanalyse eingesetzt wurde. Als Sehätzpa- 
rameter orgeben sich ein Offset y 0 . der Funktionsmittelpunkt x 0 , eine Amplitude c und die 
Funkti .nsbroite <j. Wie gui 'Ii" geschätzte Gaussfunktion die Messdatei] besehreibt, lässi 
sich über das Bestimmtheitsmaß R 2 € [0, 1] ausdrücken. Der Wert entspricht dem qua- 
dratischen Korrelationskoeffizient und nähert sich eins, je besser Mess- und Schätzdaten 
übereinstimmen. Die Halbwertsbreite (FWHM, engl, füll width at half maximum) 
einer Gaussfunktion lässt sich über den Zusammenhang FWHM 2v'2 In 2a ss 2.35(7 
ausdnicken. Durch den konstanten Offset f/ 0 der Funktion werden Störeinflüsse im 
Objektbild berücksichtigt. Der Funktionsmittelpunkt und die Amplitude c sind für 
die weitere Betrachtung nicht relevant. 




(2.10) 
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Abb. 2.16. Linienbildfunktionerj (Schlitzbreite l px bzw. 100/un) in drei verschiedenes 
Allbildungstiefen t ({15; 30: 45} mm) in der Modellstudie (a) und der Detektormessung (b). 
An jedes Linienprofil (graue Balken) wird eine Gauss funkt ion (durchgezogene Linie) angepasst, 
die eine Funktionsbreite a und, im Falle der Messdaten, einen Offset i/o liefert. Uber da.s 
Bestiniintheitsinaß R 1 lä.sst sieh beurteilen, wie gut die Funktion durch die Sehätzfunktion 
erklärt wird. Die räumliche Auflösung in einem Objektbild steht im Zusammenhang mit der 
Breite des Linienbildprofils. 

Abb. 2.16 zeigt, beispielhaft Linienbildprofile in verschiedenen Abbildungstiefen. Die 
Linienbildprofile der Modellstudie sind komplett rauschfrei. Die Störsignale in den 
Messdaten hingegen sind gerade in kleinen Abbildungstiefen nicht unerheblich, sodass 
sich die Messwerte dort nur noch schlecht durch die Gaussfunktion erklären lassen. Zum 
einen liegt das an Reflexionen in der Septenmaske und zum anderen wurde für dieses 
Beispiel der breitere 100-Jtni-Schlitz verwendet. Mit steigender Abbildungstiefe verlieren 
sich die Störeinflüsse immer mehr und die Gaussapproximation verbessert sich, wie es 
auch das Bestimmt heitsmaß zeigt. 



Jedes Linienbildprofil in einer Messreihe zur tiefenabhängigen Auflösung wird letztlich 
auf die geschätzte Linienbildbredte o reduziert. Das hat zum einen den Effekt, dass 
die Linienbildprofile geglättet werden, andererseits vermeidet man, dass der konstante 
Offset nach clor Fouriertransfönnation zu einem hohen Signal bei u," = ü führt und 
die Normierung der MTF erschwert. Für jede gaussförinige LSF y(x) f repräsentiert 
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durch die Funktionsbreite er, ist die Fouriertransformiert e >'(>') ebenfalls wieder eine 
Gaussfunktion. Ks gilt der Zusaniinenhang 



mit u,' = 27TM.. Wird |y(u;)/F(0)| als normierte MTF betrachtet, so ist sie mit der 
LSF reziprok über er verbunden, d. h. verbreitert sich die LSF, so verschmälert, sich 
die MTF im Frequenzraum und die räumliche Auflösung nimmt ab. Es sei angemerkt, 
dass aufgrund der Glättung des Linienbildprofils, der Rauschunterdrückung und der 
Vernachlässigung weitere Einflüsse wie z. B. die endlich kleine Schlitzbreite der Lini- 
enquelle, diese Methode nicht zu einer exakten MTF des optischen Detektors führt. 
Vielmehr soll hier der Zusammenhang zwischen der Breite der Linienbildfunktion und 
der räumlichen Auflösung hergestellt werden. Die tiefenabhängige Veränderung der 
MTF, die sich aus der messbaren Verbreiterung der Linienbildfunktion ableiten lässt, 
ermöglicht eine recht genaue Abschätzung des tiefenabhängigen Auflösungsvermögens 
für das Abbildungsmodell. 

Neben der Bestimmung der räumlichen Auflösung, hat sich die Profilbreite der LSF auch 
als effiziente Messgröße für die Tiefenkalibrierung des Detektors erwiesen. Während 
für die Autiösungsmessung der Abstand der Linienquelle und die Abbildungstiefe auf- 
einander abgestimmt sind, wird für die Tiefenkalibrierung lediglich die Abbildungstiefe 
variiert, nicht jedoch der Quellenabstand. Als Folge der Defokussierung vi rgrößort sich 
die geschätzte Profilbreite in den Objektbildern. Ist die Linienquelle in einem Objektbild 
fokussiert. dann nimmt die geschätzte Profilbreite einen minimalen Wert an. Insgesamt 
zeigt diese Funktion einen parabolischen Verlauf, ihnen Minimum immer dann mit der 
tatsächlichen Abbildungstiefo der Linienquelle zusammenfällt, wenn die Vorwärts- und 
Rückwärtsabbildung mit derselben Fokuslänge entstanden sind (Abb. 3.14 im Ergeb- 
nisteil). Mit der Tiefenkalibrierung lässt sich also testen, ob für die R ückwärt sabbildung 
auch die tatsächliche Fokuslänge verwendet wird, mit der das MLA die Abbildung zuvor 
erzeugt hat. Aufgrund der nicht unerheblichen Ausdehnung des MLAs wird für die Rück- 
wärtsabbildung zu diesem Zweck eine effektive Fokuslänge (Absehn. 2.1.3) verwendet, 
die sich mit dieser Methode überprüfen lässt. Die Ergebnisse zur Auflösungsmessung 
und Tiefenkalibrierung mit einer Linienquelle finden sich in Abschn. 3.1.3. 

2.3.2 Signal-Rausch- Verhältnis 

Der optische Detektor sollte eine möglichst hohe Empfindlichkeit bei gleichzeitig gutem 
Signal-Rausch- Verhältnis (SNR, engl, signal-to-noise ratio) aufweisen. Die Empfind- 




(2.17) 
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lichkeit ist das grundlegende Vermögen, ein Quellsignal überhaupt zu detektieren und 
das Signal- Rausch- Verhältnis allgemein ein Qualitätsmaß für die delektierte Signalin- 
tensität relativ zum Rauschen. Ausschlaggebend für die Signahncssung sind im Falle 
des optischen Detektors die technischen Spezifikationen des CMOS-Sensors, die in 
Tab. A.2 im Anhang zu finden sind. Während der Belichtungszeit fallen Photonen 
auf die einzelnen Photodioden des Sensors, die dort mit einer bestimmten Quanten- 
ausbeute erst in Ladungsträger und anschließend durch die Ausleseelektronik in digitale 
Grau werte umgewandelt werden. Durch jeden dieser Schritte wird dem Signal zusätzlich 
ein Rauschanteil (z.B. Dunkelrauschen. Leserauschen. Quantenrauschen) mit unter- 
schiedlichen Abhängigkeiten (z. B. Zeit, Temperatur und Photonenanzahl) hinzugefügt. 

Der Detektor wird für zwei Messvarianten mit unterschiedlichen Signalstärken eingesetzt. 
Für die Oberflächenextraktion (Abschn. 2.43) kann die Signalintensität der Leucht- 
folien an die Sensorempfindlichkeit angepasst werden. Kurze Belichtungszeiten und eine 
Signalaussteuerung auf 80 % des Sättigungswertes haben sich für einen guten Objekt- 
kontrast in den Projektionsbildern bewährt. Das SNR ist bei den höchsten Grauwerten 
allgemein am besten [43). Für die eigentliche Signalmessung (Abschn. 2.4.4) lassen 
sich die Signalstärken nicht, beeinflussen und sind in der Regel sehr schwach. Gerade 
hier müsste eigentlich umgekehrt die Empfindlichkeit und das SNR auf die vorhandene 
Signalintensität abgestimmt werden. Aufgrund fehlender Alternativen für die groß- 
flächige Lichtdetektion, die für den Detektor erforderlich ist, bietet die Sensorauswahl 
hier aber keinen großen Spielraum. Als Signalquellen dienen in dieser Arbeit kleine 
Lumineszenzquellen und das Emissionslicht eines angeregten Fluoreszenzfarbstoffs. In 
beiden Fällen ist das Signal relativ schwach. Das Sensorsignal erreicht in der Regel 
dabei nur 5% des Sättigungswertes. Speziell für diesen Fall soll der Unterschied im SNR 
zwischen einem Fokuspunktbild und dem Abbildungsmodell untersucht werden. 

Da ein Fokuspunktbild sich direkt aus extrahierten Pixeln des Sensordatensatzes zu- 
sammensetzt, ist in diesem Bild auch das SNR zu erwarten, das im Sensordatensatz 
selbst zu finden ist. Die Rückwärtsabbildung hingegen berechnet ein Voxel aus einer 
Vielzahl von Sensorpixeln. Hierfür fließt ein wesentlich größerer Inforinationsanteil aus 
den Sensordaten i in und es findet gleichzeitig eine Mittelwert bildung statt. Die Frage 
ist nun, ob sich dadurch auch tatsächlich das Signal-Rausch-Verhältnis gegenüber einem 
Fokuspunktbild verbessen) lässt. Als Vergleichsmaß der Messbilder wird unter Annahme 
eines dominanten Dunkelrauschens bei schwacher Beleuchtung das SNR in Anlehnung 
an [44| definiert als 
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SNR = H — — ^. (2.18) 

Dabei ist </,„ der maximale (Jrauwert der Sigualquelle im Bild, r/ u der Mittelwert und 
o n die Standardahwoichung des Bildrauschens. Beurteilt werden soll damit die resul- 
tierende Signalhöhe der schwach leuchtenden Quellen relativ zum BÜdrauschen nach der 
Bilderzeugung. Interessante Variationsgrößen sind hierfür vor allem die Abbildungstiefe 
und die Belichtungszeit. Die Messergebnisse sind in Absch. 3.1.4 zusammengefasst. 

Das signalinvariante Dunkelrauschen des Sensors kann bereits in der Ausleseelektronik 
durch einen einstellbaren Schwellwertfilter abgeschnitten werden. Ein hoher Schwellwert, 
setzt aber gleichzeitig die Empfindlichkeit herab. Ein Vorteil ergibt sich hieraus vor allem 
dann, wenn die Quelleniutensität steuerbar und auf den Sensor abgestimmt werden 
kann. Dies ist. z.B. der Fall, wenn der Sensor in der digitalen Radiographie oder wie in 
dieser Arbeit für die Oberflächenextrakt iou eingesetzt wird. Ein zu niedriger Schwellwert 
andererseits, führt zu einem erhöhten Rauschinittelwert g a , der dann wiederum den 
Aussteuerbereich einschränkt. Für die Signalmessungen ist der Schwellwert gemäß 
Herstellerangaben für jeden Detektor so eingestellt worden, dass in jedem Dunkelbild 
(Belichtungszeit 1 s) ein Bildrauschen mit Mittelwert // = 20 und Standardabweichung 
o = 2 enthalten ist. 4 In der Regel wurde dann nach einer Signalmessung der Rausch- 
mittelwert, von den Sensordaten abgezogen. 



2.4 Experimenteller Scaimeraufbau mit mehreren Detektoren 

In diesem Abschnitt soll nun der voll rotierbarer Seunneraufbau mit mehreren Detek- 
tormodulen vorgestellt werden. Durch die kompakten Abmessungen der Detektoren 
und dem unmittelbar vor dem MLA beginnenden Objektraum können die Detektoren 
bis auf wenige Millimeter an das Untersuchungsobjekt herangebracht werden, sodass 
sich ein sehr platzeffizienter Scanneraufbau für die Kleintierbildgebung realisieren lässt. 
Das Ziel des Experimentalaufbaus ist es, das Zusammenwirken mehrerer Detektoren zu 
koordinieren, dabei das Volumen einer Labonnaus möglichst vollständig zu erfassen und 
schließlich die eingesetzten Methoden für die OberHächensignalmessung zu optimieren. 

Nachfolgend wird in Abschn. 2.4.1 der experimentelle Messaufbau sowie im Detail die 
Detektoranordnung auf der Rotationsplattform. die integrierte Objektbeleuchtung und 
die Scannersteuerung beschrieben. Die Koordinierung der Detektormodule erfordert 



Die Angaben beziehen sich auf den 12-Bit- Wertebereich des Sensor, der somit eiueu Umfang 
von 0 bis 41195 hat. 
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«'ine exakte Scannerkalibrierung, die in Abschn. 2.4.2 demonstriert wird. Schließlich 
benötigt eine fokussierte Rückwärtsabbildung die Abbildungstiefen im Objektraum. 
In Abschn. 2.1.3 wird ein Verfahren präsentiert, mit dein die Objektoberfläche über 
mehrere Winkelprojektionen sehr viel genauer bestimmt werden kann, als es aus einem 
einzelnen Sensordatensatz möglich ist. Abschließend wird in Abschn. 2.4.4 das komplette 
Messprotokoll für die Oberfl&cbensignalmessung vorgestellt. 



2.4.1 Scanneraufbau 




Steuersoftware 
(ImageJ-Plugin ) 



Energieversorgung 
für Leuchtfolie 




Abb. 2.17. Experimenteller Scanneraufbau im Überblick. Die vier Detektoren bilden das 
Kernstück und sind auf einer Rotationsplatt form ringförmig mit einem Durchmesser von etwa 
Sein befestigt. Energieversorgung und Steuerelektronik befinden sich außerhalb des Scanners. 
Der gesamte Messablauf wird über eine hierfür entwickelte Software (realisiert als Image J- 
Plugin) gesteuert. In dieser Abbildung sind lediglich zwei Sensoren an die Ausleseelektronik 
angeschlossen. 



Der Scanneraufbau ist in Abb. 2. 17 dargestellt und zeigt die Rotationsplattform 8 mit den 
vier montierten Detektormodulen sowie die externe Steuer- und Ausleseelektronik. Die 
Detektoren lassen sich fast vollständig um einen Objekthalter (Objektbett) drehen. Die 
genaue Detektoranordnung wird in Unterabschn. 2.4.1.1 behandelt. Für die Oberflächen- 
extraktion ist eine homogene Objektbeleuchtung notwendig, die in den Scanneraufbau 
integriert wurde und in Unter abschn. 2.4.1.2 erläutert wird. Die Scannersteuerung 



Der mechanische Aufbau der Rotatiunsplatiforiu ist in der Diplomarbeit von K. 
Schneider |4Fj] entstanden. 
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(Unterabschn. 2.4.1.3) umfasst. das .synchrone Auslesen aller vier Detektoren sowie die 
darauf abgestimmte Ansteuerung des Schrittmotors. Die Ablaufprotokolle sind in einem 
IniageJ-Plugin implementiert, mit denen die automatische Datenaufhahme z.B. für die 
Oberflächenextrakt ion oder Signalmessung durchgeführt werden kann. 

Für die Auswertung der Sensordaten aus einer kompletten Scannerrotation fallen große 
Mengen an Daten an; so ist bei 153 Winkelstellungen für einen einzelnen Detektorda- 
tensatz mit 306 MB Datenvolumen zu rechnen. Um kürzere Rechenzeiten zu erzielen 
wird für ein Objektbild nur noch jedes achte Voxel in x- bzw. y-Richtung berechnet, 
sodass sich eine Bildgröße von 64 x 128 px ergibt. Der Voxelabstand vergrößert sich 
dann auf 0,384 nun und die Grenzfrequenz in einem Objektbild verringert sich ent- 
sprechend auf f g = 1,3 'p/mm. Dieser Schritt lässt sich durch die Auflösungsmessungen 
aus Abb. 3.13(b) rechtfertigen, da schon ab Abbildungstiefen von t > 30 mm kaum 
noch höhere Raumfrequenzen aufgrund der MLA-Geometrie dargestellt werden können. 
Für die Oberflächenextraktion wird später eine konstante Abbildungstiefe t = 38 mm 
verwendet. Die Oberflächensignalmessung ließe sich durchaus mit einem leicht kleineren 
Voxelabstand (z.B. 0.192mm) durchführen, da die Abbildungstiefen hier im Bereich 
von 1T> mm - 30 mm liegen und eine etwas bessere räumliche Auflösung zulassen. Für die 
abschließende Beurteilung des Scanneraufbaus ist. diese Option jedoch nicht ausgenutzt 
worden. Es bleibt zu bemerken, dass auch ein Objektbild mit dem o.g. Voxelabstand 
nach wie vor ein besseres räumliches Auflösungsvermögen und SNR als ein vergleichbares 
Fokuspunktbild zeigt, da über die Rückwärtsabbildung ein Großteil der Sensorpixel 
für die Bilderzeugung verwendet wird. Für die Objektbildberechnung wird lediglich 
die einfache Rückwärtsabbildung verwendet, da durch die iterative Anwendung des 
Abbildungsmodells weder für die Oberflächenextraktion noch für die eigentliche Signal- 
messung ein Vorteil erwartet wird. 6 

Der Scannerobjektraum, um den alle vier Detektoren positioniert sind, wird nun analog 
zum Objekt räum des Abbildungsmodells durch ein diskretes Volumen repräsentiert 
und nachfolgend schichtweise dargestellt. Im Zusammenhang mit dem Scanneraufbau 
wird der Seannerobjektraum im Weiteren auch einfach als Ohjt ktmum bezeichnet. Die 
Voxelintensitäten berechnen sich je nach Kontext aus einem oder mehreren Detektor- 
datensätzen. Ein Detektordatensatz bestellt dabei aus allen Sensordatensätzen eines 



Als Vorgriff auf den Ergebinsteil (Abschn. 3.1.2.2) sei hier angemerkt, dass die iterative 
Anwendung des Abbildungsmodclls lediglich in kleinen Abbildungstiefen Vorteile bringt, 
die den zusätzlichen Rechenaufwand rechtfertigen würde. Für die Abbildungstiefen, die im 
Scanneraufbau zu erwarten sind, ist dies nicht mehr der Fall. 



Copyrighted material 



42 



2. Material und Methoden 



Detektors, die während einer Scaimerrotation in ver.schiedenen Winkelstelinngen aufge- 
nommen wurden. Für die Bezeichnung der Schichten werden in späteren Darstellungen 
folgende Konventionen verwendet: 

• Transversalebene (TE)/z-Schicht: x-y-Ebene. adressiert über z-Position in [px] 

• Sagittalobene (SE)/x-Schicht: y-z- Ebene, adressiert über x-Position in [px] 

• Koronarebene (KE)/y-Schicht: x-z-Ebene, adressiert über y-Position in [px] 

Weiterhin werden im Scanneraufbau für jeden Detektor die Objektbilder im Koordi- 
natensystem des Abbildungsmodells berechnet, das dann mit dem Detektor rotiert. 
Aus diesem Grund werden im weiteren Textvorlauf sowohl die winkelabhängigen Sen- 
sordaten als auch die direkt daraus berechneten Objektbilder als Detektorprojektionen 
bezeichnet. In einer Detektorprojektion ändert sich dann auch die Bezeichnung der y- zur 
z-Achse und der x- zur s-Achse gegenüber dem Objektbild, um für die winkelabhängige 
Rückprojektion konform mit dem Koordinatensystem des Scannerobjektraums zu sein. 
Die Intensitäten des Scannerobjektraumes werden für jeden rotierenden Detektor im 
jeweiligen Detektorktwrdmatenaysteme.n (DKS 1 bis 4) separat berechnet. Uber die 
Scannerkalibrierung lassen sich alle DKSe in das Hauptkoordinatensystem (HKS) über- 
führen, das durch die OberHächenberechnung und Detektor 0 vorgegeben wird. Die 
Konvention und Orientierung der Detektoren und des Scannerobjektraumes sind im 
Anhang in Abb. A.3 und Abb. A.4 graphisch skizziert. 

2.4.1.1 Detektoranordnung 

Die Rahmenkonstruktion der Detektormodule soll einen möglichst robusten Zusam- 
menbau der justierten Komponenten gegenüber Erschütterungen gewährleisten, wie sie 
bei dei S> anneiT >tation zu ovw irten sind. 1 ür eine raumsparende D< tekt Bauordnung 
auf der Rotationsplattform bietet sich ein regelmäßiges Polygon an. dessen Kantenlänge 
durch die Rahmenbreite bestimmt ist und dessen umschlossenes Volumen ausreichend 
Platz für eine Labormaus bereithält. Eine Anordnung der vier Detektoren zu einem 
Hexagon, skizziert in Abb. 2.18(a), ist hier ein guter Kompromiss. Die offene Seite 

iles Detektorrings eignet sich zudem für eine einfa. ln> H - ■ t . • 1 1 1 ; i "■ ■ i : 1 1 ^ ab/, ibi'.l l.-n 

Objektes und für die spätere Anbringung eines Anregungslasers. Mit einer zentrierten 
Anordnung der Detektoren ist es möglich, innerhalb einer kompletten Umdrehung 
redundante Projektionen mit verschiedenen Detektoren aufzunehmen. Alternativ kann 
bei dieser Anordnung schon über eine Teilrotation eine komplette Winkelabdeckung des 
Objektes erreicht werden. Das Gesichtsfeld 7 (FOV, engl, field-of-view) der zentrierten 

7 Das FOV entspricht hier dein Bereich, der aus allen möglichen Projcktionswinkelu durch 
alle Detektoren gemeinsam abgedeckt wird. 
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Del ektoranordnung en( spricht im Querschnil I einem Kreis, der als Durchmesser die 
Breite der aktiven Detektorfläche hat. Nachteilig erweist sich diese Detektoranordnimg 
bei der Oberflächenberechnuug für Objekte, die diesen Durchmesser überschreiten 
(Abschn. 2.4.3). Es sei an dieser Stelle erwähnt, dass sich diese Anordnung durchaus zur 
Signalmessung eignet, falls die Objektoberfläche a-priori bekannt ist. da für die Bildbe- 
rechnung mit dem Abbildungsmodell keine komplette Objektabdeckung erforderlich ist. 

Eine FOV- Vergrößerung lässt sich erreichen, indem zwei Detektoren aus der Hexagonal- 
struktur so herausgedreht werden, dass sich ihre individuellen FOVs nur noch teilweise 
überschneiden. Diese Anordnung ist in Abb. 2.18(1)) skizziert. Beide Projektionen 
können mittels der Methode in Absehn. 2.4.3.1 zu einer einzigen zusammengesetzt 
werden, mit der dann auch eine Labonnaus für die Oberflächenberechnuug abgedeckt 
wird. Diese spezielle Detektoranordnung ermöglicht somit eine Oberflächenextraktion 
für Labormäuse, die alleine mit dem Scanner durchgeführt werden kann. Dadurch kann 
der Scanner autark für die Oberflächensignalmessungen eingesetzt werden. Nachfolgend 
werden die Detektoren 2 und 3 als Obcrflächendctektorcn und die Detektoren 1 und 4 
als Sifinnldctektotvn bezeichnet. 

Die Sensoren sind über eine Kabelhaltorung auf der Plattform zugfrei an die Ausle- 
se •lektruuik angeschlossen. Für diesen Fxperimentalaufbau ist noch akzeptabel, dass 
die vier Sensorkabel bei der Platt formdrehung teilweise manuell nachgeführt werden 
müssen, jedoch erfordern zukünftige Anwendungen eine verbesserte Kabelkonstruktion. 
Angesteuert wird die Ausleseelektronik über einen digitalen Framegrabbor (PXDlüüü, 
CyberOptics Semiconductor. Boaverton, OR. USA). 

2.4.1.2 Objektbeleuchtung 

Die Objektbeleuchtung ist ein wichtiger Bestandteil für die Oberflächenextraktion mit 
dem Ziel einer möglichst homogenen Ausleuchtung der Objektsilhouette. Als Licht- 
quelle werden Elektrohmiineszenzfolien (Conrad Electronic SE. Hirschau. Deutschland) 
eingesetzt, die auch schon zuvor als Leuchtmittel für die Test- und Schlitzschablone 
zur Auflösungsmessung Verwendung fanden. Eine Leuchtfolie ist eine gleichmäßig ab- 
strahlende Kaltlichtquelle, in der eine Schicht z. B. aus Phosphorverbindungen über ein 
wechselndes, elektrisches Feld zum Leuchten angeregt wird (Elektrolumineszenz) . Die 
Folie kann in beliebiger Form zugeschnitten werden, ohne ihre Funktion zu verlieren. 
Dadurch lässt sich die Objektbeleuchtung direkt in den Scanneraufbau integrieren, 
indem die Leuchtfolie passend auf die Detektorrahnien zugeschnitten und um die aktive 
Detektorfläche angebracht wird. Damit wird eine diffuse Beleuchtung erreicht, ohne 
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(a) (b) 



Abb. 2.18. Detektoranordnung auf der Rotationsplattforni: (a) Zentriert«' Detektoranordnung 
mit einem FOV (0 - 2.5cm). das durch die Breite der aktiven Detektorfläche vorgegeben ist. 
(b) Durch die gedrehten Detektoren 2 und 3 lässt sich das FOV vergrößern (0 > 4 cm). um für 
die Oberflächenextraktion eine Labormaus abzudecken. Die zusammengesetzten Projektionen, 
die das vergrößerte FOV erzeugen, werden in Absehn. 2.4.3.1 einem virtuellen Detektor I) 
zugeordnet. Nach wie vor können alle Detektoren für die Signalniessung eingesetzt werden. 

dass die Obcrllächendetcktoren (Det. 2 und 3) geblendet sind. Die Signaldetektoren 
(Det. 1 und 1) hingegen kommen bei der Oberflächenextraktion nicht zum Einsatz. 
Durch die gegenüber angeordneten Leuchtfolien können hier keine vernünftigen Auf- 
nahmen entstehen. Die Leuchtfolien sind in Abb. 2.1?)(a) im ausgeschalteten Zustand 
auf den Detektoren zu erkennen und in (b) im Einsatz für die Oberflächenextraktion 
einer Nacktmans gezeigt. Durch die seitliche Holouchtung ist sichergestellt, dass die 
Ob jektkonturen in einer Projektion deutlich sichtbar sind. Die Leuchtfolien werden 
parallel über eine Weehselfrequenz (GOÜ Hz) außerhalb der Plattform betrieben, zu dem 
die dünnen Anschlusskabel entlang eines Sensorkabels geführt werden. 

2.4.1.3 Scannersteuerung 

Die Rotationsplattform wird durch einen rückseitig befestigten Schrittmotor bewegt, 
der sich über eine programmierbare Schrittmotorkarte (Elektroniklabor DKFZ. Hei- 
delberg. Deutschland) ansteuern lässt. Die Motordrehimg wird im Verhältnis 1:4 auf die 
Plattform übertragen, sodass eine komplette 3G0 "-Drehung 1600 Schritten entspricht. 
Als inkrementelle Winkelrotation hat sich eine Drehung um 10 Schritte (= 2.25°) als 
hinreichend genau bezüglich der mechanischen Toleranz und als ausreichend für die Ob- 
jektrauinabtastung zur Oberflächenrekonstruktion erwiesen. An der Rotationsplattform 
angebracht ist ein Motorschalter, durch den eine feste Anfangs- und Endposition vor- 
gegeben ist . Zwischen beiden Positionen bestellt konstruktionsbedingt eine Lücke von 
etwa 10°. Die deswegen fehlenden sieben Projektionen lassen sich für die Oborfiachenex- 
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(a) (b) 

Abb. 2.19. (a) Anordnung der vier Detektoren auf der Rotat ioiisplattfonn. Für die Ober- 
Hächenextraktion werden ausschließlich die Detektoren 2 und 3 (Oberflächendetektoren) 
verwendet, die nicht auf das Drehzentrum zentriert sind. Nach dem Zusammensetzen beider 
Projektionen lässt sich dadurch das FOV auf Mausgröße anpassen. Zur anschließenden Signal- 
niessung können wieder alle Detektoren eingesetzt werden. (1>) Auf jedem Detektorralnnen sind 
zugeschnittene elektrische I.enchtfolieii befestigt, die das Objekt (hier eine Nacktiiiaus) nahezu 
homogen und diffus beleuchten. In diesem Zustand werden nur die Oberflächendetektoren 
betrieben. 

traktion jedoch einfach interpolieren (s. Abschn. 2.4.3.1). Für die Oberflächenextraktion 
werden somit gewöhnlich 153 Winkelpositionen angefahren, die entsprechende Anzahl 
an Detektorprojektionen aufgenommen und automatisch gespeichert. Die Winkelab- 
weichung nach dein Anfahren aller Positionen beträgt ±1 Schritt. Die Aufnahmezeit 
der Sensoren kann entsprechend der Vorgaben gewählt werden und richtet sich nach 
der Anwendung. Typisch erweise dauert das Erzeugen eines Detektordatensatzes für die 
Oberflächenextraktion etwa 2min. Plattformrotation und Datenaufuahme werden dabei 
automatisch über Ablaufprotokolle koordiniert, die synchron den Schrittmotor und die 
Sensorelektronik steuern. 

2.4.2 Scannerkalibrierung 

Der Ilauptaspekt der Scannerkalibrierung liegt auf der genauen Positionsbestimmung der 
einzelnen Detektoren relativ zur Rotationsachse des Scanners. Erst mit den bekannten 
Lageparametern lassen sich die einzelnen DKSe ineinander überführen. Zusätzlich 
werden auch Konstruktionstoleranzen ausgeglichen, die durch den St-anuerzusnmmenbau 
entstanden sind. Mit einem dünnen Stäbchen werden in jedem DKS die Lage der Rota- 
tionsachse und die Winkelposition des Detektors auf der Plattform bestimmt. Dafür 
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wird das Stäbchen im Scannerobjekt räum fixiert, die Leuchtfolie des zu kalibrierendem 
Detektors eingeschaltet, und ein kompletter Datensatz (Kalibrierungsdatensatz) über 
eine volle Umdrehung aufgenommen. Dieser Vorgang wird für jeden Detektor einzeln 
wiederholt. 

Für die Scannerkalibrierung (wie auch später für die Oberflächenberechnung) wird 
dann mit dem Abbildungsmodell aus jeder Winkelstellung je ein Objektbild mit einer 
festen Abbildungstiefe / 38mm berechnet. Der Objektfokus in dieser Projektion 
ist somit nicht auf die Oberfläche sondern auf die Kontur des beleuchteten Objektes 
gerichtet, die sich dadurch gut abgrenzbar gegenüber dem Hintergrund darstellt. Diese 
Detektorprojektion wird dann näherungsweise als senkrechte Winkelprojektion des 
Objektes, d. h. ohne perspokt ivische Verzerrung, auf die 1 Mekt .rlläche bei rächtet. 

Ein Kalibrierungsdatensatz besteht aus 153 Projektionen, in denen dann die Lage 
des Stäbchens sichtbar ist. In Abb. 2.21) sind beispielhaft vier Detektorprojektionen 
dargestellt. Um die genaue Position der Rotationsachse daraus zu bestimmen, ist es 
vorteilhaft, den Datensatz in Sinogramme umzusortieren. Dafür wird jede Zeile einer 
Detektorprojektion in Abhängigkeit der Winkelst eilung 0 in der entsprechenden z- 
Schicht neu angeordnet. Als Resultat beschreibt das Stäbchen dann in jeder z-Schicht 
einen Sinusverlauf, der ebenfalls beispielhaft in Abb. 2.20 dargestellt ist. Aus diesen 
Sin. grammeu wird jetzt per geeigneter Schätzfunktion sowohl die Rot nl ionsachse als 
auch die Winkelposition des Detektors bestimmt. Hierfür wird jedes Sinogramm in ein 
Binärbild umgewandelt und mittels Regressionsrechnung an die Funktion 



angepasst. Für jede z-Schicht lässt sich so die Position der Rotationsachse s g in den Pro- 
ji kl i< "h u m zen in' 1 1 i 'in ■ ,i I i-v vii' ■ ' W in 'Ij;i <h1 i> n . ;v m'.i 1 < 'Ii; . Der Radialal ist und < 
des Stäbchens zum Rotationszcntruin ist hier nebensächlich. Wiederum beispielhaft 
hierfür sind in Abb. 2.20 für zwei z-Schichten sowohl die Messwerte (graue Punkte) 
als auch die dazugehörige Schätzfunktion (durchgezogene Linie) im selben Diagramm 
dargestellt. In der Schätzung wird eine volle Umdrehung mit Hit) Projektionen berück- 
sichtigt, obwohl im Datensatz selbst nur 153 Projektionen vorhanden sind. 

Da der Verlauf der Rotationsachse entlang der z-Aehse durchaus variieren kann, wird 
die Position der Rotationsachse $- in zwei Sinograminen (10px; 100 px) bestimmt und 
linear interpoliert mit 




(2.10) 
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Abb. 2.20. Projektionen eines dünnen Stäbchens aus unterschiedlichen Winkelstellungen von 
Detektor 4. Das Streulicht der Objektbelcuchtung wird über einen Schwellwert abgeschnitten. 
Die Lage der Rotationsachse lässt sich schichtweise aus den Sinogranmieii durch die Anpassung 
an eine Sinusfunktion und den n ÜfFsrt schätzen. Dn konstruktionsbedingt nur 153 von 160 
Projektionen aufgenommen werden, fehlt ein kleiner Teil der Messdaten. Die Winkelposition 
des Detektors berechnet sich aus dem Mittelwert beider Schätzungen zu ^4 = 2,0 rad. Ist das 
Rutationszentruni korrekt bestimmt, so ist die Lage des Stäbchens im Sca nnerub jok t ra 1 1 in 
nach der Rückprojektion sichtbar. 

•*(«) = * 10 ° ~ * W (z - 10) + *io. (2-20) 

Mit der nun bekannten Lage der Rotationsachse für einen Detektordatensatz lässt 
sich der Scannerobjektrauin über eine gefilterte Rückprojektion berechnen, die in 
Abschn. 2.4.3.2 noch eingehend beschrieben wird. Die Kalibrierung war erfolgreich, 
wenn das Stäbchen im Objektraum dann schichtweise wieder als Punkt, erscheint, wie 
in Abb. 2.20 auf der rechten Seite zu sehen ist . 

Die W inkelposition <p eines Detektors ist mit dieser Schätzung nur relativ grob zu 
bestimmen. Sie w r ird als Mittelwert aus den beiden Sinogrammschätzungen berechnet 
und als invariante Größe für den jeweiligen Detektor festgesetzt. Zusammen mit der 
Lage der Rotationsachse und einem z-Achsen-Offset 8 können dann die verschiedenen 
DKSe mit den bekannten Winkelpositionen ineinander überführt werden. Die Übersicht 
der ermittelten Werte und deren Toleranzen für alle vier Detektoren sind in Tab. A. 1 
zusammengefaßt. Auf die Besonderheil bei der Kalibrierung der Oberflächendetektoren 
wird in Abs« -lt. 2. 1.3. 1 oiiigrg.iugi'n. 



8 Der Offsct Az für die Z- Achse wird über die gleichzeitige Abbildung einer Kante in allen 
vier Detektoren ermittelt. 
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2.4.3 Oberflächenberechnung 

Obwohl die Tiefen in form ation schon aus einer einzelnen Sensoraufnahme gewonnen 
wen Ion kann, v\ \< in U>schn. 2.2.G bes. Iniobeit, so sind die daraus abgeleiteten Uj 
bildungstiefen für die Signalrückprojoktion noch zu grob. Zudem befindet sich das 
Objekt im Scanner in einer Entfernung zu den Detektoren, in der die Anwendung 
dieser Methoden zusätzlich ersehwert wird, sodass deren Einsatz hier nicht sinnvoll 
ist (s. Ergebniskapitel Abschn. 3.1.1.2). Dagegen besteht nun mit dem Scanneraufbau 
die Möglichkeit mehrere Winkelprojektionen eines Objektes für die Bestimmung der 
Abbildungstiefe zu verwenden. In diesem Abschnitt wird gezeigt, wie und unter welchen 
Voraussetzungen die beleuchtete Objektoberfläche berechnet werden kann. 

Eine Bedingung für die Oberflächenberechnung ist, dass jode Detektorprojektion die 
gesamte Objektbreite abdeckt. Da ein einzelner Detektor diese Anforderung für eine 
Labonnaus nicht erfüllt, wird in Abschn. 2.4.3.1 gezeigt, wie sich ein ausreichend großes 
FOV aus zwei Projektionen zusammensetzen lässt. Für die Oberflächenberechnung wird 
dann eine gefilterte Rückprojektion analog zur Rekonstruktion in der Computertoino- 
graphie (CT) eingesetzt, die in Absch. 2.4.3.2 beschrieben wird. Der Haupt unterschied 
dabei ist. dass den Signalintensitäten in den Detektorprojektionen die Lichtreflexion an 
der Objektoberfläche zugrunde liegt und nicht die Schwächung eines Transmissionss- 
trahls wie in der CT. Die Segmentierung erfolgt im Objektraum nach der Rückprojektion 
mittels eines Schwell wertes, der eine Isooberfläche (engl, iso surface) des Objektes liefert. 
Mit der nun bekannten Objektoberfläche lässt sich dann für jede beliebige Projektion 
mit einem Referenzabstand tnz die Abbildungstiefe zum Sensor berechnen und für die 
anschließende Signalmessung verwenden. 

2.4.3.1 Vergrößerung des Gesichtsfeldes 

Eine wichtige Voraussetzung für die erfolgreiche Oberflächenberechnung ist. dass mit 
jeder Detektorprojektion die gesamte Objektbreite abgedeckt wird. Mit der zentrierten 
Detektoranor In in° isl ler Durchmesse] des I OV in eine] transversalen Ob jekt mimi- 
sch icht auf die Breite der akt iven Detektorfläche von 2,5 cm begrenzt (s. Abb. 2.18(a)). 
Für entsprechend kleine Objekte isl das ausreichend, jedoch nicht für eine Labormaus. 
Um das FOV an Mausgröße anzupassen sind im Scanneraufbau die Detektoren 2 und 
3 etwas aus der zentrierten Anordnung herausgedreht (s. Abb. 2.18(b)), sodass sich 
aus diesen beiden Oberflächendetektoren eine vergrößerte Detektorprojektion zusam- 
mensetzen lässt. Diese zusammengesetzte Detektorprojektion wird einem virtuellen 
Detektor 0 zugeordnet, der ebenfalls über eigene Lageparameter verfügt. 
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Abb. 2.21. Lageparameter der übeifiächcudetektoren: Die Detektoren 2 und 3 decken durch 
ihre Anordnung (s. Abb. 2.18(1:»)) jeweils einen separaten T il des Objektraiuncs al<. dei sich 
nur wenig Überschneidet. Das Kalibrierstabchcn ist somit nur teilweis*" in den Projektionen 
sichtbar. Der dadurch ebenfalls in < Um i Siiiograiuinen lückenhafte Siuusverlauf kann dennoch 
zur Positionsschätzung der Rotationsachse und der Winkelposition der Detektoren verwendet 
werden. 

Wie für die übrigen Detektoren im Seanneraufbau, so lassen sieh auch für die Ober- 
flächen detektoren die Lageparaineter über ein dünnes Stäbchen linden. Beispielhafte 
Detektorprojektionen und Sinogramme sowie der Verlauf der Mess- und Schätzdaten 
sind in Abb. 2.21 dargestellt. Da beide Detektoren jeweils getrennte Bereiche des 
Objektraumes abbilden und nicht mehr auf die Rotationsachse zentriert sind, ist das 
Kalibrierstäbchen nur in einem Teil der Detektorprojektionen zu sehen. Folglich ist 
dadurch dann der Siuusverlauf in den Sinogrammen lückenhaft. Dennoch können auch 
hier über Gl. 2.19 die Lage der Rotationsachse und die Winkelpositionen der Detektoren 
geschätzt werden. 

Mit den bekannten Lngeparametorn werden nun die beiden Detektordatensätze als Sino- 
gramme zusammengefügt. Das Vorgehen ist in Abb. 2.22 in einer z-Schicht demonstriert. 
Die beiden einzelnen Sinogramme lassen sich dabei durch eine einfache Translationsver- 
schiebung zusammensetzen, indem die Positionen der Rotationsachse auf der s-Achse 
zur Deckung gebracht und die Winkeldifferenz der Detektorpositionen y 3 - ip-} s» 28 0 
auf der 0-Achse berücksichtigt wird. Für jede /.-Schicht können derart die disjunkten 
Objektraumbereiche der beiden Oberflachendetektoren kombiniert werden. Durch den 
periodischen Verlauf lassen sich die Sinogrammdaten innerhalb einer z-Schicht. nach 
einer vollen Umdrehung umsortieren. Da zu einer vollen Umdrehung bei diesem Scanner 
konstruktionsbedingt sieben Projektionen fehlen (Abs< Im. 2.4.1.3), entstellt ein Lücke 
in jedem zusammengesetzten Smogramm, die mit dieser Methode sogar noch linear 
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interpoliert werden kann. 9 Die zusammengesetzten Sinogramme lassen sich anschließend 
wieder in Detektorprojektionen umsortieren. 



s[px] 




(a) (b) (c) 



Abb. 2.22. Zusammensetzen der Sinogramme: (a) Die disjiinkten Objektrauminforinationeii 
der Obcrflächcndetektoren werden über die Position der Rotationsachse und ihren Winkelunter- 
schied zusammengefügt, (b) Zusätzlich lassen sich Sinograininzeilen aufgrund der Periodizität 
umsortieren. Konstruktionsbedingt fehlen sieben Projektionen, deren Bereich sich hier mm in 
jeder z-Schicht interpoliert lässt. (c) So entstehen Sinogramme mit jeweils 128 X 1 53 px und 
einem dadurch vergrößerten Gesichtsfeld. Die Projektionen werden einem virtuellen Detektor 
0 zugeordnet. 

Der zusammengesetzte Datensatz von Detektor 0. dargestellt in Abb. 2.2M. besteht jetzt 
aus L53 Detektorpro jektionen mit 128 x 128 px. Hier ist nun in joder Detektorprojektion 
<las Kalibrierstäbchen zu sehen und die Sinogramme sind lückenlos. Die festgelegte Lage 
der Rotationsachse s Hot 63 ± 0,5 px und die Winkelposit ion tpa - 0,3 db 0.1 rad von 
Detektor 0 werden durch die Parameterschätzung bestätigt. Die Rückprojektion in den 
Objektraum zeigt dann auch das Stäbchen wieder punktförmig in jeder /-Schicht. Das 
Koordinatensystem dieses Objektraumes definiert zugleich das HKS. da aus diesem 
Datensatz später die Objektoberfläche segmentiert wird. Mit diesen zusammengesetzten 
Detektorprojektionen lassen sich nun nachfolgend transversale Objektraumschichten 
mit 128 x 128 px berechnen. Mit einem Kinzeldotektor waren zuvor lediglich 64 x 64 px 
möglich. 



Grundsätzlich ließe ein geschicktes Einfügen und Interpolieren der Sinogramme zu. dass 
ein komplettes Sinograniiu mit 128 x IGOpx entsteht. Nachfolgend werden jedoch nur 
Sinogranime mit 128 x 153 px erzeugt, die für die Oberflächeiibestiinmiuig ausreichend sind. 
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Abb. 2.23. Die zusammengefügten Sinogramme bzw. Detektorprojektionen (Detektor 0) 
zeigen nun vollständig das Kalibrierstäbchen. An den Uherlappbereichen kommt es zu Kanten, 
die für die Oberflächenextrakt ioii später aber nicht hinderlich sind. Die Paraineterschätzuiig 
bestätigt die zuvor festgelegten Lagepa rameter. über die Rückprojektion in den Objektraum 
können die Lageparameter überprüft werden. Eine transversale Schicht umfasst nun einen 
Bereich von 128 x 128 px. 

2.4.3.2 Gefilterte Rückprojektion und Abbildungsticfcn 

Die Oberflächenextraktion erfolgt unabhängig davon, ob es sich um eine zusammen- 
gesetzte Projektion oder einen Einzeldetektor handelt, immer auf dieselbe Art: Wie 
schon bei der vorangegangenen Scannerkalibriorung wird mit dem Abbildungsmodell 
eine Detektorprojekt ion mit fester Al>l>ildungstiefe t = 38 mm berechnet, um die Ob- 
jektkontur deutlich darzustellen. 10 Die aus unterschiedlichen Winkeln aufgenommenen 
Detektorprojekt Ionen werden anschließend orthogonal in den Objektraum Zurückpro- 
jiziert, nachdem durch eine Filterung der Konturenkontrast verstärkt wurde. Damit 
später aus dem Objektraum erfolgreich eine Isooberfläche segmentiert werden kann, 
müssen folgende Voraussetzungen erfüllt sein, die experimentell gefunden wurden: 

• Bekannte Projektionsgeometrie: Die Winkelstellung jeder Projektion und die Lage der 
Rotationsachse müssen genau bekannt sein. Die Geometrie kann Über die Kalibrierung 
in Abschn. 2.1.2 bestiminl werden. 

• Homogene und diffuse Ohjekth ■ leiichiung: Eine Oberfläche im Objektraum soll 
möglichst gleichmäßig beleuchtet werden, ohne dass die beteiligten Detektoren durch 
die Lichtquellen geblendet werden. Zur Beleuchtung werden elektrische Leuchtfolien 
verwendet, die an den Detektorrahmen befestigt sind (s. Abschn. 2. 1.1.2). 

• Homogene Lichtreflexion an der Objektoberfläche: Die Reflexionseigenselmftou an der 
Objektoberfläche sind die Grundlage für die Oberflächenextraktion. Eine gleichmäßig 

111 Die Objektbilder sind für diesen Fall somit immer auf die Objektsilhouette fokussiert. 
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reflektierende Oberfläche lässt sich nach der gefilterten Rückprojektion über einen 
Schwellwert segmentieren. Eine Nacktmaus weist diese Eigenschaft typischerweise 
auf. 

• Transversales FOV in Objektgröße: Ist ein Objekt breiter als der Sensor, so ist die 
Kontur nur teilweise in den Projektionen sichtbar. Dies führt zu Artefakten nach der 
Rückprojektion und verhindert eine korrekte Oberflächensegmentierung. Das trans- 
versale FOV für den Scanneraufbau lässt sich über die Methode in Abschn. 2.4.3.1 
an die Mausgröße anpassen. 

Der Objektraum wird entlang der z- Achse für jede transversale Schicht ft r (x,y) aus den 
entsprechenden Projektionszeilen p e über eine digitale gefiltert Rückprojektion (FBP, 
engl, filtered backprojeetion) berechnet. Die Methode ist analeg zur Bildrekonstruktion 
in der CT [37. 46] implementiert. Es sei darauf hingewiesen, dass die FDP in dieser 
Arbeit auf Reflexionsprojektionen anstatt auf Schwächung eines Transinissionstralds 
angewendet wird. Das heißt, das Verfahren entspricht nicht mehr einer Radontransfor- 
mation, wie sie für die CT gültig ist. Die Anwendung der FBP beschränkt sich hier 
nur auf die Oberflächenextraktion unter den o.g. Bedingungen. Zudem ist für die FBP 
mit den Reflexionsprojektionen, anders als bei einer Transmissionsmessung. eine volle 
Drehung um das Objekt notwendig. Die Rückprojektion erfolgt unter Annahme einer 
Parallelst .rahige* un« -t rie. 

Für die Filterung hat sich die Funktion von Shepp und Legan, dargestellt in Abb. 2.24, 
ids geeignet erwiesen Mit D>«tektorabstand As ist der digital* Faltungskern durch 



h SL (kAs) = -^-3 • T7J5 T (2-21) 



J 

n*(Asf 4fe2-l 

beschrieben. Die gefilterte Projektionszeile p© wird über die Faltungssumme 



K 

f) 0 {nAs) As ^ i> e (nAs — kAs)h SI (kAs) (2.22) 
k —K 

mit einer Filterlänge von K = 1 berechnet, die sich als adäquater Kompromiss zwischen 
Konturverstärkung und Artefakterzeugung im Objektraum erwiesen hat Anschließend 
wird jede gefilterte Projektionszeile orthogonal zur Detektorebene in den Objektraum 
unter dem entsprechenden Winkel O rückprojiziert: 

Uix, y) = ^ p B (x cos 0; + y sin G, ) . (2.23) 
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Abb. 2.24. Filterfunktion von Sliepp lind Logaa im Ortsraum (a). berechnet als analoges 
hsi(s) und digitales Filter hsL(kAa) mit Detektorabstand As = 1, und dargestellt im 
Frequenzraum (b). 



Die Pixelwerte von p& worden entlang einer Geraden s j-cosÖ + ysinÖ für ein festes 
(-) zu den Pixelwerten der Objektraumschicht f tr addiert und hierfür entsprechend 
bilinear interpoliert. Auf diese Weise wird mit allen M Projektionen verfahren. 



Entscheidend für die Oberflächenextraktion ist nun. dass sich in den rückprojizierten 
Objektraumschichten die Objektkontur über einen Schwellwert segmentieren lässt. Es ist 
hierfür hilfreich, wenn schichtweise die Randartefakte ent fernt und auf den verbleibenden 
Maximalwert normiert wird. Nachfolgend werden die Eigenschaften der vorgestellten 
Methode anhand verschiedener Studien erläutert. 



Das Gehäuse einer Meeresschnecke, dargestellt in Abb. 2.25, weist eine relativ komplexe 
Oberfläche auf und genügt weitestgehend der geforderten gleichmäßigen Lichtreflexion. 
Für die Objektgröße ist das Gesichtsfeld eines einzelnen Detektors ausreichend. Die Ober- 
flächenkontur ist in den normierten Objektraumschichten nach der FDP gut erkennbar. 
Uber einen Schwellwert lässt sich mittels Renderingsoftware |47J eine Isooberfläche dar- 
stellen, die das segmentierte Objekt repräsentiert. Die dunkleren Stellen des Gehäuses 
führen zu Eindellungen der Oberfläche. Durch die Normierung treten die Haltestäbe 
des Objektbettes an den Außenseiten stärker hervor. 

Das Mausphantom in Abb. 2.26 zeigt die Grenzen der optischen Oberflächenextraktion. 
Aufgrund der farblich unterschiedlichen Oberfläche lässt sich diese nicht hinreichend 
bestimmen. Wiederum ist ein einzelner Detektor ausreichend für die Objektgröße, doch 
schon in den Detektorprojekt ionen ist die Objektsilhouette nur ungleichmäßig abge- 
bildet. Nach der FDP dominiert die besonders hell reflektierende Unterseite, sodass 
sich nicht einmal innerhalb einer Transversalschicht ein einheitlicher Schwellwert linden 
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Abb. 2.25. Ohcrfl/i cliencxt rakt iou eines kleinen Schnerkeugeliäiises mit einem Detektor. Die 
Objektkoatur ist nach der FBP in den Objektraumschichten im Querschnitt zu erkennen. 
Über einen Schwellwert lässt sich daraus die Isooberfläche reudern. Durch die schichtweise 
Normierung sind die Haltestäbe des Objektbettes um das Objekt herum als Artefakte zu 
seilen. 
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Abb. 2.26. Das Mausphantom, für das ebenfalls ein einzelner Detektor ausreicht, ist unge- 
eignet für die optische Oberflächenextraktion, da die unterschiedliche Oberflächenfärbung 
zu einer ungleichmäßigen Lichtreflexion führt. Die helle Unterseite dominiert die Silhouette, 
sodass keine passende Isooberfläche gefunden werden kann. 



lässt. durch den die Objektkontur passend segmentiert werden kann. 



Die Oberfläehenextraktion aus zusammengesetzten Detokturprojektionen ist in Abb. 2.27 
demonstriert. Der Würfel hat eine Kantenlänge von 20 nun und überschreitet, in der 
Diagonalen damit das FOV eines einzelnen Detektors. Durch das Zusammenfügen sind in 
den Projektionen bzw. Sinogrammeu sichtbare Kaulen in der Mitte entstanden, die nach 



2.4. Experimenteller Scanneraufbau mit mehreren Detektoren 



55 




0[rad] tt/2 3tt/4 z[px]: 40 80 



■s, 




z[P*l 40 80 Isooberflache Phantom 

Abb. 2.27. Oberflächeuextraktioii aus zusammengesetzten I )etektorprojektionoii. Die Ab- 
messung des Papicrwiirfels können niclit mehr mit einem einzelnen Detektor abgedeckt werden. 
Der Unterschied der hier dargestellten Sonogramme zu einer Transmissionsini •ssung Lsl im 
Vergleich ZU den Simulationsdaten in Abb. 2.28 erkennbar. Neben den Würfelecken sind auch 
die Haltest ähe des Objektbettes deutlich zu Sehen; die Würfelseiten hingegen lassen sich nur 

undeutlich abgrenzen. 

der FBP zu Artefakten in den Objektraumschichten innerhalb der Würfelkontur führen. 
Solange diese Kanten die Kontur nicht durchbrechen, stören sie die Segmentierung nicht. 
Auffallend ist, dass die scharfkantigen Ecken des Würfels deutlich in den transversalen 
Objektraumseliichten hervortreten, wohingegen die glatten Seilenflächen nur schlecht 
abgrenzbar sind. Der Grund hierfür ist, dass die Wütfelnberfläehe in jede Winkclpro- 
jektion nahezu gleichmäßig hell reflektiert. Durch diesen Umstand unterscheiden sich 
die resultierenden optischen Reflexionsprojektionen deutlich von denen einer Transmissi- 
onsmessung. Weiterhin sind die unteren Seitenflächen teilweise durch die Haltestäbe des 
Objektbettes verdeckt und dadurch für die Segmentierung beeinträchtigt. Zumindest für 
den oberen Teil des Würfels ist eine Oberflächenextraktion mittels einer Isooberfläche 
möglich, die das Objekt erkennen lässt. Trotz der erschwerten Abgrenzung der Kontur 
bleibt zu erwähnen, dass die sichtbaren Reken die Quersehnittsgeoiiiel rie des Würfels 
maßstabsgetreu wiedergeben. Zum Vergleich ist in Abb. 2.28 eine simulierte Transmissi- 
onsmessung durch einen dichten Würfel gleicher Größe gezeigt. Die Projektionsdaten 
sind als Sinograunn dargestellt und die daraus resultierende Objektraumschiclit nach 
der FBP abgebildet. Die Studie ist mit denselben Parametern wie für die optische 
Oberflächenextraktion durchgeführt. 

Im nächsten Schritt können mit der bekannten Isooberfläche die notwendigen Ab- 
bildungstiefen für die Signalmessung berechnet werden. Abb. 2.29 zeigt ein größeres 
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Abb. 2.28. Simulation einer Transmissionsmessung durch einen gleichmäßig dichten Würfel 
(Dichte 1) mit Haltestäben (Dichte 0.75). Der Originalschicht (a) ist ein mittelwertfreier 
Rauschauteil (er - 0.5) hinzugefügt worden. Das Siiiogramm (l>) hestehl aus 153 Projektionen 
und die FBP (c) wurde mit denselben Parametern berechnet, wie sie für das Ergebnis 
in Abb. 2.27 verwendet wurden. Der Würfelquerschnitt ist hier erwartungsgemäß besser 
abgrenz bar als mit der optischen Reflexionsprojektionen. 

Schneckengehäuse, von dem die Oberfläche bestimmt wurde Diese Oberfläche lässt 
sieh nun in alle DKSe transformieren, wo dann für jede Winkelstellung eines Detektors 
eine Tiefenkarte des Sichtbaren Bereichs dieser Oberfläche erstellt werden kann. In der 
Abbildung sind vier Tiefenkarte des Objektes beispielhaft dargestellt, die die Abstände 
/iir ( >borflä« he aus demselben Absolul w inkel /< -ig« -u. Kür die Leiden Signal. I. -lekt «an 
(Det. 1 und 4) ist das Objekt zentriert, wogegen die Oberflächeudetektoren (Det. 2 und 3) 
lediglich einen Teil des Objektes erfassen. Eine Tiefenkarte enthält, dann die genauen 
Abbildungstiefen, d.h. den senkrechten Abstand von der Sensor- zur Objektoberfläche. 
in Bezug zu einer Referenztiefe t HZ , die den tatsächlichen Abstand der Sensoroberfläche 
eines Detektor zur Rotationsachse angibt, (s. Abb. A.3). Mit diesen Tiefenkarten sind 
schließlich die erforderlichen Abbildungstiefen determiniert, die für eine fokussierte 
Rücfcwfirtsabbildung der zu messenden Signalintensitäten auf der Objektoberfläche 
notwendig sind. Für die Anwendung des Abbildungsmodells ist es nun nicht mehr 
notwendig, dass die Objektbreite durch einen Detektor abgedeckt wird. Vielmehr ist für 
die Signalmessung mit den jetzt bekannten Abbildungstiefen die zentrierte Anordnung 
der Signaldetektoren von Vorteil. Die hier vorgestellte Methode der optischen Oberflä- 
chenextraktion ist im Ergebnisteil in Abschn. 3.2 eingehend mil einer Labormaus und 
verschiedenen Zylindoi phantomen untersucht worden. Km die segmentierte Oberfläche 
zu validieren, sind einige Studien mit kombinierter CT-Bildgebung durchgeführt worden. 

2.4.4 Signalmessung 

Die Zielstellung dieser Arbeit ist eine möglichst akkurate Signalmessung auf der Ob jekt- 
oberfläche. Typischerweise sind diese Oberflächensignale relativ schwach und es gibt 
keinerlei Information über die Tiefe, aus der die vom Detektor gemessene Intensität 
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Abb. 2.29. Tiefenkarten der vier Detektoren aus derselben absoluten Winkelstellung. Die 
Oberflächendetektoren erfassen durch ihre spezielle Anordnung nur kleinere Bereiche des 
Objektes, während für die Sigiialdetcktoreu ddS Objekt zentriert ist. Die min bekannten 
AbbildungBtiefen ermöglichen eine fokusaierte Rückwärt sabbildung der Signalmessung direkt 
auf die Objektoberfläche. Die Abbildungstiefen sind entsprechend der Skala farbkodiert. 

sinmnit. Mit der im vorigen Abschnitt gezeigten optischen Obernaeheuextraktion 
lassen sich jedoch die Abbildungstiefen für eine Objektoberfläche in jeder beliebigen 
Winkelstellung eines Detektors bestimmen. Unter der Annahme, dass das Objekt 
zwischenzeitlich nicht bewegt wurde, können dann damit die Oberflächenintensitäten 
aus einer anschließenden Signalmessung eindeutig berechnet werden. Ein kompletter 
Messdurchlauf gliedert sich somit in zwei Abschnitte: 

1. Oberflächenbestiinmung des beleuchteten Objektes 

2. Signalmessung ohne Objoktbeleuchtung 

Der Ablauf der eigentlichen Signalmessung hängt in erster Linie davon ab, wie das 
Signal anschließend tomographisch rekonstruiert werden soll und welche Stärke da.s 
Signal aufweist. Im Ergebnisteil in Abschn. 3.2.3 sind selbstleuchtende und extern 
angeregte Signalquellen mit Zy linder phantomen gemessen worden. Als effizient erwiesen 
sich hier die zentriert ausgerichteten Signaldetektoren. Zwar können die Oberllächende- 
tektoren grundsätzlich ebenfalls für die Signalmessung eingesetzt werden, jedoch liefern 
sie weniger Signalinformation, da sie nur einen kleineren Teil der Oberfläche abdecken. 11 
Ein direkter Vergleich der Detektoren für die Signalmessung ist im Ergebnisteil in 
Abschn. 3.2.3.1 durchgeführt worden. 

Die Signaldaten lassen sich sowohl in den einzelnen Signalprojektion als auch summiert, 
auf der Oberfläche auswerten. Eine Signalprvjektion (G4 x 128 px) zeigt dann die Signal- 
intensitäten, die aus den gemessenen Sensordaten aus der jeweiligen Winkelstellung 
direkt auf die zuvor ermittelte Objektoberl lache rückabgebildet wurde. Für ausreichend 

11 Anders als bei der Oberflächenextraktion, in der lediglich die Objektkonturen deutlich 
abzubilden waren, ist das Zusammenfügen zweier Detektorprojektioneu mit variablen 
Abbildungstiefen innerhalb einer Projektion nicht mehr ohne Weiteres möglich. 
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Abb. 2.30. Ablaufprotokoll für die Obernächensiguahnessung. Die komplette Messung kann 
unabhängig von der anschließenden Datenauswertung erfolgen. Die Basis für eine genaue 
Signalberechnung ist eine möglichst akkurate Oberflächenbestimmung. Der Ablauf der Signal- 
uiessung selbst richtet sieh nach der Anwendung und der Signalstärke. 

intensive Signalqucllcn kann eine relativ schnelle Signalmcssung aus nur einer Winkel- 
stellung durchgeführt werden. Für schwächere Signale erfolgt die Signalmessung im 
Ergebnisteil über 153 Winkelstellungen wie schon für die Oberflächexibestimnning. Eine 
typische Aufnahmezeit liegt bei 3 s pro Winkelstellung. Somit erhöht sich folglich die 
Messzeit, aber gleichzeitig auch die gemessene Signalinformation. 



Die Aufnahme der Oberflächen- und Sigiialdatensätze kann unmittelbar hintereinander 
erfolgen. Aus beiden Datensätzen zusammen lassen sich dann nachträglich die Oberllä- 
chensignale berechnen und auswerten. Abb. 2.30 stellt ein typisches Messprotokoll und 
den Ablauf der Datenauswertung dar. Der Oberjlächendatenmtz enthält die segmentierte 
Objektoberfläche im HKS von Detektor 0. Die erforderlichen Reflexionsprojektionen 
werden ausschließlich über die Oberflächendetektoren aufgenommen. Nach dem Zusam- 
menfügen und der FBP wird die Oberfläche hierfür mittels Schwellwert festgelegt. Mit 
der bekannten Oberfläche lassen sich individuell für alle Detektoren und deren Winkel- 
Stellungen Tiefenkarten berechnen, mit denen die Rüekwärtsabbildung nun eindeutig 
bestimmt ist. Ein Signaldatensatz umfasst dann alle Signalprojektionen eines Detektors. 
Typische Messdauern für Oberflächen- bzw. Signaldatensätze betragen etwa 2 min bzw. 
7 min. 
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Ergebnisse 



In diesem Kapitel sind die experimentellen Resultate zusammeugcfasst, die mit den 
zuvor präsentierten Methoden und Aufbauten erzielt wurden. Das Kapitel gliedert 
sich dabei in zwei Teile: In Absehn. 3.1 wird zuerst die Bildgebung mit einem Ein- 
zeldetektor behandelt und die Abbildungseigenschaften charakterisiert. Anschließend 
widmet sich Absehn. 3.2 der Bildgebung mit denn Scanneraufbau und der eigentlichem 
Oberflächensignalmessung. Innerhalb dieser beiden Abschnitte sind die Ergebnisse den 
verschiedenen Phantomstudien zugeordnet, mit denen jeweils spezielle Fragestellungen 
untersucht worden sind. 

3.1 Bildgebung mit einem Detektor 

Der spätere Scanneraufbau besteht aus vier Einzeldetektoren, die nach demselben 
Prinzip funktionieren: Aus dem dreidimensionalen Objektraum wird eine Intensitätsver- 
teilung auf ein zweidimensionales Sensorbild reduziert, in dein über die individuellen 
Mikrolinsenprojektionen die Tiefeninfonnation inhärent enthalten ist. Nachfolgend 
werden die Abbildlingseigenschaften eines einzelnen optischen Detektors experimentell 
untersucht. 

Die Funktionsweise des Abbildungsniodells bei unbekannten Abbildungstiefen wird 
anhand der dreidimensionalen Oberfläche eines Mausphantoms in Unterabschn. 3.1.1 
gezeigt. ( dei< hzeitig stelll sieh hier dh Frage, "I, eine einzelne Sensoraufnahtiie für die 
adäquat e Tiefenrekoustruktion ausreicht. 

Ist die Abbildungstiefe eines Sensurdatensatzes bekannt, so lassen sieh damit unmit- 
telbar fokussierte Obje'ktbilder berechnen. Für diesen Fall sind die Abbildungseigen- 
schaften eines Detektors mit Blick auf die spätere Oberflächensignalmessung im Scan- 
neraufbau bestimmt. Anhand einer Testschablone (Unterabschn. 3.1.2) lässt sich die 
räumliche Auflösung für kleine Abbildungstiefen (Snim - 20 mm) untersuchen und der 
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Effekt der iterativen Anwendung des Abbildungsmodells beurteilen. Für die gesamte 
Tiefe (10 mm — CO mm) des Scannerobjekt raumes eignet sich eine Liiiienc)uelle (l'nter- 
abschn. 3.1.3) besser, mit der sieh die Tiefenabhängigkeit der räumlichen Auflösung 
bestimmen und die korrekte Tiefenkalibrierung eines Detektors prüfen lässt. 

Für die Signalmessung schwacher Lichtquellen wird das Signal-Rausch- Verhältnis des 
optischen Detektors in Abschn. 3.1.4 beurteilt. Anhand einer Lumineszenzquelle wird 
die Signalqualität bei variierender Abbildungstiefe und Belichtungszeit untersucht, die 
mii dem Abbildungsmodell bzw. Fokuspunkt bildern erzieh werden kann. 

3.1.1 Mausphantom 

Eine typische Messanordnung mit einem Einzeldetektor ist in Abb. 3.1 gezeigt. Unmit- 
telbar vor der Detektorfläche ist ein Mausphantoni positioniert, das aus unterschiedlichen 
Richtungen angeleuchtet wird. Das Mausphantoni, etwas kleiner als eine echte Maus, 
bietet eine dreidimensionale Oberfläche mit unterschiedlichen Reflexionseigenschafieu 
und Strukturell, anhand der sich die Bildgebung mit dem optischen Detektor gilt stu- 
dieren lässt. Die resultierende Sensoraufnahme des Detektors ist zusammen mit einem 
Nonnierungsbild in Abb. 3.2 dargestellt. Diese beiden Datensätze sind die Grundlage 
für die Ergebnisse in den folgenden drei Unterabschnitten. 



Abb. 3.1. Messanordnung mit einem optischen Detektor und einem Mausphantoni, das schräg 
vordem Detektor positioniert und angeleuchtet ist. Iiier Ist noch d;is vorherige Rahineiikonzept 
zu sehen, bei dem Sopteninaske und MLA über Justierschrauben gehalten und kalibriert wurden 
(vgl. Abschn. 2.1.2). Die Abbildung ist entnommen aus [13]. 

3.1.1.1 Digitale Refokussierung des Objektraumes 

Ohne Kenntnis der Ob jektoberfläche sind die Abbildungstiefen nicht eindeutig bestimmt. 
Der Objektraum lässt sich aber aus einer einzelnen Sensoraufnahme mit dem Abbil- 
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Abb. 3.2. Sensui aufnahmen (512 x 1000 px): (a) Normierungsbild. aufgenommen mit einer 
homogenen Leuchtfolie, und (b) Mauspltantom aus Abb. 3.1. In den vergrößerten Ausschnitten 
sind die einzelnen Mikrolinsenprojektioiien zu erkennen, die jeweils die individuelle Sicht in 
den Objektraum wiedergeben. Die Abbildung ist entnommen aus [13]. 

dungsmodell sukzessive durchfokussieren. Wie in Abschn. 2.2.5 beschrieben, wird clor 
Objektfokus bzw. die Abbildungstiefe während der Rückwärtsabbildung schrittweise 
variiert, die entsprechenden Objektbilder erzeugt und zu einem Objektdatensatz zusam- 
mongefasst. Die Berechnung eines solchen Datensatzes, der hier aus CO Objektbildern 
mit jeweils 512 x 1000 px besteht, benötigt auf einem Standard-PC (2,13 GHz Intel Gore 
2 Duo Prozessor; 3,25 GB RAM) ungefähr 25 min. Die Abbildungstiefen liegen dabei in 
einem Bereich von 8mm — 38mm und haben konstante Abstände von Ai - 0.5 mm. 
In Abb. 3.3 sind exemplarisch drei Objektbilder in unterschiedlichen Abbildungstiefen 
({11: 18: 28} nun) aus diesem Datensatz gezeigt. Zum direkten Vergleich ist ein Fokus- 
punktbild in der linken oberen Ecke dargestellt, das aus demselben Datensatz extrahiert 
wurde (vlg. auch Abb. 1.2). Der Unterschied in der räumlichen Auflösung ist deutlich 
sichtbar. Weiter lässt sich beobachten, dass die in diesem Sensordatensatz enthaltene 
Tiefeninforutation in den Objektbildern nach der Rückwärtsabbildung wiederzufinden 
ist, und zwar in den tiefenabhängig ((»kassierten Oberflächenstrukturen. Diese Aussage 
war bereits mit der Modellstudie dreier Punktquellen in Abb. 2.12 möglich und zeigt, 
sich nun auch im Experiment. 

3.1.1.2 Tiefenbestimmung einer komplexen Oberfläche aus einer einzelnen 
Sensoraufnahme 

Um die Abbildungstiefen einer einzelnen Sensoraufnahme zu rekonstruieren, sind in 
Abschn. 2.2.6 zwei Ansätze vorgeschlagen, mit denen der Objektdatensatz analysiert 
werden kann. Ziel dabei ist es, die fokussierten Oberflächenbereiche zu identifizieren 
und diese einer Abbildungstiefe zuzuordnen. Hier soll experimentell anhand des Maus- 
phantoms untersucht werden, ob sich dies*? Methoden für eine dreidimensionale Oberllä- 
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Abb. 3.3. Objektbilder (512 x 1000 px) in verschiedenen Abbildungstiefen, berechnet mittels 
Rückwärtsabbildung aus einer einzelnen Sensoraufnahme (Abb. 3.2). Der Objektfokus variiert 
entsprechend und liegt einmal im Kopf- (a). Torso- (b) und Hinterteilbereich (c). Zum Vergleich 
ist in der oberen linken Ecke ein Fokuspiinktbild (, r >l X 10t) px) mit deutlich niedrigerer 
rtUimlicfaer Auflösung vergrößert dargestellt. Der komplette Datensatz ist. iils Mediendatei zu 
Referenz [13] verfügbar, aus der auch diese Abbildung entnommen ist. 

chenbestimmung eignen. 

Die Suche nach lokalen Voxelmaxima ist zur Tiefenbestimmung für isolierte Punktquellen 
durchaus anwendbar, wie es sich in der Modellstudie in Abb. 2.12 zeigt. Der Objeklda- 
tcnsaiz des Mausphantmus ist mit diesem Ansatz analysiert worden und beispielhaft 
für eine x-Sehieht in Abb. 3.4 dargestellt. Exemplarisch sind hier vier Prolildiagranime 
zu sehen, die den Intensitätsverlauf der Voxul entlang der Abbildungstiefe wiedergeben. 
Zwei gut identili zier bare Maxinm (1; 3) lassen sich in dieser Schicht linden. Werden mit 
diesen Abbildungstiefen dann die korrespondierenden Objektbilder dargestellt, so zeigen 
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Abb. 3.4. Ticfcnhcstüiuiuing über lokn.|<> Intensitätsniaxima: In einer x-Schicht wird entlang 
der Abbildungstiefe nach dem Inleusiuüsinaxhnuni in den Objekt vuxeln gesucht. Diese 
Methode führt nur in wenigen Fidlen zum Erfolg (1: 3), der überwiegende Teil der Oberflä- 
chenintensitäten (z.B. 2; 4) lässt sich nicht eindeutig einer Abbildungstiefe zuordnen. 

sich die Oberflächenbereiche um diese Funkte tatsächlich fokussiert. Anders verhält es 
sich mit dem überwiegenden Teil der Voxeliutonsi täten dieser x-Schicht (z.B. 2: 1). der 
nicht mehr so eindeutig einer Abbildungstiefe zugeordnet werden kann. Festzuhalten 
bleibt, dass helle Stellen mit ..Punktquellen-Charakter" auf der Oberfläche in diesem 
Objektdatensatz in der Tat ihrer ursprünglichen Abbildungstiefe zugeordnet werden 
können. I ber markante Punkte, die zudem noch genügend weit voneinander entfernt sein 
müssen, lässt sich also eine grobe Tiefenbestimmung mit dieser Methode durchführen. 

Mit dem zweiten Ansatz wird versucht, fokussierte Oberflächenbereiche in einem Ob- 
jektbild über ein maximales Varianzmaß zu finden. Diese Methode basiert auf der 
Annahme, dass fokussierte Oberflächenstrukturen einen maximalen Kontrast bieten. 
Abb. 3.5 zeigt zwei Tiefenkarten, die mit unterschiedlichen Blockgrüße (20 x 20 px: 
40 x 40 px) aus dem Objektdatensatz berechnet und mit einem Medianfilter in gleicher 
Blockgrüße geglättet wurden. Diese Tiefenkarten bieten schon umfangreichere Tiefenin- 
formationen als der erste Ansatz, sind aber von einer akkuraten Oberflächenbestimmung 
noch deutlich entfernt, da sich einige lückenhafte Stellen und Tiefensprünge ergeben. 
Probleme bereiten diesem Verfahren kontrastreiche Bildartefakte in den defokussierten 
Objekt bereichern Da diese Strukturen besonders in kleinen Abbildungstiefen vorhanden 
sind, beginnt die Tiefenanalyse des Objektdatensatzes erst ab 13 min. Ebenso verhalt es 
sich an den Objekt rändern, die hier nicht eindeut ig bestimmbar sind. Der Hintergrund 
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Abb. 3.5. Tiefenbesthninung über die Intensitätsvariauz der Objektvoxel: Resultierende Tie- 
fenkarten aus dem Objektdateusatz mit unterschiedlicher Blockgröße (20 X 20 px; 40 x 40 px). 
Für jedes Voxcl in einem Objektbild des Olijektdatensatzes wird die Voxclvarianz in einer zuvor 
bestimmten BlockuiURobting berechnet. Über alle Objektbilder wird dann entlang der z- Achse 
der Maximalwert ermittelt und darüber die Abbildungstiefe lest gestellt . Die Tiefenschritte 
der Karte sind durch den Abstand der Objektbilder (Af — 0.5min) vorgegeben. 

könnte noch mittels geeigneter Segmentierung entfernt werden. 

Als vorgezogenes Fazit zur Tiefenbestimmung einer Mausoberfläche aus einer einzelnen 
Sensoraufnahme sei bemerkt, dass grob die Abbildungstiefen unter bestimmten Bedin- 
gungen durchaus gefunden werden können. Allerdings sind die erzielten Ergebnisse mit 
beiden Analysemethoden noch zu ungenau, als dass sie für die exakte Obcrflächonsignal- 
messung eingesetzt werden könnten. Ein interessanter Aspekt für die Optimierung des 
ersten Ansatzes wäre, ein spezielles Louchtpunktmuster auf die Objektoberfläche zu 
projizieren. Dann könnte man versuchen, die Abbildungstiefen dieser Punktquellen zu 
linden. Mit der Varianzmethode werden die Abbildungstiefen über die Oberflächen- 
Struktur des Objektes bestimmt. Grundvoraussetzung hierfür Ist, dass diese Strukturen 
in den Objektbildern hinreichend aufgelöst werden können. Als Nachteil erweist sich 
dabei das abnehmende Auflösungsvermögen in den Objektbildern bei ansteigender 
Abbildungstiefe. Somit wird die Tiefenbestimmuug aus einer einzelnen Sensoraufnahme 
schwieriger, je weiter ein Untersuchungsobjekt entfernt ist. Hinzu kommt, dass die 
erzielbare Tiofenauflösung nochmals geringer ist als die räumliche Auflösung innerhalb 
eines Objektbildes. Diese Umstände bereiten bereits bei dem hier analysierten Ob- 
jektdateusatz mit einem Tiefenbereich von 8 nun — 38 mm Probleme. Deshalb ist die 
Tiefeiibestiiniining uns lediglich einer Sensoraufnahme für den Scannerobjektraiiin nicht 
geeignet, da dieser sich las zu einer Tiefe von 60 min erstreckt. 

3.1.1.3 Iterative Anwendung des Abbildungsinodells auf ein Objektbild 

Die iterative Anwendung des Abbildungsinodells nach Absehn. 2.2.7 kann unmittelbar 
für die Bildberechnung des Mausphantoms eingesetzt werden. Abb. 3.6 zeigt zwei 
Objektbilder in gleicher Abbildungstiefe (t — 18 mm), einmal nach einfacher Rück- 
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(a) (b) 

Abb. 3.6. Objektbild (/ = 18 jiiiii) nach einfacher Rückw&ftS&bbildung (a) und nach iterativer 
Anwendung (-1 Iterationen) des Abbildungsmodells (b). 



wärtsabbildung (a) und nach iterat iver Anwendung des Abbildungsrnodells (b), Nach 
vier Iterationen sind die Oberflächenstrukturen des Phantoms tatsächlich deutlicher 
zu erkennen. Vorteile für die Tiefenbestimmung aus einer einzelnen Sensoraufnahnie 
ergeben sich jedoch nicht aus dieser Methode, da unweigerlich auch die kontrastreichen 
Strukturen der defokussierten Objektbereiche verstärkt werden. 



3.1.2 Testschablone 

Die planare Testschablone aus Abb. 2.15 kann, anders als die dreidimensionale Ober- 
fläche des Mausphantoms. in genau bekannten Abständen vor dein Detektor positioniert 
werden. Damit lässt sieh dann die Abbildung einzelner Raumfrequenzen in einem Ob- 
jektbild in Abhängigkeit der Abbildungstiefe untersuchen. Da die Testschablone ebenso 
als digitale Vorlage verwendet werden kann, ist für jede Detektorinessung zum Vergleich 
eine Modellstudie möglich. Beide Ergebnisse lassen sich dann noch den entsprechenden 
Fokuspunktbildern gegenüberstellen. Das tiefenabhängige Auflösungsvermögen ver- 
schiedener Raumfrequenzen und ebenso die iterative Anwendung des Abbildungsmodells 
ist Gegenstand der zwei folgenden Unterabschnitte. 

3.1.2.1 Räumliche Auflösung in kleinen Abbildungsticfcn 

Die Linienpaare der Testschablone sind so gewählt, dass sie sich nicht mehr in einem 
Fbkuspunktbild darstellen lassen, sondern nur noch in den Objektbilden] des Abbil- 
dungsniodells. Anhand dieser Rauinf'requenzen ist dann abschätzbar, bis zu welcher 
Abbildungstiefe mit der Rückwärtsabbildung feine Strukt uren noch aufgelöst werden 
können [48]. Abb. 3.7 zeigt das Ergebnis der Linienpaannessung. Von der Testschablone 
sind hier für drei Abbildungstiefen / ({8: 11: 20} mm) künstliche und gemessene Sonsor- 
datensätze erzeugt und daraus die Ob jektbilder für die entsprechenden Abbildungstiefen 
berechnet worden. Aus der Detektorinessung wurden zusätzlich die Fokuspunktbilder 
extrahiert. Vergleicht man die Resultate in einer Tiefe, so zeigt sich die beste Auflösung 
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stets für die rauschfreien Sensordatensat /e der Modollstudie. Per Augenmaß lassen sich 
folgende gerade noch unterscheidbare Linienmuster erkennen: {3.5; 2,1: l t 5} 'p/mm in der 
Modellstudie und {2.6; 2.1: 1,3} 'i'/mm in der Detektormessung. Im Gegensatz dazu ist in 
den Fokuspunktbildern bestenfalls das Mustor mit 1 'p/mm erkennbar, das aufgrund der 
Grenzfrequenz eines Fokuspunktbildes gerade noch aliasingfrei abgetastet werden kann. 
Speziell in kleinen Abbildungstiefen sind die Vorteile des Abbildungsmodells gegenüber 
der Fokuspuiiktbilder somit sehr deutlich. 
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Abb. 3.7. Räumliche Auflösung bei ansteigender Abbildungstiefe in den Objektbildern, 

berechnet uns künstlichen (a) und gemessenen (h) Sensordaten. Die Fokuspunkthikler (c) sind 
aus der Detektormessung extrahiert und zum Vergleich entsprechend vergrößert dargestellt. 
Die Abbildung ist entnommen aus [48]. 




Die tiefenabhängige Abnahme der räumlichen Auflösung in den Objektbildern ist bei- 
spielhaft in Abb. 3.8 für ein Linienpaarmusters (lj'p/mm) der Testschablone gezeigt. 
In eleu Prolildiagrannnen ist zu erkennen, dass sich mit ansteigender Abbildungstiefe 
zunehmend der Rechteckkontrast verringert. Da der Rechteckkontrast über Korrek- 
turtenne mit der MTF verbunden ist. kann aus der Kontrastverringerung auf eine 
Abnahme der räumlichen Auflösung geschlossen werden. 

Im nächsten Abschnitt wird über die iterative Anwendung des Abbildungsmodells 
der Kontrast in den Objektbildern wieder verbessert. Das Vorgehen wurde 1 hierfür 
bereits in Abb, 2.14 demonstriert. Der Kontrast eines Redlteckmusters wird dafür 
als Differenz zwischen dem Mittelwert der drei Maxime und dem Mittelwert der zwei 
Minima berechnet. Für die Betrachtung der Kontraständerung wird der Kontrast 
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in einem Objektbild (Ausgangskontrast) ins Verhältnis zum Kontrast der Vorlage 
(Eingaiigskont rast.) gesetzt . 




Abb. 3.8. Kontraste hnahnie eines Rechteckmusters (1,7 'p/mm) der Tcstschnblonc in unter- 
schiedlichen Abbild u lißstiefen t ({8: 11; 20} mm): LinienproHl des Rechteckmusters in der 
Schablone (a) und der entsprechende Kontrast im Objektbild der Modellstndie (b) und der 
Detektonnessung (c). 



3.1.2.2 Iterative Anwendung des Abbildungsuiodells in kleinen 
Abbildungstiefen 

Anhand der Testschablone lässt sich die Kontrastvorbossorung dnreh die iterative 
Anwendung des Abbildungsmodells quantitativ beurteilen. Dafür ist die Kontrast Ver- 
änderung verschiedener Liiiienpaarinuster der Testschablone während der einzelnen 
Iterationen in zwei Abbildungstiefen / ({ 14: 20} tum) ermittelt und in Abb. :3.9 gezeigt. 
Auffallend ist. dass für die künstlichen Sensordaten der Ausgangskontrast den Ein- 
gangskontrast teilweise übersteigt. Die Ursache liegt hier an Kantenübersehwingern, die 
diese Methode speziell an ohnehin schon kontrastreichen Strukturen erzeugt. Erwar- 
tungsgemäß schwächer stellt sich die Kontrastverbesserimg bei den Detektormessdaten 
dar. Grundsätzlich gilt in beiden Fällen, je weniger Ausgangskontrast nach der einfachen 
Rückwärtsabbildung (Iteration 0) vorhanden ist. desto flacher verläuft die Verbesserung. 
Nach fünf Iterationen sind schließlieh auch bei anfangs steilen Kontraständerungen 
nur noch gering«' Verbesserungen möglich. Die positive Auswirkung der iterativen An- 
wendung des Abbildungsmodells zeigt sich eher bei kleinen Abbildungstiefen. Da der 
Effekt der Auflösungsverbesserung bei größeren Abbildungstiefen an Einfluss verliert, 
lässt diese Methode keinen Vorteil für die Bildgebung im Scanneraufbau erwarten. 

3.1.3 Linienquelle 

Um den Zusammenhang zwischen Abbildungstiefe und Auflösungsveränderung in einem 
Objektbild auch über einen größeren Bereich zu untersuchen, eignet sich die Linien- 
bildfunktion (s. Abschn. 2.3.1) besser als die Testschablone. Wenn man also weiß, wie 
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Abb. 3.9. Koutrastäuderiing der Linienpaarmuster (Beschriftung der Verläufe in ('n/mm]) 
nach iterativer Anwendung des Abbildungsmodells in zwei Abbildungsticfen t ({14; 20} nun) 
jeweils für künstliche (a) und gemessene (b) Sensordaten. Die Kontraständcrung bezieht sich 
auf den Eingangskontrast der Schablone. 



sich das Linionbild in Abhängigkeit der Tiefe verbreitert, dann kann daraus auf das 
tiefenabhäugige Auflösungsvermögen geschlossen werden. Uber die Profilbreite der Lini- 
eubiklfuiiktion lässl sich so in Modells! udien ein funktioneller Zusammenhang zwischen 
der raumlichen Auflösung und der Abbildungstiefe herstellen. Weiterhin können auch 
Einflüsse der einzelnen Modellparameter auf die räumliche Auflösung untersucht werden. 
Die Ergebnisse lassen sich dann mit einer Detektormessung unter gleichen Bedingungen 
überprüfen. In den nachfolgenden zwei Unterabschnitten werden Linienquellen für die 
Auflösungsmessnng und die Tiefenkalibrierung eingesetzt. 



Die Prolilbteite einer Linienbildfuuktion wird dabei wie in Abschn. 2.3.1.2 ermittelt. 
Für die Ergebnisdarstellung in den nachfolgenden Diagrammen sind ungefüllte Symbole 
für Modellstudien und gefüllte Symbole für die Detektormessung gewählt. Für die 
Detektormessung ist die Proiiibreite an mehreren Stellen im Objekt bild bestimmt und 
als Mittelwert mit der Standardabweichung als Fehler im Diagramm eingetragen. 
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Abb. 3.10. Modellstudie zur tiefenabhängigen Linienbikibreite (Symbole) bei variierender 
Fokuslänge. Der lineare Zusammenhang zwischen Abbildungstiefe / und Halbwertsbreile des 
Liuienbildprofils lä.ssl sieh näheruugsweise durch Gl. 3.1 beschreiben (durchgezogene Linien). 

3.1.3.1 Räumliche Auflösung 

Der Einfluss verschiedener Modellparameter auf die räumliche Auflösung lässt sieh 
am Einfachsten über eine Modellstudie untersuchen. Abb. .3.10 zeigt die Halbwerts- 
breiten der Linienbildfunktionen für unterschied liehe Fokuslängen / ({1: 1,5; 2; 2.5; 
3} mm) in ansteigenden Abbildungstiefen. Die Ordinatenskala verläuft in umgekehrter 
Reihenfolge, d.h. je breiter die Linienbildfunktion. desto tieler liegt ein Wert. Der 
Zusaminenhang zwischen der Verbreiterung der Linienbildfunktion bei unterschiedlichen 
Abbildungstiefen / und einer festen Fokuslänge / lässt sich auf Grundlage dieser Studie 
näherungsweise durch eine lineare Funktion beschreiben mit 



Die Halbwertsbreile ist dabei auf den Pixelahstaud <i bezogen. Die durchgezogenen 
Linien im Diagramm zeigen die entsprechenden Verläufe. Iu Abb. 3.11 ist der Einfluss 
weiterer Modellparameter bei fester Fokuslänge / = 2 mm dargestellt mit 

(a) Mikrolinsenabstand d ({O.'Mi: II. 18: 0.2 1} mm) und angepasstem Bohrlochdurchmesser 
D ({0,94: 0,46; 0,22} mm), 

(b) Bohrlochdurchmesser im Verhältnis zum festen Mikrolinsenabstand % (1; 2/3; 1/2; 
V-i) mit d — 0.48 mm. 



FWHM,,„(0 _t-f 



(3.1) 



/ 
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Zwar variieren die Resultate 1 etwas, folgen aber in der Tendenz dem Verlauf nach 
Gl. 3.1. Als Ergebnis dieser Modellstudien lässt sich dann auch für die Parameter des 
optischen Detektors die Verbreiterung der Liniölbildfunktion in guter Näherung durch 
diesen Zusammenhang beschreiben. 
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Abb. 3.11. Modellstudie: Einthiss weiterer Modellparameter auf die tiefcnahhängige Linien- 
bildbreite bei fester Fokuslänge / - 2 nun: (a) Mikrolinsenahstand d und (b) Bohrlochgröße als 
Verhältnis D/d mit d = 0.48 min. Trotz stellemveiser Abweichungen folgt die Halbwertsbreite 
des Linienbildprofils weiterhin dem Verlauf von Gl. 3.1. 



Die Detektormessung der Linienbildbreite gegen die Abbildungstiefe ist in Abb. 3.12 
gezeigt. In bekannten Abständen t (±2 mm) ist eine durchleuchtete Schlitzschablone 
parallel zum Detektor positioniert. Für zwei Schlitzbreiten (100 /nn; 25 //m) erfolgten 
5 bzw. 3 Messdurchgänge. Die Halbwertsbreiten der Linienbildfunktion sind daraus 
entsprechend gemittelt. Die Messwerte verlaufen unterhalb der Ideallinie nach Gl. 3.1 
für den Detektor (gestrichelte Linie). Das bedeutet, dass die Linienbilder der De- 
tektormessung um etwa 0,2 mm breiter sind als theoretisch möglich. Somit ist die 
räumliche Auflösung erwartungsgemäß schlechter als im Idealfall. Der mittlere Fehler 
der Messwerte ist für die kleine Schlitzbreite (25 /an) insbesondere bei kleinen Abbil- 
dungstiofen (f < 15 mm) geringer und bestätigt dort weiterhin den Idealverlauf der 
Halbwertsbreite. In größeren Abbildungstiefen sind beide Schlitzgrößen ausreichend 
kloin gegenüber der Halbwertsbreite der Linienbildfunktion. Das Bestimmthoitsmaß 
wird mit ansteigender Abbildungstiefe besser, da dort die Linienbildprofile glatter 

1 Es sei angemerkt, dass nicht jede hier untersuchte Parameterkombination mit der Paraxialnä- 
herung des Abbildungsmodells vereinbar ist. Bei einem Mikrolinseuabstaud d = 0,96iinu 
und einer Pokuslänge / = 2 nun lässt sieh die MI.A- Abbildung nicht mehr passend durch 
das Abbildungsmodell erklären, da der größere Akzeptanzwinkel einfallender Zentralstrahlen 
die Paraxialnäheruug nicht mehr erlaubt. 
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sind. Insgesamt liegen die Werte weitestgehend über 90% und sprechen für akzeptable 
Schätzwerte. Auffallend sind Ausreißer in bestimmten Abbildungstiefell. In diesen Ab- 
ständen haben Objektbilder eine besonders „pixelige" Struktur aufgrund der diskreten 
Objektrauinberechnung. che speziell in diesen Abbildungstiefen einsieht. Das führt zu 
Linienbildprofilen. an die die Gaussfunktion nur unzureichend angepasst werden kann 
(niedriges Bestimmtheitsmaß, größerer Messfehler) und somit abweichende Funktions- 
breiten liefert. Diese Ausreißer lassen sich auch in Abb. 3.11 finden. 

Abbildungstiefe t [mm] Abbildungstiefe t [mm] 




(a) (b) 



Abb. 3.12. Detektormessung der tiefenabh&ngigea Verbreiterung der Linienbildfunktion. Als 
Linienquelle dienen zwei Schlitzschablonen mit 100 um (a) und 25 /im (b) breiten Schlitzen. 
Oben sind die gemittelteu Messwerte aus mehreren Messdurchläufen im Vergleich zum Verlauf 
von GL .'LI (gestrichelte Linie) dargestellt. Die mittlere Standardabweiehim« der einzelnen 
Messwerte (Mitte) und das mittlere Bestunmtbeitsmaß für die Regressionsrechnung (unten) 
geben die Anpassungsgüte wieder. 

Aus den Halbwertsbreiten aus Abb. 3.12(b) lässt sich jetzt mit der Methode aus 
Abschn. 2.3.1.2 eine MTF für die entsprechende Abbildungstiefe angeben und mit 
einer idealen MTF aus der Modellstudie vergleichen. In Abb. 3.13 sind für verschiedene 
Abbildungstiefen ({ 10; 20; 30; 40; 50; 60} mm) die MTF- Verläufe einmal für den Idealfall 
(a) und einmal für die Detektormessung (b) dargestellt. Durch die reziproke Beziehung 
der Linienbildbreite im Ortsraum und der entsprechenden MTF-Breite im Frequenzraum 
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wird die normierte MTF immer schmaler, gleichbedeutend mit der Verschlechterung der 
räumlichen Auflösung. Der Abstand zwischen den idealen und gemessenen Verläufen ist 
anfangs noch deutlieh, wird aber mit zunehmender Abbildnr.^stiof,. ti; '-n:^ t. Kur die 
Detektormessung ist schon ab einer Abbildungstiefe von t = 10 mm die Übertragung 
von Raumfrequenzen höher als 4 , P/ m m auf null gedrückt, obwohl die Grenzfrequenz des 
Sensors bei f g = 10,4 lp /mm liegt. Im Vergleich dazu könnte mit einem Fokuspunktbild 
(f ;l = 1,04 •p/mm) eine Detektor-MTF aber erst ab einer Abbildungstiefe (t > 50mm) 
aliasiugt'rei abgetastel werden. Das Vbbildungsmodell liefert also über den gesamten 
relevanten Bereich des Scannerobjektraumes eine bessere räumliche Auflösung als ein 
Kokuspun i\ l)ild. 




1 2 3 

Raumfrequenz f [Ip/mm] 

(a) 



1 2 3 

Raumfrequenz f rjp/mm] 

(b) 



Abb. 3.13. Tiefenabhängige, normierte MTF: Gegenüberstellung der Verläufe in verschiedenen 
Abbildungstiefcu (10 nun — 60 nun) im Idealfall (a) und für den optischen Detektor (b). 
Grundlage der dargestellten Funktionen sind die tiefenabhängigen binienbildbreiten nach 
Gl. 3.1 bzw. der Detektormessimg aus Abb. 3.12(b). 



3.1.3.2 Tiefenkalibrierung 

Die Halbwertsbreiten der Linienbilder lassen sich zusätzlich effizient für die Tiefenkali- 
brierung verwenden. Wird die Abbildungstiefe für einen festen Quellenabstand variiert, 
dann lässt sich testen, ob die Linienquelle mit der gegebenen Fokuslänge auch in der 
richtigen Tiefe fokussiert wird (Absehn. 2.3.1.2). Derart kann z.B. die Annahme der 
effektiven Fokuslänge für das MLA überprüft werden. Für Abb. 3.14(a) wurden mit 
der Vorwärtsabbildung (/ = 2,1 nun; t = 20mm bzw. 10 mm) zwei künstliche Sensor- 
datensätze erzeugt. Daraus sind dann die Profilbreiten der Linienbildfunk tion über 
die Rückwärtsabbildung mit drei verschiedenen Fokuslängpn / ({2.1: 2,4; 2,7} mm) 
und in Tiefenschritten von At = 2 mm ermittelt. Beim Defokussieren vergrößert sich 
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die Prolilbreite. und insgesamt gesehen folgt sie in etwa einein parabolischen Verlauf 
(durchgezogene Linie). Lediglich für die richtige Fokuslänge fällt die minimale Linien- 
bildbreite mit der tatsächlichen Abbildungstiefe zusammen. Auf die gleiche Weise ist in 
Abb. 3.14{b) der Detektor für dieselben Fokuslängen getestet. 
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Abb. 3.14. Tiefenkalibrierung der Fokuslünge: In der Modellstudie (a) und für einen ge- 
iness. neu Sensordatensatz (b) worden die Abbildung iefen si hrittweise mit drei unt' rschie I- 
lichen Fokuslängen {2.1: 2,4: 2,7} mm variiert. Nur für die ursprüngliche Fokuslänge fallen 
die tatsächliche Abbildungstiefe (t - 20 nun bzw. 40 mm) der Linienquelle und die minimale 
Linienbildbreite zusammen. 



In Abb. 3.15 ist die Tiefenkalibrierung für / = 2,4 mm in verschiedenen Abbildungs- 
tiefe)) ({ 15; 10: 20; 30; 40; 50: 60} mm) einmal als Modellstudie (a) und einmal mit den 
Detektormessdaten (b) gezeigt. Die Parabeln werden mit ansteigender Abbildungstiefe 
flacher, d. h. dass die fokussierte Linienquelle gar nicht mehr so genau lokalisiert werden 
kann. Dennoch liegen die Miniina der Parabeln immer nahe den entsprechenden Ab- 
bildungstiefen. Somit ist die effektive Fokuslänge der geeignete Parameter für die 
Rückwärtsabbildung, um mit dem Detektor tiefenkalibrierte Objektbilder zu berechnen. 
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Abbildungstiefe t [mm] Abbildungstiefe t [mm] 



10 20 30 40 50 60 10 20 30 40 50 60 




(a) (b) 



Abb. 3.15. T iefcn kalibrierte Detektorabbildung mit der effektiven Fokuslänge. Für ver- 
schiedene Linienquellenabstände wurden Objektbilder einer Modellstudie (a) und einer Detek- 
tormessung (b) mit dein 10(>-/iin-Schhtz ausgewertet. Mit der effektiven Fokuslange / = 2,4 min 
wird die Linienquelle in den richtigen Abbildungstiefen fokussicit. 

3.1.4 Lumineszenzquelle 

Für die Signalmessung mit dem Scanneraufbau werden kleine, selbstleuchtende Quellen 
eingesetzt, die mit einer konstanten Stärke Licht emittieren. Bei diesen relativ schwach 
leuchtenden Luuiineszenzqucllen handelt es sich um kleine, luftdicht verschlossene 
Glaszylinder (0.9mm Durchmesser. 2.5mm Länge), die an der Innenseite mit Phosphor 
beschichtet und mit Tritiumgas (/^-Strahler) befüllt sind [49]. Die Zerfallselektronen des 
Tritiums werden durch den Phosphor absorbiert und in Emissionsphotonen bestimmter 
Energie umgewandelt. Anhand dieser Quellen, die beliebte Signalgeber für Experimente 
zur Biolumineszenztomographie |50| sind, soll die Bildgebung schwacher Signale mit 
dem Abbildtuigsmodell untersucht werden. Sensorsignale liegen bei 3 s Belichtungszeit 
etwa bei 5% des Sättigungswertes. Die Frage ist, ob das Abbildungsmodell das Signal- 
Rausch- Verhältnis des Sensors bzw. Fokuspunktbildes bei schwachen Signalstärkeu 
gemäß Absehn. 2.3.2 verbessern kann. 

3.1.4.1 SNR und Abbildungstiefe 

Für die Modellstudie sind von einer digitalen Vorlage (400 x 400 px) mit einer rechte- 
ckigen, einheitlichen Flache (10 x lÜpx), die die Luminoszcnzquelle nachbilden soll, über 
die Vorwärtsabbildung in verschiedenen Abbildungstiefen ({10; 20; 30; 40: 50; 00} nun) 
künstliche Sensordaten erzeugt worden. Jedem Sensordatensatz wurde dann Rauschen 
mit ft — 0.15 und o — 0.5 hinzugefügt, das ungefähr der Größenordnung in vergleichbaren 
Messdaten entspricht. Abb. 3.16 zeigt Ausschnitte der resultierenden Objekt- und Fo- 
kuspunktbilder aus den Sensordaten. In den Objektbildern nimmt mit ansteigender 
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Abbildungstiefe sowohl das Maxhnalsignal aufgrund der Abbildungsverbreiterung ab. 
zugleich jedoch auch che Rauschvarianz. Im Gegensatz dazu bleibt das Rauschen in den 
Fokuspunktbildern in jeder Abbildungstiefe bei abnehmender Signalstärke gleich. Die 
Detektormessung der Lumineszenzquelle für dieselben Abbildungstiefen ist in Abb. 3.17 
zu sehen. Bei einer Belichtungszeit von 3 s zeigt die Signalhöhe einen vergleichbaren 
Verlauf wie in der Modellstudie, und auch die Rauschvarianz ist in den Objektbildern 
schwächer als in den Fokuspunktbildern. In kurzen Abbildungstiefen sind besonders 
in den Fokuspunktbildern BlendefFekte um die Quelle herum zu erkennen, die auf 
«las unvorteilhaft reflektierende Material der Septenmaske zurückzuführen sind. Diese 
Artefakte nehmen mit ansteigender Abbildungstiefe ab und sind in den Objektbildern 
durch die Rüekwärtsabbildung generell abgeschwächt. 
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Abb. 3.16. Modellstudie: Objekt- und Fokuspunktbilder aus verrauschten Sensordaten in 
ansteigender Abbildun&stief'e. I);is l'rorildiagrninm zeigt in relativem Verhältnis die Signalhöhe 

zur Rausch varianz. Die Fbkuspunktbilder sind zum besseren Vergleich vergrößert dargestellt. 



Die resulti< reu den Verläufe des SM) Im die Modellstudie und die I ). i >kt.< 'lmessung 
sind in Abb. 3.18 dargestellt. Für den theoretischen Fall ist mit dem Abbildungsmodell 
für ansteigende Abbildungstiefen sogar eine Verbesserung möglich, während durch 
die gleichbleibende Rauschvarianz und eine abnehmende Signalintensität das SNR in 
den Fokuspunktbildern langsam abfällt. Die Messdaten zeigen jedoch, dass das SNR 
tendenziell mit ansteigender Abbildungstiefe abnimmt, aber für das Abbildungsmodell 
jederzeit einen besseren Wert liefert, als es mit den Fokuspunktbildern möglich ist. 
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Abb. 3.17. Detektormessung: Objekt- und Fokuspunktbilder aus verrauschten Sensordaten 
in Husteigeuder Abbildungst iefe im Vergleich zu Abb. 3.1(i. Die Signalhöhen zeigen einen 
ähnlichen Verlauf wie in der Model Ist udie. In kleinen Abbilduugstiefen sind hier zusätzlich die 
Störeinflüsse der reflektierenden Septenmaske deutlich in den Fokuspunktbildern zu erkennen. 

Modellstudie Detektormessung 




10 20 30 40 50 10 20 30 40 50 60 

Abbildungstiefe i [mm] Abbildungstiefe t [mm] 

(a) (b) 



Abb. 3.18. SNR bei variierender Abbüdungstiefe: Gegenüberstellung der Modellstudie (a) 
und der Detektonuessung (b). Die Werte sind nach Gl. 2.18 berechnet und in logaritlunischcr 
Skalierung dargestellt. Das SNR in den Objektbildern (Obj.) ist stets besser als in den 
entsprechenden Fokuspunktbildern (Fok.). 

3.1.4.2 SNR und Belichtungszeit 

Bei einer konstanten Abbildungstiefe t — 20 nun der Luinineszenzquelle ist nun die 
liH teht in igszeil nach Möglichkeiten des Sensors wariieri worden. Abb. 3.19 zeigt den 
erwartungsgemäß linearen Anstieg der Signalhöhe und gleichzeitig des Rausclnnittel- 
werte> in den Objekt- und Fokuspunktbildern. Die Rauschvarianz ist aber auch hier 
in den Objektbildern niedriger. Das SNR in Abhängigkeit der Belichtungszeit ist in 
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Abb. 3.20 dargestellt und verbessert sich für das Abbildungsmodell mit zunehmender 
Belichtungszeit leicht. 
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Abb. 3.19. Detektormcssimg: Objekt- und Fokuspunktbilder bei / - 2ömm und zunehmender 
Belichtungszeit. Die Fokuspunktbilder sind wiederum vergrößert dargestellt. 
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Belichtungszeit [s] 

Abb. 3.20. SNR bei variierender Belichtungszeit: Die Werte sind nach Gl. 2.18 berechnet 
und in logarithinischer Skalierung dargestellt. Wie zuvor bei den Abbildungstiefen, ist auch 
hier das SNR in den Objektbildern (Obj.) ausnahmslos besser als in den entsprechenden 
Fokuspunktbildern (Fok.). 



3.2 Bildgebung mit mehreren Detektoren 

Nachdem die Bildgebung eines Einzeldetektors mit dem Abbildungsmodell charakte- 
risiert ist, soll der experimentelle Seannernufbau mit vier Detektoren untersucht werden. 
Der abzubildende Scannerob jektraum, bezogen auf einen einzelnen Detektor, hat eine 
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Tiefe von etwa lOinni - 60mm, für die die Abbildungseigenschaften nun bekannl sind. 
Grundsätzlich muss vor einer Signahnessung die Objektoberfläche bestimmt werden, 
um für die Rüekwärtsabbildung eine eindeutige Abbildungstiefe zu erhalten. Die 01»t- 
flächonext rakt ion mit dem Scanneraufbau ist in Absehn. 2.4.3 detailliert erklärt und 
wird nachfolgend auf alle Untersuchungsobjekte augewendet. Um die Ergebnisse zu 
überprüfen, sind die Studien mit kombinierter OT-Bildgebung durchgeführt worden. 
Gleichzeitig lässt sich hiermit auch der Ansatz für eine multiinodale Bildgebung demons- 
trieren, da der kompakte, optische Scanneraulbau in einem CT-Scanner in Patientengröße 
Platz findet. Die beiden resultierenden Datensätze lassen sich über eine rigide Transfor- 
mation fusionieren. Derart sind nachfolgend die Oberfläche einer nackten Labormaus 
extrahiert und die Lumineszenzquellen auf und in Z>-lin< b 'i ph;tii^ mo -n gemessen wurden. 
In einer zusätzlichen Studie wurde der Scannerautbau um einen optischen Filtersatz 
und einen Anregungshuser erweitert, sodass abschließend die Fluoreszenzbildgebung 
getestet werden kann. 

3.2.1 Kombinierte optische und CT-Bikigcbung 




Abb. 3.21. Experimentalanordnung für die optischen Phantommessungen in einem Compu- 
tertomographen. Der optische Scanneraufbau ist während der kombinierten Messung auf der 
Patientenliege fixiert. 



Für die kombinierte optische und CT-Bildgebung ist der experimentelle Scanneraufbau 
selbst in einen größereu CT-Scanner gestellt worden. Bei dem Gerät handelt es sich um 
einen „Hat-panel Volunie-CT" (fpVCT) Scanner (Siemens Modical Solutions, Forchhehn, 
Deutschland) mit Cone-Beam-Geometrie. der als Prototyp am DKFZ installiert ist [51]. 
Die Rotationsplattform war hierfür auf der Patientenliege während der zuerst durch- 
geführten optischen Aufnahme und der anschließenden Röntgenmessung an gleicher 
Stelle fixiert. Abb. 3.21 zeigt die Vorsuehsaiiordnung. die für die optische Aufnahme 
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komplett, abgedunkelt wurde. Die CT-Messung ist jedes Mal mit einer Rölirenspamiuiig 
von SOkV, einem Anodenstroui von 50 mA und mit einer Ausleserate von 30 fps hei 
zwei vollen Umdrehungen durchgeführt worden. Die Vohunendatensätze sind dann für 
jede Modalität unabhängig voneinander rekonstruiert. Die Voxelgrößc des optischen 
Datensatzes entspricht (0,384)'* nnn' J und die des CT-Datensatzes (0.2)' 1 imu J . Für die 
Registrierung sind lediglich Rotation und Translation notwendig, um die beiden Da- 
tensätze zur Deckung zu bringen. Die Überlagerung ist mit der Software AMIDE [52] 
über Bezugspunkte durchgeführt worden. Abi). 3.22 zeigt das schon zuvor verwendete Ka- 
libriorstäbchon (Abschn. 2.4.2) in drei Schichten der beiden registrierten Datensätze. In 
der CT- Rekonstruktion (Grauwerte) sind die optischen Detektoren sichtbar. Aufgrund 
der Materialeigenschaften verursachen die Detektoren zusätzlich Artefakte im CT- 
Datensatz. die auch in den nachfolgenden Studien nicht vermeidbar sind. Der optische 
Datensatz (Falschfarben) liegt innerhalb des Detektorrings und die Randartefakte deuten 
das vergrößerte FOV an. Das Kalibrierstäbchen deckt sich nach der Registrierung voll- 
ständig in beiden Modalitäten. Später gezeigte überlagerte Datensätze konzentrieren 
sich nur noch auf den Ausschnitt, den die optische Oberflächenextraktion einnimmt. 
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Abb. 3.22. ÜbcrlaRertc Darstellung der CT-Rckonstruktion (Grauwerte) mit dem optischen 
Datensatz (Falschfarhen) des Kalibrierstäbchens. Die Rotationsplattfonn befindet sich in 
der AnfangswinkeLstellung. Das vergrößerte FOV des optischen Datensatzes ist aufgrund der 
Raiidartefakte erkennbar. Da.s Stäbchen deckt sicli in allen drei Raiunan.siehten der beiden 
Datensätze. 



Es sei an dieser Stelle angemerkt, dass die Röntgenstrahlung, die während der CT- 
Messungen appliziert wurde (10 Messdurchgänge), zu einem bleibenden Effekt bei den 
vier Sensoren geführt hat. So war zu Beginn der Studie über die Ausleseelektronik 
ein Dunkelrauschoffset von // = 20 bei 1 s Belichtungszeit eingestellt, der sich nach 
den Messungen unter gleichen Bedingungen verdoppelt hatte. Als logische Konsequenz 
verringert sich dadurch die Empfindlichkeit der optischen Detektoren. Es bleibt zu 
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Isooberfläche Nacktmaus 



Abb. 3.23. Optische Oberflächenextraktion einer Nacktmaus: Dargestellt sind je zwei Projek- 
tionen. Sinogramme und Transversalebenen des Objektraumes. Zusätzlich ist die extrahierte 
Oberfläche im Vergleich zur Nacktmaiis abgebildet. Zum Zeitpunkt der Messung war die Maus 
nicht lebendig und zuvor für die bessere Aufbewahrung tiefgefroren. 

Untersuchen, Wie sich unter diesen Umständen die optischen Detektoren zukünftig mit 
einer Röntgenbildgebung kombinieren lassen. 

3.2.2 Oberfläche einer Nacktinaus 

Der experimentelle Scanneraufbau soll zukünftig in der Kleintiorbildgebung für La- 
bormäuse eingesetzt werden. In dieser St udie ist die optische Oberflächenextraktion 
an einer toten Nacktmaus mit kombinierter CT-Bildgebung durchgeführt. Von der 
beleuchteten Nacktmaus ist ein Oberflächendatensatz aus 153 Projektionen mit jeweils 
0.5 s Belichtungszeit berechnet worden. Die Vorstufen zur Oberflächenextraktion sowie 
die resultierende Isooberfläche sind in Abb. 3.23 gezeigt. Zwar kann die Maus mit dem 
vergrößerten FOV in ihrer kompletten Breite erfasst werden, jedoch ist sie in der Länge 
nicht vollständig durch die optischen Detektoren abgedeckt. 

Die registrierten Datensätze der kombinierten optisc hen und CT-Bildgebung sind in 
Abb. 3.21 gezeigt. Für jede der drei Ansichten ist eine CT-Schicht einmal alleine und 
einmal mit überlagerter Schicht aus dem optischen Datensatz zu sehen. Die optisch 
extrahierte Mausoberfläche ist in jeder Schicht als Binärbild (weiße Bereiche) dargestellt 
und der CT-Datensatz weiterhin in Graustufen. Lokale Abweichungen bis maximal 
2.7 mm zwischen den Oberflächen kommen besonders im Bereich der Haltestäbe des 
Objektbettes auf der Mausunterseite vor. Auf der Oberseite lässt sich dagegen eine 
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Abb. 3.24. Mausoberfläche in kombinierter optischer und CT-Bildgebung: Von jeder Ansicht 
ist zum Vergleich die CT-Schicht (CJraustufen) einmal alleine und einmal überlagert mit der 
Schicht aus dem optischen Datensatz (Binärbilder aufgrund des Schwellwertes, weiß) dargestellt. 
Da die Maus nicht in ihrer gesamten Länge auf die optischen Detektoren pa.sst. ist die optisch 
extrahierte Oberfläche in der Länge beschnitten. In dem gemeinsamen Objektraumbereich 
decken sich die Mausoberflächen der beiden Modalitäten aber gut. 

maximale Oberfl&chendifferenz von etwa 1.9 mm abschätzen. 

Mit der nun bekannten Mausoberfläche kann für jede Detektorposition die Abbil- 
dungstiefe bestimmt werden, wie sie für die Signalmessung mit dem Abbildlingsmodell 
benötigt wird. Dazu wird die Isooberfläche in das entsprechende DKS transformiert und 
der Abstand vom Sensor zu den sichtbaren Oberflächenbereichen berechnet. Dabei ist der 
Abstand 1 uz der Sensoroberfläche zum Rotationszentrum die Bezugsgröße. In Abb. 3.25 
sind die Abbildungstiefen zur Mausoberfläche für verschiedene Winkelstellungen von De- 
tektor 1 gezeigt. Für die Signalmessung ist eine komplette Objektabdeckung nicht mehr 
notwendig. Signalintensitäten können jetzt direkt auf der Mausoberflache berechnet 
werden, falls sich die Mausposition nicht ändert. 

3.2.3 Zylinderphantome 

Anhand von Zylinderphantomen soll nun die Signalmessung mit dem Scanneraufbau 
demonstriert werden. Zylinderphantome eignen sich für die Untersuchung optischer 
tomographischer Verfahren, da sie eine einfache geometrische Form haben und mit spe- 
zifizierten optischen Eigenschaften hergestellt werden können |12]. In den nachfolgenden 
Studien werden die Phantome jedoch verwendet, ohne dass die genauen optischen 
Kigonsehaflen des Materials im Vordergrund stehen. Vielmehr soll die Signalmessung 
für typische Messsituationen getestet werden. Bei zwei Zylinderphantomen werden die 
Lumineszenzquellen als Signalgeber eingesetzt und jeweils auf bzw. unter der Oberfläche 
positioniert. Beide Studien sind mit kombinierter CT-Bildgebung durchgeführt. Für 
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Detektor 1 




Abb. 3.25. Tiefenkarten der Mausober flache für Detektor 1: Die auf Grundlage der seg- 
mentierten Oberfläche berechneten AbbUdungstiefen können direkt für die fokussierte Rück- 
wärtsabbildung verwendet werden. Die Abbildungstiefen sind entsprechend der Skala am 
rechten Rand farbkodiert dargestellt. Unter der Voraussetzung, dass die Maus unverändert 
liegen bleibi, lassen sich beliebige Sigualinessungen du ich füll reu und mit diesen Tiefenkarten 
auswerten. 

das dritte Zylindcrphantom wird als Signalquelle ein angeregter Fhuiieszenzfarhstoff 
verwendet. Für diese Messung ist der Aufbau um einen optischen Filiersatz und einen 
Anregungslaser erweitert worden. 

3.2.3.1 Lumineszenzquellen auf der Oberfläche 

In dieser Studie sind auf der Oberfläche eines Zyliuderpliantoins (0 24 mm) au zwei 
Stellen Luinineszenzquellen (s. Abschn. 3.1.4) angebracht. Um die Signalquantifizierung 
aller vier Detektoren zu vergleichen, sind au einer Position zwei Quellen dicht neben- 
einander gesetzt. Dabei gilt dann die Annahme, dass dort auch die doppelte Intensität 
zu messen ist. 

Die extrahierte Oberfläche (153 Proj., je 0.5 s Belichtungszeit) und die Vorstufen sind 
in Abb. 3.26 dargestellt. Die kleinen Glaszylinder sind mit durchsichtigem Klebeband 
fixiert, das ebenfalls in den Projektionen und Sinogrammen an dessen Lichtreflexion zu 
sehen ist. Die gerenderte Isooberfläche ist dem Zylinderphantom (hier ohne Quellen) 
gegenübergestellt. Mit dieser Oberfläche werden dann die Abbildungstiefen für die 
Signalim 'ssung here< linet . 

In Abb. 3.27 sind fusionierte Schichten aus der kombinierten optischen und CT- 
Bildgebung gezeigt. Jeder CT-Schicht (Graustufen) ist die entsprechende optische 
Schicht als Binärbild (weiß) überlagert. Die Transversalebene zeigt einen Querschnitt 
durch die Position, an der sich die zwei dicht nebeneinander liegenden Signalquellen 
befinden, während in der Sagittalebene der Abstand der beiden Positionen zueinander 
zu sehen ist, Die maximale Abweichung der optisch extrahierten Oberflächen gegenüber 
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Isooberfläche Phantom 



Abb. 3.26. Optische Kxtraktion der Zylindernberrlächc. nuT > U i- ;m zwei Stellen Lumincs- 
zeir/quellen mit Klebestreifen befestigt sind. Die Positionen der Quellen selbst sind in den 
Projektionen. Sinogrammen und auf der Isooberfläche schon hier erkennbar. Das Zylinder- 
phautom ist. ohne aufgeklebte Quellen abgebildet. 




TE KE 



Abb. 3.27. Zyliuderoberfläche in kombinierter optiseher und C'T-Bildgebung: Die CT- 
Schiehten (Graustufen) sind einmal separat und einmal mit der optisch extrahierten Iso- 
oberfläche (weiß) überlagert dargestellt. Die auf der Oberfläche angebrachten Lumineszenz- 
queüen sind im CT-Datensatz gut erkennbar. 

der sichtbaren Zylinderoberfläche in der CT-Rekonstrnktion beträgt auf der Unters* it. 
des Phantoms um die Haltestäbe etwa 1,9mm. Für den übrigen Oberflächenbereidi 
liegen die Abfltaildsdifl'erenzen unter 1,4 mm. 

Die Signalprojektionen von Detektor 1 (153 Proj.. je 3s Belichtungszeit) sind in Abb. 3.28 
dargestellt. Zusätzlich sind die Maximalwerte beider Quellenpositionen winkelabhängig 
im Diagramm dargestellt. Das fehlende Stück der Messwerte resultiert aus den fehlenden 
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0[TT/8Orad]: 1 
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Abb. 3.28. Signalprojektionen von Detektor I: An der oberen Position befindet sieb eine 
Luniiiieszenzquelle, während darunter zwei Quellen dicht nebeneinander sitzen. Die Maxi- 
iiialiuteusitäteii sind im Diagramm \v i nk^ -l i ) > l i;i 1 1 dargestellt. Die Lin k'' resultier! aus den 

konstruktionsbedingt fehlenden Winkelstelhingen des Datensatzes. Das erwartete 2: 1- Verhältnis 
ist in den Signalprojektionen gut zu erkennen. 



Winkelpro jektiouen des Scanners. Schon in dieser Darstellung ist das 2:1- Verhältnis der 
beiden Quellenpositionen zu erkennen. 



CT optisch / CT 




Abb. 3.29. Optische Signalmcssung mit kombinierter (T-Bildgebung: Die CT-Schichten 
(Graustufen) des CT-Datensatzes uns Abb. ii.'iT sind auf der linken Seite separat und rechts mit 
den aufsuinniierten Signalintensitäten von Detektor 4 (Falschfarben) überlagert dargestellt. 



Um ein stärkeres Signal als in einer einzelnen Signalprojektion zu erhalten, können alle 
153 Signalprojektionen eines Detektors direkt auf der Oberfläche aufsummiert werden. 
Abb. 3.29 stellt die Oberflächensumme von Detektor 1 an beiden Quellenpositionen 
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Detektor Detektor 



Abb. 3.30. Auswertung der aufsuininierten Oberflächensignale der einzelnen Detektoren: 
Die Mittelwerte basieren auf einer Statistik über die Voxelintensitäten, die mindestens 50% 
des Maximalwertes in einer kleinen Umgebung um die jeweilige Quellenposition aufweisen. 
Der linke Balken repräsentiert die Musswerte aus dein Bereich der Doppelquelle. Das nahezu 
2:1 -Verhältnis ist in jeder Messung aber mit Variation zwischen den Detektoren sichtbar. Das 
Messvolumen der beteiligten Voxel ist für beide Quellenpositionen ähnlich groß. 

(Falschfarben) im Vergleich zum zuvor aufgenommen CT-Datensatz (Graustufen) dar. 
Im optischen Datensatz ist die jeweilige Quellenposition über den Signalschwerpunkt 
bestimmt worden. Im CT-Datensatz wurde dagegen der Mittelpunkt der gut erkennbaren 
Glaszylinder ermittelt. Der euklidische Abstand dieser Positionen beträgt 0,4mm für die 
Einzelquelle und (1.22 mm für die Doppelquelle, sodass der Lokalisierungsfehler in beiden 
Fällen kleiner als die Quellenausdehnung ist. Für die optische Signalquantilizierung in 
Abb. 3.30 ist um jede Quelle eine 3D-Box gelegt und alle Voxel innerhalb dieser Box. 
die mindestens 50% des Maximalwertes erreichen, in die Statistik einbezogen worden. 
Die so berechneten Mittelwerte stehen wie erwartet für alle vier Detektoren nahezu im 
Verhältnis 2:1. variieren jedoch zwischen den Detektoren. Der Grund hierfür ist, dass 
die Detektoren 2 und 3 durch ihre Anordnung nur einen Teil des Phantoms erfassen 
und somit di< Quellen in weni,u;i t l'i o j. .ii> i n zu sehen sind. Das Messvolumen diese! 
Auswertung, d.h. der von den beteiligten Voxeln gefüllte Raum, unterscheidet sich nur 
unwesentlich zwischen den Quellenpositionen und liegt im Schnitt aller Detektoren bei 
etwa 2,5mm 3 . Das Gesamtvolumen einer Lumineszenzquelle liegt bei etwa 1,6mm 3 . 

3.2.3.2 Lumineszenzquellen unter der Oberfläche 

In dieser Studie sollen Oberflächensignale von vier Lumineszenzquellen (zwei grüne, zwei 
etwas schwächer leuchtende rote) gemessen werden, die nun wenige Millimeter unter 
der Oberfläche eines Zyliuderphantoms (0 = 28 mm) sitzen. Diese Versuchsanordnung 
entspricht den Voraussetzungen für die Biolumineszenztoinographie. 

Die extrahierte Oberfläche (153 Proj., je 0,5 s Belichtungszeit) des Zylinderphantoms 
ist in Abb. 3.31 gezeigt. Der Zylinder liegt hier nicht mit seiner kompletten Länge im 
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Scannerobjektraum. Der Teil des Objektraumes, der nicht vom Zylinder eingenommen 
ist, wurde vorder Segmentierung abgeschnitten, da hier in den normierten Projektionen 
Störsignale dominieren. Die auffälligen Kreisartefakte im Innenbereich des Zylinders 
resultieren aus Spiegelungen der Leuchtfolie an der Phantomoberfläche und sind nicht 
weiter relevant. 




Isooberfläche Phantom 



Abb. 3.31. Optische Oberflächenextraktion des Zylinderphantoms. Die Luinineszenzquellen 
sitzen unter der Oberfläche und sind deswegen hier nicht sichtbar. An der Phantomober- 
fläche Spiegelt Sich teilweise die I.enchtfolie. Dies führt zu den auffälligen Kreisartefakt en 
innerhalb der Oberfläche in den Olijektrauinschicliteu. die jedoch keinen EinfluSS auf die 
Oberfiächenseginentierung haben. 

Die überlagerten Datensätze der kombinierten optischen und OT-Bildgebung sind in 
Abb. 3.32 dargestellt. Die vier Quellen sitzen in einer Transversalebene (sichtbar als 
weifie Punkte und besser zu erkennen in Abb. 3.34) in einiger Entfernung zum Zylin- 
derende. In diesem Bereich liegen die maximalen Oberflächenabweiehuugon in derselben 
Größenordnung wie in der vorangegangenen Phantomstudie. Am Zylinderende dagegen 
finden sich deutlich stärkere lokale Abstandsunterschiede (> 3 mm) als bei den zuvor 
untersuchten Objekten. Die zusätzlichen Lochst rukturen. die im CT-DatensatZ im 
Phantominneren zu sehen sind, haben für diese Studie keine Bedeutung. 

Die Signalprojektionen von Detektor 1 (153 Proj.. je 3 s Belichtungszeit) sind in Abb. 3.33 
dargestellt. Die vier Signalprojektionen stehen immer senkrecht zu einer der vier Lu- 
mineszenzquellen. Die entsprechenden Signalintensitäten sind zum Vergleich als Pro- 
Bldiagramine abgebildet. Der Maximalwert einer grünen Lumineszenzquelle erreicht 
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Abb. 3.32. Zylinderoberfläche in kombinierter optischer und CT-Bildgebung: Die CT- 
Schichten (Graustufen) sind separat und mit der optisch extrahierten Isooberflächc überlagert 
dargestellt. Die scharfe Schnittkante ain Zylinderende kann nicht gut mit der optischen 
Methode segmentiert werden. In der Transversalschicht, in der die vier Quellen sitzen, ist 
die Übereinstimmung jedoch vergleichbar mit den vorherigen Studien. Die Quellen {weiße 
Punkte, besser sichtbar in Abb. 3.31) sitzen in oberflächennahen Löchern. Die übrigen im 
CT-Datensatz sichtbaren Hohlräume haben hier keine Relevanz. 

dabei lediglich die halbe Intensität im Vergleich zur selben Quelle an der Oberfläche 
in Abb. 3.28. Die Signalintensität variiert mit der Leuchtstärke der Quellen und deren 
Tiefenpositionen im Phantom. Durch das schwächere Messsignal treten deutliche Stör- 
signale am Zylinderrand auf. verursacht durch Streulicht während der Messung, das 
nicht komplett unterbunden werden konnte. 

Für die Transversalebone. in der die Quellen sitzen, sind in Abb. 3.34 die Oberflächenin- 
tensitaten über alle 153 Projekt ionen von Detektor 1 aufsummiert (Falschfarben) und mit 
der entsprechenden CT-Schicht (Grauwerte) überlagert. Anhand der Quellenpositionen 
(weiße Punkte) lassen sich die resultierenden Oberflachenintensitäten erkennen. Zur 
besseren Darstellung sind die Intensitäten zusätzlich in einem Kreisdiagramm gezeigt, in 
dem zu jedem Winkel der Oberflächenmaximalwert als Radialkomponeute eingetragen 
ist. Die Signalintensität korreliert zum einen mit der Tiefe im Phantom, zum anderen 
leuchten die roten Lumineszenzquellen deutliche schwächer als die grünen Quellen. 
Letzteres ist allerdings der entscheidende Faktor, mit dem sich die unterschiedlich 
starken Signalintensitäten erklären lassen. 

3.2.3.3 Fluoreszenz unter der Oberfläche 

Für die Oberflächensignalniessung mit einem angeregten Fluoreszenzfarbstoff ist der 
Scanneraufbau um einen Filtersatz und einen Diodenlaser erweitert worden. Solch eine 
Anordnung wird typischerweise für die Signalmessung in der Fluoreszenztomographie 
verwendet. In dieser Studie sind für die Signahnessnng lediglich die beiden Signalde- 
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Abb. 3.33. Signalprojektionen von Detektor 4 jeweils senkrecht zu einer Lumineszeuzquelle 
unter der Oberfläche. In den Profildiagraminen sind die Sigitalinteiisitaten der vier Quellen 
dargestellt und direkt vergleichbar mit der Messung aus Abb. 3.28. Die hellen Bereiche mit 
Zylindcrrand sind durch Streulicht, während der Messung verursacht und treten aufgrund lies 
schwachen Messsignals so deutlich hervor. 





0,3 Int [a u ] 

Abb. 3.34. Optische Signahnessung mit kombinierter CT-Bildgebung: Dem CT-Datensatz 
(Grauwerte) aus Abb. 3.32 sind die aufcu initiierten Signalintensitäten von Detektor I (Falsch- 
farben) überlagert. Die roten (r) und grünen (g) Lumineszenzcptellen sind als weiße Punkte 
markiert. Im Krei.sdiagratnm sind die Maximalwerte der Oberflächeninteiisitäten als Radial- 
komponente dargestellt, in dem die vier Quellensignale zu erkennen sind. 



Lektoren eingesetzt worden. Zwei Langpasslilter RG715 (Schott AG, Mainz, Deutschland) 
aus Farbglas mit je 3 mm Dicke und einer Durchlässigkeit ab 715 ±9 um sind als Fenstrt 
in einer lichtdichten Abdeckung eingesetzt, die nach der Oberflächenexüaktion vor die 
Detektoren geschoben wird. Dadurch sind die Detektoren 2 und 3 komplett abgedunkelt, 
während die Signahnessung der Detektoren 1 und 4 hinter den Filtern durchgeführt wird. 
Zusätzlich ist auf der Oberseite der Abdeckung ein Anregungslaser (Laser Components 
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GmbH. Olching, Deutschland) mit einer Wellenlänge von 035 nm und einer Leistung 
von 1 mW angebracht. Der Laser und die beiden Signaldetektoren sind durch den 
Einschub stets rechtwinklig zueinander angeordnet. Diese Abdeckung wird dann mit der 
Ruial ionsplat t form um das Objekt gedreht, um Signalini cnsilüleu aus verschiedenen 
Wiukelstellungen messen zu können. Die Scannererweiterung ist in Abb. 3.35 gezeigt. 




(a) (b) 



Abb. 3.35. Einschub für die Fluoreszenzmessimg als schein atisrhe Darstellung (a) und 
während des Experiments im Scanncraufbau (b). Die Signaldctcktorcn (Detektoren 1 und 
4) sind durch die Langpassfilter abgedeckt. Der Einschub wird nach der Messung für die 
Oberflächencxtraktion eingesetzt . 

Die optisch extrahierte Zylinderoberfläche (153 Proj.. je 0.5 s Belichtungszeit) ist in 
Abb. 3.36 gezeigt. Die Messung für die Oberflächencxtraktion ohne den Einschub erfolgte 
vor der Signalmessung. Der Schwellwert für die Oberflüchensegnientieruug wurde so 
gewählt, dass die resultierende Oberfläche den bekannten Zylinderdurchmesser (24 min) 
aufweist. In der Phantommitte befindet sich eine Bohrung {0 - 2 mm), die über die 
gesamte Phantomlange mit Cy5.5 (Ameraham Bioscience, Piscataway, N.J. USA) in 
einer Konzentration von 8' Mno, /i befüllt war. Über das Laserlicht wird der Farbstoff 
mit einer bestimmten Wellenlänge angeregt und emittiert als Folge Licht mit einer im 
Spektrum verschobenen Wellenlänge. Maximalintensitäten werden mit Cy5.5 für die 
Anregung bei 673 um und für die Emission bei 692 nm erreicht [53]. Somit ist die in 
diesem Aufbau verwendete Laset -Filter-Kombination hinsichtlich der Signalausbeute 
nicht optimal gewählt, aber für diese Tcstincssung ausreichend. 
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(*) (b) 

Abb. 3.36. Optisch extrahiert«' Isooherfiache (a) des Zylinderphantoins (h). Das Phantom 
hat. in der Mitte ein Bohrloch, in das der FluoreszenzfarhstoH' eingefüllt ist. 

Die Signalprojektionen von Detektor 4 { 153 Pro j.. je 3 s Belichtungszeit) sind in Abb. 3.37 
dargestellt. Das Phantom ist zwar über die gesainte Länge befällt, der Farbstoff wird 
jedoch vom Laser nur in einem bestimmten Bereich angeregt. Aus diesem Grund er- 
Streckt sich das Signal nicht über die gesamte Länge dt >s Phantoms. Bei der Rotation um 
das Phantom beeinflussen die Haltestäbe des Objektbettes die Signalmessung zweifach: 
einmal können sie das Anregungslicht blockieren und zum anderen das Emissionslicht 
venleekon. Diese Effekte sind deutlich in den Sinogrammen zu erkennen: Die waag- 
rechten Auslöschungen sind durch das blockierte Anrogungshiserlicht verursacht, das 
immer in derselben Winkelstellung auf die Haltestäbe trifft. Die schräg verlaufenden 
Abschwächungen ergeben sich ebenfalls aus den Haltestäben, die an diesen Positionen 
das Emissionslicht verdecken. Ohne Einfluss der Haltestäbe ist das Fluoreszenzlicht 
deutlich in den beiden linken Projektionen messbar. Blockieren die Haltestäbe da.s 
Anregungslicht, so verschwinden die Signale wie in der mittleren Projektion fast voll- 
ständig. In den beiden rechten Projektionen sind die Oberflächenmessung stellenweise 
durch die Haltestäbe beeinträchtigt. Die delektierten Maximalintensitäten, die aus 
den vertikalen Profildiagraininen zu entnehmen sind, liegen etwas über den Werten 
der Lu m ineszenzq uellen in Abb. 3.28. Für eine tomographische Rekonstruktion dieser 
Fluoreszenzverteilung wäre zusätzlich noch die Intensität sinessung des Anregungslichtes, 
das aus der Phantomoberfläche heraustritt, erforderlich. Diese Messung ist nicht durch- 
geführt worden, da das Anregungslicht vielfach stärker ist als das Emissionslicht und 
den Sensor selbst bei kürzester Belichtungszeit (367 ms) noch sättigt. Eine mögliche 
Lösung wäre, das Laserlicht soweit abzuschwächen bis das Anregungslichl schließlich 
gemessen werden kann. 
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Abb. 3.37. Bmiasionslichl des angeraten Fluoreszenz farbstoffs, gemessen mit Detektor 1. 
Der Farbstoff wird von der punktförmigen Laserquelle nur in einem bestimmten Bereicli 
;iu;j,i'jv -i . Di.' Sijjualabsciiwäciunigeii in den Sinogranin* n und Projektionen sind veriirsaelu 
durch <lie Haltestäbe des Objektbettes, die einmal das Anregungslicht und ein andermal das 
Kmissionsliclit blockieren. 
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Diskussion 



Der Schwerpunkt dieser Arbeit liegt auf der Bildgehung mit Mikrolinsen-Detektorcn im 
Hinblick auf die optische Signahnessung. wie sie für die tomographische Rekonstruktion 
im Kleintierbereich erforderlich ist. Als optisches Element des Detektors ermöglicht 
das MLA einerseits einen sehr flachen Detektoraufbau, andererseits jedoch setzt die 
M LA- Abbildung eine spezielle Bearbeitung der Seusordaten voraus, um Objektbilder 
vergleichbar mit denen einer klassischen Kamera zu erhalten. Bisher sind aus den 
Seusordaten Fokuspunktbilder extrahiert worden, für die lediglii Ii I /< der Seusorpixel 
verwendet wurden. Die verworfenen Pixelintensitäten eines Fokuspunktbildes enthalten 
allerdings wertvolle Informationen über den Objektraun i. Das in dieser Arbeit entwi- 
ckelte .V bil lungsmodell hingegen l>hn!>i /us.'ii/'i-.-h.- S.-n- rpixel (etwa 64 'X ). die nicht 
von der Septenmaske abgedeckt werden, gleichzeii in die Bildbereclniung mit ein. 
Dadurch ist die Bildqualität inessbar gesteigert, und die Tiefeninformation zeigt sich in 
der tiefenabhängigen Intensitätsverteilung der Objektbilder. 

Das Detektor-Abbildungsmodell ist allgemein formuliert und kann auf jede Detektorgeo- 
metrie übertragen werden, solange die paraxiale Näherung für die Mikrolinsen gültig ist. 
Uber die Rückwärtsahhilduug lassen sirli aus den Seusordaten normierte Objektbilder 
in beliebigen Abbildungstiefen berechnen. Die Vorwärtsabbildung kann für umfang- 
reiche Modellstudien eingesetzt werden, mit denen sich die Abbildungseigenschaften 
theoretisch untersuchen lassen. Pixel- Voxel- Beziehungen werden über den vereinfachten 
Strahlenverlauf durch die optischen Zentren der Mikrolinsen bestimmt. Die Einflüsse 
von Reflexion. Diffraktion und der Wellenlängen des Lichtes sind vernachlässigt. Die 
Dicke des MLAs wird näherungsweise über die effektive Fokuslänge berücksichtigt. Die 
Modellstudien und experimentellen Untersuchungen in dieser Arbeit zeigen, dass das 
Abbildungsmodell auf den optischen Detektor angewendet werden kann. 
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Der Detektoraufbau ist vergleichbar mit dem der Compound-Eye-Kamera [15], die 
jedoch eine zehnmal kleinere Fläche als der in dieser Arbeit eingesetzte optische De- 
tektor aufweist. Zudem ist in der Compound-Eye-Kamera jede Mikrolinse lediglich 
mit einem Sensorpixel assoziiert. Dadurch lassen sich direkt zusammengesetzte Bilder 
aufnehmen, die somit den Fokuspunkt bildern des optischen Detektors entsprechen. 
Konzepte, die wie das Abbildungsmodell mehrere Sensorpixel mit einer Mikrolinse in 
Verbindung stellen und darüber tiefenabhängige Objektbilder berechnen, basieren auf 
Experimentalanordnungen mit größeren Abmessungen als die des optischen Detektors. 
Ng et al. integrieren das MLA in eine tragbare Handkamera samt Objektiv, sodass 
sich die Aufnahmen nachträglich digital refokussieren lassen [17]. Hong et al. hingegen 
verwenden eine separate CCD-Kamera hinter dem MLA. um die Elementarbilder (engl. 
elemental images) eines nahen Objektes zu erfassen [21]. Lee et al. erweitern diesen 
Aufbau um eine Einzellinse, damit auch weiter entfernte Objekte abgebildet werden 
können [54]. Im Vergleich zu diesen Aufnahinesyst einen, hat der optische Detektor ei- 
nerseits kompakte Abmessungen, kann dabei aber gleichzeitig die Vorteile der variablen 
Abbildungstiefe nutzen. In dieser Arbeit ist gezeigt, wie sich dieses Detektorprinzip 
effizient für die Kleintierbildgebung einsetzen lässt. 

Die Abbildungseigenschaften des optischen Detektors sind in einem Tiefenbereich von 
10 mm - GOmm untersucht, der durch den Scanneraufbau vorgegeben ist. Modell- und 
Experimentalstudien haben gezeigt, dass die räumliche Auflösung in den Objektbildern 
in erster Line von der festgelegten MLA-Fokuslänge und der variablen Abbildungstiefe 
abhängt, Dabei verschlechtert sich di< r.'iuiihUhe Aullösung mit ansteigender Abbil- 
dungstiefe. Im Vergleich zu einem Fokuspunkt bild liefert das Abbildungsmodell auf 
Grundlage derselben Sensordaten dennoch Objektbilder mit besserer räumlicher Auf- 
lösung über den gesamten Tiefenbereich des Scannerobjektraumes. Besonders deutlich 
ist diese] Unterschied unmittolbai vor dem Detektor 1/ • 20 mm). AI- Konsequenz des 
abnehmenden Auflösungsvermögens werden die Detektorprojektionen im Scanneraufbau 
mit einer achtfach verringerten Auflösung gegenüber der Sensoroberfläche berechnet, da 
ab einer Abbildungstiefe t > 30 mm kaum noch höhere Raumfrequeuzeu rekonstruierbar 
sind. Auf diese Weise lässt sich die Datenanalyse der Detektorprojektionen erheblich be- 
schleunigen. Die Signalprojektionen könnten durchaus mit einer etwas höheren Aullösung 
berechnet werden, da die typischen Abbildungstiefen hier zwischen 15 mm - 30 mm 
liegen. Das Signal-Rausch- Verhältnis ist auf Grundlage derselben Sensordaten je einmal 
im Objekt- und Fokuspunktbild mit einer schwachen Lichtquelle bestimmt, wobei der 
Quellenabstand und die Belichtungszeit variiert wurden. Als Ergebnis findet sich in 
den Objektbilden] grundsätzlich das bessere SNR. das in kleinen Abbildungstiefen und 
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mit längeren Belichtungszeiten höchste Werte annimmt. Die iterative Anwendung des 
Abbildungsmodolls kann in kleinen Abbildungstiefen den Bildkontrast in den Objekt- 
bildern zusätzlich verbessern, bietet aber keine Verteile im relevanten Tiefenbereich des 
Scannerobjektraums, der den gesteigerten Reehenaufwand rechtfertigen würde. 

Der optische Detektor soll für die kontaktlose Oberflächensignalmessung anstelle einer 
CCD-Kamera eingesetzt werden [55], um mit minimalem Platzaufwand zukünftig 
optische Tomographie für eine Maus zu ermöglichen. Über das Abbildungsmodell 
kann die Signalmessung direkt auf der Objektoberfläche durchgeführt werden. Sowohl 
für die exakte Rückwärtsabbildung der Signalmessung als auch für die zukünftige 
tomographische Rekonstruktion ist die Kenntnis der Objektoberfläche notwendige Vor- 
aussetzung [12, 50, 56]. In dieser Arbeit ist die Tiefenbestimmung anfangs aus einer 
einzelnen Detektoraufnahme und schließlich über mehrere Detektorprojektionen im 
Scanneraufbau untersucht worden. 

Bereits aus einer einzelnen Sensoraufnahme lässt sich ein Objektdatensatz berechnen, in 
dem die Tiefeninformation über das Objekt enthalten ist. Basierend auf der Analyse der 
Vi ■xeliiitonsitäten bzw. Varianzen sind experimentell Tiefenpunkte bzw . s< h >n ganze 
Tiefenkarten berechnet worden. Limitierend wirkt hier in jedem Fall die tiefenabhängige 
räumliche Auflösung in den Objektbildern, sodass die Tiefenbestimmung bei steigender 
Objektentfernung immer schwieriger wird. Ebenso ist hier eine passenden Objektbe- 
leuchtung notwendig, um die erforderlichen Strukturen für die Tiefenanalyse sichtbar 
zu machen. Anders als bei digitalen Hologrammen, von denen die Varianzmethode 
adaptiert ist. zeigt der Objektdatensatz zudem eine zu geringe Tiefenauflösung, mit 
der die resultierenden Tiefenkarten nicht an die Ergebnisse von McElhinney et al [32] 
herankommen. Grundsätzlich sind die Resultate aus einer Sensoraufnahme deutlich 
schlechter als die Ergebnisse der Oberflächenextraktion mit mehreren Detektorprojek- 
tionen. 

Die Stärken des Detektors liegen nicht so sehr in der Einzelanwendung, sondern aufgrund 
der kompakten Abmessungen in der effizienten Kombination mehrerer Detektoren. Für 
die Ringanordnung der vier Detektoren mit etwa 8 cm Durchmesser auf der Rotati- 
onsplattform ist ein Kompromiss erarbeitet, mit dem clor Scanneraufbau für die zwei 
komplementären Aufgaben der Oberflächenextraktion und Signalmessung autark ein- 
gesetzt werden kann. Jeder der Detektoren kann als eigenständige Kamera betrachtet 
werden, sodass bereits ohne Rotation das Objekt über vier Winkelprojektionen gleich- 
zeitig erfasst wird. Durch die Scannerrotation kann das Objekt in vollem Umfang 
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abgedeckt werden. Wichtig für die Koordinierung der Detektoren ist die exakte geome- 
trische Seannerkalibrierung. die mit Hilfe eines dünnen Stäbchens demonstriert ist. Über 
die bekannten Lageparameter der Detektoren lassen sich die DKSe der Oberflächen- 
und Signaldetektoren einfach ineinander überführen. 

Für die Oberflächenextraktion mit dem Scanneraufbau ist gezeigt, wie aus Reflexions- 
projektionen unter bestimmten Bedingungen die Objektoberfläche segmentiert und die 
Abbildungstiefen für die fokussierte Rückwärtsabbildung berechnet werden können. Um 
die Voraussetzungen für eine Nacktmaus allesamt zu erfüllen, sind zwei Detektoren 
im Scanner speziell angeordnet und deren Projektionen für eine FOV-Vergrößerung 
zusammengefügt. Dabei können einige konstruktionsbedingt fehlende Detektorprojek- 
tionen sogar effizient interpoliert werden. 153 Winkelpositionen und eine Messdauer 
von etwa 2 min liefern ein adäquates Ergebnis, das gute Übereinstimmung mit einem 
simultan aufgenommenen (T-Datensatz zeigt. Schartkantige Bereiche und die Halte- 
rungen des Objektbettes führen zu etwas größeren lokalen Abweichungen der Oberfläche. 
Die Segmentierung erfolgt über einen manuell abgestimmten Schwellwert für jeden 
Oberflächendatensatz. Hier lässt sich das Verfahren mit effektiveren Segmentierungs- 
methoden [25] durchaus noch verbessern und automatisieren. Zudem ist die optische 
Oberflächenextraktion bisher nur an unbewegten Objekten getestet und muss sich 
im nächsten Schritt erst noch als robust gegenüber Bewegungen (z. B. Atmung und 
Herzschlag) erweisen. 

Über die nahezu 360 °- Abdeckung des Objektes kann eine geschlossene Oberfläche 
im vollen Umfang des Objektes gefunden werden, anders als bei nicht rotierbaren 
Systemen [50, 56, 57]. Die Oberfläche wird direkt im Koordinatensystem des kali- 
brierten Scanners segmentiert. Somit ist eine Transformation in die einzelnen DKSe 
ohne zusätzliche Bezugspunkte möglich [56] und kann allein über die bekannten Lagepa- 
ranieter erfolgen. Perspektivische Verzerrungen, wie sie bei der Oberflächenextraktion 
mit einer rotierbaren CCD-Kamera und Weitwinkelobjektiv auftreten [12], sind mit 
den Oberflächendetektoren nicht zu beobachten. Für die Detektoren kann tatsächlich 
eine Parallelstrahlgeometrie für die gefilterte Rückprojektion angenommen werden, 
ohne dass große Abweichungen in der segmentierten Oberfläche entstehen. Die Oberflä- 
chensegmentierung ist in dieser Arbeit mit kombinierter CT-Bildgebung unabhängig 
validiert worden und zeigt lediglich im Bereich der Haltestäbe und scharfkantiger Ecken 
größere lokale Abweichungen. Ansonsten lässt sich eine Genauigkeit wie in vergleichbaren 
Studien erreichen [12. 50, 56]. Ein Nachteil der hier verwendeten Detektoren ist, dass 
eine Nacktmaus nicht in ihrer ganzen Länge erfasst werden kann. Eine direkte Verbes- 



Copyrighted material 



97 



seruug des Sennneraufbaus wäre möglich, wenn die aktive Detektorfläche sowohl in der 
Breite als auch in der Länge (z. B. 4 x 10 cm 2 ) eine Nacktmaus komplett abdeckte. Die 
Einzeldetektoren ließen sieh dann alle zentriert anordnen und vollwertig sowohl für die 
Oberflächensegmentierung als auch für die Signalmessung einsetzen. Als Folge könnten 
entweder bei einer vollen Umdrehung redundante Informationen zur Ergebnisverbes- 
serung genutzt oder durch die lediglich erforderliche Teilrotation die Messdauer verkürzt 
werden. 

Der Scanneraufbau ist deutlich kompakter gegenüber Systemen mit CCD-Kameras [12, 
50, 56, 57]. Zudem sind gleichzeitig vier Detektoren für die Bildgebung einsatzbereit. 
Es muss jedoch noch gezeigt werden, dass mit den gemessenen Signalintensitäten eine 
tomographische Rekonstruktion im gleichen Maße möglich ist. Die Signahnessung ist 
unter den Voraussetzungen, wie sie bei der BLT oder FMT zu erwarten sind, getestet 
worden. Beide Arten sind im Continous-Wave-Modus [9] durchgeführt, bei dem die 
Lichtintensität einer konstanten Lichtquelle gemessen wird. Ein kompletter Messablauf 
sieht zwei Schritte vor: zuerst die optische Oberflächenbestimmung und anschließend die 
Signalmessung. Die Datenaufnahme läuft automatisiert, wohingegen die Datenanalyse 
stellenweise einen manuellen Eingriff fordert (z.B. Rauschfiltor. Normierung. Schwellwert 
für die Segmentierung) und somit noch Verbossorungspotential bietet. 

Mit schwach leuchtenden Lumineszenzquellen sind Signalmessungen durchgeführt 
worden, die den Voraussetzungen für die BLT entsprechen. Au der Oberfläche sitzend, 
sind diese Quellen deutlich in den Signalprojektionen eines Detektors lokalisierbar und 
können relativ zueinander quantifiziert werden. Sitzen die Quellen wenige Millimeter 
unter der Oberfläche, dann sind die Signale unter gleichen Bedingungen noch eindeutig 
messbar. Der maximale Signalwert lä.sst sich etwa um das Zehnfache steigern, wenn die 
Intensit äten aller 15:3 Signalpro jekti .neu eines 1 letektors an di r < >berflä< he aufsummiert 
werden. Im direkten Vergleich erweist sich erwartungsgemäß die zent rierte Anordnung 
der Signaldetektoren als effektiver für die Signahnessung gegenüber der unzentrierten 
Anordnung der Oberflächendetektoren. Der Einsatz mehrerer Detektoren gleichzeitig 
bietet sich zudem für eine spektral aufgelöste Messung [50. 58] an, indem jeder Detektor 
mit einem unterschiedlichen Bandpassfilter ausgerüstet wird. Bioluniineszenz ist tat- 
sächlich deutlich schwächer zu detektieren als Fluoreszenz [9]. So bleibt zu untersuchen, 
ob die Detektorseusitivität grundsätzlich hierfür ausreichend ist. 

Ein Messdurchgang, wie er für die FMT erforderlich ist. wurde anhand einer heraus- 
nehmbaren Detektorabdeckung mit Filtereinsatz und Anregungslaser demonstriert. Die 
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Maxiiualintonsiläten in den Signalprojektionen liegen in der Größenordnung der zuvor 
untersuchten Luinineszenzquellon, obwohl «1er Fluoreszenzfarbstoff sehr viel tiefer im 
Phantom eingebracht ist. Aufgrund der suboplhnalen Lascr-Filter-Konfiguration war 
dennoch eine vergleichsweise hohe Cy5.. ^-Konzentration [12. 56. 59] notwendig, um ein 
messbares Signal zu erzeugen. Der Anregung.slaser lässt sich zwar mit dem Einsehub 
rotieren, ist darin aber fixiert und sollte zukünftig für eine bessere Quellen-Detektor- 
Kombination ebenfalls in z-Richtung beweglich sein [12]. Weiter muss noch eine Lösung 
gefunden werden, um eine Signalmessmig des Aiiregmigslicht.es durchzuführen, da das 
Laserlieht sehr viel stärker ist als die Fluoreszenz und den Sensor selbst bei minimaler 
Belichtungszeil bisher sättigt. Auch hier bietet sich dann aber die Möglichkeit, das 
Anregungs- und Emissionslicht gleichzeitig über mehrere Detektoren zu erfassen. 

Als Resultat aus den Signalmessungen ist die Lichtempfindlichkeit, des eingesetzten 
Sensors grundsätzlich kritisch zu bewerten. Zusätzlich einschränkend wirken sich die 
minimale Auslese- und die maximale Belichtungszeit aus. die im Vergleich zu den 
hochsensitiven CCD-Kameras besonders lang bzw. kurz sind. Längere Messzeiten ließen 
sich über mehrere hintereinander aufgenommene Sensordatensätze durchaus realisieren, 
wobei der 1 lewinu beim SM! no< Ii genau zu untersuchen wäre. Kürzere Belichtungszeiten 
sind durch die minimale Auslesezeit (367 ms) limitiert. Abgesehen davon lassen sich 
aber mit dem Sensor ohne Weiteres die Hachen Detektoren und der Scanneraufbau um- 
setzen, sowie die für die MLA-Bildgebung entwickelten Methoden eingehend untersuchen. 
Darüber hinaus kann der Sensor problemlos für die Überflächenextraktion eingesetzt 
werden, da hier die Lichtintensität über die Beleuchtung gesteuert ist. Letztlich wird 
die Anwendung tomographischer Verfahren aber davon abhängen, mit welcher Empfind- 
lichkeit sich die notwendigen Signalintensitäten detektieren und quantifizieren hissen. 
Xichtsdestotrotz können die in dieser Arbeit entwickelten Methoden auch auf Detektoren 
mit lichtsensitiveren Sensoren gleichermaßen angewendet werden. 

Der vielversprechende Vorteil, den kompakten Scanneraufbau in eine zusätzliche Mo- 
dalität zu integrieren, wurde in mehreren Experimenten mit kombinierter optischer 
und CT-Bildgebung demonstriert. Für die Registrierung der beiden Datensätze waren 
lediglich Rotations- und Translationsverschiebungen notwendig. So konnte einerseits die 
optische Oberflächenextraktion validiert und andererseits die strukturelle Darstellung 
der CT-Bildgebung mit der optischen Oberflächensignalmessung ergänzt werden. Extrem 
nachteilig hat sich jedoch die applizierte Röntgenstrahlung auf die Sensoren ausgewirkt, 
bei denen nach wenigen Messdurchläufen bereits ein irreversibel erhöhter Rauschoffset 
festzustellen war, der wiederum zur Verringerung der optischen Messempfindlichkeit 
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führt. Die Vorteile l>eider Modalitäten lassen sich nur dann sinnvoll nutzen, wenn die 
Nachteile auf ein Minimum reduziert werden und ein längerfristiger Betrieb gewährleistet 
ist. Gerade der platzeffiziente Scanneraufbau auf Grundlage der optischen Detektoren 
bietet großes Potential auch für die Kombination mit der PET [6] und der Magnetreso- 
nanztomographie (MRT) [(>()]. 

Abschließend lässt sich nun mit Blick auf die eingangs definierte Zielstellung (Abschn. 1.1) 
folgendes Resümee ziehen: In dieser Arbeit ist gezeigt, dass mit den optischen De- 
tektoren Oberflächensignale gemessen werden können. Basierend auf dem Prinzip eines 
Einzeldetektors und dessen charakterisierten Abbildungseigenschaften ist ein optischer 
Kleintier Scanner aufgebaut und optimiert worden. Einerseits kann mit diesem Gerät 
die optische Oberflächensignalmessung unabhängig durchgeführt werden, andererseits 
erlaubt der platzefhziente Aufbau tatsächlich eine kombinierte Bildgebung mit zusätz- 
lichen Modalitäten. In experimentellen Phantomstudien ist die kontaktlose Oberflächen- 
signalmessung für Lumineszenz- und Fluoreszenzsignale demonstriert worden. Diese 
Signalmessung und die extrahierte Oberfläche bieten nun die Grundlage für eine an- 
schließende tomographische Rekonstruktion. Anstelle einer herkömmlichen Kamera ist 
in dieser Arbeit die optische Bildgebung erfolgreich mit den eigens hierfür entwickelten 
Mikrolinsen-Detekroren realisiert worden. 
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Zusammenfassung 



In dieser Arbeit wird die Anwendung sehr Macher Mikrolinsen-Dotektoron für die optische 
Signahnessung untersucht, mit denen sich ein kompakter Scanneraufhau für Labormäuse 
realisieren lässt. Der Scanner soll zukünftig für optische Tomographie eingesetzt werden. 
Für die Bildgebimg mit den Detektoren ist ein Abbildungsmodell entwickelt worden, 
das deutlich bessere Bilder liefert als die bisherige Fokuspunktbildmethode. Im Unter- 
schied zu einer klassischen Kamera können Objektbilder mit dem Abbildungsmodell 
in beliebigen Abbildungstiefen nach der Aufnahme erzeugt worden. Für unbekannte 
Abbildungstiefen lässt sich so der Objektraum bereits auf Grundlage einer einzelnen 
Detektoranfnahme durch fokussieren und auf die enthaltene Tiefeninformation ana- 
lysieren. Mit bekannten Abbildungstiefen können direkt ein fokussiertes Objektbild 
oder im Falle einer dreidimensionale Oberfläche die dort lokalisierten Voxelintensitäten 
berechnet werden. Die räumliche Auflösung in den Objektbildern nimmt aufgrund der 
MLA-Ceometrie mit ansteigender Abbildungstiefe ab. Die Abbildungseigenschaften der 
Detektoren sind für den relevanten Tiefenbereich (15 mm — 50 mm) des Scanneraufbaus 
in Modellstudien und experimentellen Untersuchungen eingehend charakterisiert. Die 
Objektbilder des Abbildungsmodells zeigen für diesen Bereich eine höhere räumlich*» 
Auflösung und ein besseres Signal- Rausch-Verhältnis für schwache Lichtquellen als die 
entsprechenden Fokuspunktbilder. Für jeden Detektor sind passgenaue Rahmenteile 
entworfen worden, die die Detektorkomponenten (MLA, Septenmaske und CMOS- 
Sensor) justiert und robust gegen Erschütterungen zusammenhalten. Der rotierbare 
Scannerautbau besteht aus vier Einzeldetektoren in einer Ringanordnung mit 8 cm 
Durchmesser. Mit einem dünnen Stäbchen lassen sich die Lageparameter der Detektoren 
aus deren Sinogrammen bestimmen, sodass alle vier Detektoren gleichzeitig eingesetzt 
werden können. Es wird gezeigt, wie sich für eine Nacktmaus die beleuchtete Oberfläche 
über eine gefilterte Rückprojektion unter bestimmten Voraussetzungen ermitteln lässt. 
Für das Verfahren sind zwei Oberflächendetektoren im Scanner speziell angeordnet, 
die zusammen ein vergrößertes FOV liefern, sodass eine Mausoberfläche erfolgreich 
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segmentiert werden kann. Die Kenntnis der Objektoherfläche i-i einmal für die ('.«- 
kassierte? Rückwärtsabbildung der Sensordaten and zum anderen fär eine zukünftige 
tomographische Rekonstruktion der Signahnessung notwendig. Mit der optischen Ober- 
flächenoxtraktion and den daraas determinierten Abbildungstiefen ist anschließend die 
Oberflächensignalmossang in mehreren Phantomstudien demonstriert. Zum Einsatz 
kommen selbstleuchtenden Quellen und das Emissionslicht eines angeregten Fluores- 
zenzfarbstoffes. Die Segmentierung der Mausoberfläche sowie die Signalmessung mit 
den l.euclitquellen wurden mit kombinierter CT-Hildgobung durchgeführt, um die Er- 
gebnisse zu validieren und eine mögliche multimodale Bildgebung mit dem optischen 
Scaiineraufbau zu zeigen. 
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Ausblick 



Der Scann« Taufbau biotut die Möglichkeit bei minimalem Platzbedarf Oberfläehon- 
signale auf einer Nacktmaus zu messen. Mit der bekannten Objektuberftäclie und den 
dort gemessenen Signalintensitäten ist der nächste Schritt nun die tomographische 
Rekonstruktion des Objektinnern. Die Signalm essung ist mit dem Scanneraufl>au für 
Lumineszenz- sowie Fluoreszenzsignale demonstriert worden. Es muss sich hier noch 
zeigen, ob die Sensorempfindlichkeit der Detektoren für eine adäquate Rekonstruktion 
ausreichend ist. Grundsätzlich ist der Messablauf mit der zukünftigen Rekonstrukti- 
onsmethode abzustimmen, wobei das in dieser Arbeit vorgeschlagene Protokoll und 
der Scanneraufbau als Grundlage dienen können. Weiter muss untersucht werden, ob 
die vorgestellten Methoden ebenso für die m-two-Bildgebung einer narkotisierten La- 
bormaus anwendbar sind. 

Die kombinierte optische und CT-Bilclgebung zeigt bereits die Vorteile und Probleme, 
die sich aus der Integration des Scanneraufbaus mit einem zusätzlichen bildgebenden 
Verfahren ergeben. So können einerseits die Signalintensitäten von Leuchtquellen relativ 
zueinander quantifiziert und gleichzeitig mit der Strukturinformation aus der CT- 
Rekonstruktion überlagert werden. Andererseits muss aber sichergestellt sein, dass sich 
die beiden Modalitäten auch während eines dauerhaften Einsatzes nicht beeinträchtigen, 
wie es in dieser Arbeit der Fall war. Dabei wird die Integration mit den unterschiedlichen 
Modalitäten wie MRT. PET. CT und SPECT jeweils für sich spezielle Anforderungen 
an die optischen Detektoren und den Scanneraufbau stellen. Diese Herausforderungen 
gilt es zukünftig zu meistern, um von den Vorteilen der multimodalen Bildgebung zu 
profitieren. 
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Tab. A.l. Spezifikation des MLAs nach Herstellerangaben und eigenen Berechnungen (a). 
Der Krümmungsradius und die Position der Hauptebene ist per Definition positiv, falls das 
Kreiszentrum bzw. die Ebene rechts des Scheitelpunktes liegt [30]. Die Brennweite vom Scheitei- 
bis zum Fokuspuukt auf der konvexen Linsenseite entspricht bei einer plankonvexen Linse der 
Äquivalentbrennweite, da hier Scheitelpunkt und Hauptebene zusammenfallen. 





Toi 


MLA Dicke 


[mm] 


1 


±0.03 


Krümmungsradius 1 


[mm] 


00 




Krümmungsradius 2 


[mm] 


-1,81 


±0,05 


Äquivalentbrennweite 
(bei 633 nm) 


[mm] 


2.15 


±0,06 


Pos Hauptebene 1 |a> 


[mm] 


0,54 


±0,04 


Pos. Hauptebene 2 |a > 


[mm] 


0 




Brechungsindex 


1,82925 




Glassubstrat 


S-TIH53 
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Tab. A.2. Spezifikation des CMOS-Sensors mit dor Ausleseelektronik aus [27, 29J. Der 
mittlere Dunkelstrom verdoppelt sich etwa alle 8°C. Die Ausleseperiode entspricht der 
erforderlichen Übertragungszeit eines Bildes vom Sensor zum Speicher. Die minimale Bildrate 
ist die Empfehlung des Herstellers für die verwendete Ausleseelektronik und erlaubt eine 
Sensorbeliditung von maximal (17 s. 



Mittlere Dunkelstrom 

(bei 23° C) '* 1 


4000 


Leserauschen 

(rms, 1 fps) 16 1 


150 


Sättigung [e] 


2.810« 


Conversion gain [e/ADU] 


SOG 


AD-Wandlung [bit] 


12 


Quanteneffizienz 

(500-700 nm) 1 1 


>30 


Ausleseperiode [ms] 


367 


Max. Bildrate [fps] 


2,7 


Min Bildrate [fps] 


0,15 



Tab. A.3. Parametersatz für das Abbildungsmodell. Falls nicht anders angegeben, gelten 
für die Vorwärts- bzw. Rückwärtsabbildung diu Standardwerte, die sieh auf den optischen 
Detektor beziehen. 



Bezeichner 


Beschreibung 


Standardwerte 


N [px] 


Breite des Sensordatensatzes 


512 


M [px] 


Höhe des Sensordatensatzes 


1000 


a [mm] 


Abstand der Sensorpixel 


0,048 


d [mm] 


Abstand der Mikrolinsen 


0,48 


f [mm] 


Fokuslange des MLA 


2,4 


D [mm] 


Lochdurchmesser der Septenmaske 


0.46 


fe. 


Koordinaten-Offset (x-, y-Richtung) der optischen Zentren in 


0 




Relation zum Abstand der Mikrolinsen 


D 
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50.06 




Abb. A.l. Vereinfachte Konstruktionsskizzc des oberen Detektorrahinens (hier Detektor 4, 
Maße in [min]). Bezogen auf den Ursprung sind die Mittelpunkte für die auszufräsenden 
Rechtecke für das MLA (1) und die Septenmaske (2) angegeben. Beide Komponenten werden 
linksbündig eingesetzt, da sich auf der rechten Seite später die Onboard-Elektronik des 
Sensors befindet, für di. eine zusätzliche Aussparung (3) vorgesehen ist. Die Passtolerauz 
beträgt auf eine Länge etwa 0,05 nun. Die Maße für die Zapfen- bzw. Schraublöcher sind 
weggelassen. Septenmaske und MLA sind nach dem Fixieren gegeneinander justiert. Über die 
Zapfenverbindung wird der untere Detektorrahinen (Abb. A.2), in dem der Sensor sitzt, exakt 
gegen das MLA und die Septenmaske ausgerichtet. 
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70,16 




Ansicht 




Abb. A.2. Vereinfachte Konstruktionsskizze des unteren Detektorrahmens (Maße in [mm]). 
Dur Mittelpunkt des auszfräsenden Rechtecks für den Sensor (1) bezieht sieh auf den Ursprung, 
der über die Zapfenverbindung deckungsgleich mit dein des oberen Detektorrahiuens ist. Die 
Maße für die Zapfen bzw. Schraubenlöcher und für den Kabeldurehgaug (2) sind weggelassen. 
Zusätzlich ist hier noch ein Winkelclement für die spätere Montage an der Rotationsplattform 
vorgesehen. Über die Zapfenverbindung lässt sich dann ein robuster, justierter Detektor 
zusammensetzten und mittels der Schraubverbiudungeu fixieren. 
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Abb. A.3. Winkelpositionen der Detektoren auf der Rotationsplattform: Analog zur An- 
ordnung in Abb. 2.18(b) .sind alle vier Detektoren inklusive der zusammengefügten Detektor- 
fläche (Detektor 0) dargestellt. Die einzelnen Winkelpositionen <p werden über die Kalibrierung 
ermittelt und als Relativpositiou A<p auf den Startwiukel i^» bezogen. Eine komplette Platt- 
formrotation Aifio = 344,25° endet konstruktionsbedingt bei ^5. Der Abstand tttz bezeichnet 
die Referenzdistanz jeder SensoroberHäche zum Rotationszentruni (RZ). Die ermittelten 
Parameter sind in Tab. A.4 zusammengefaßt. 

Tab. A.4. Übersicht der Kalibrierungsparameter der Detektoren im Scanneraufbau. Pro- 
jektionen werden mit achtfachem Pixelabstand a aus dem Abbildungsmodell berechnet. Die 
Koordinatensysteme der Detektoren sind durch die Lage der Rotationsachse n Rnt (CA. 2.20) 
in den Sinograminen, die Winkelposition y? sowie einen z-Achsen-Otfset gegeben und können 
in das HKS von Detektor 0 überführt werden. Bedingt durch die Schätzung ergeben sich 
für 4 He Winkelstellung relativ große Toleranzen, die speziell für die Oberflächendelektoren 
aufgrund der lückenhaften Information stark variieren können, womit sich dann auch Ay > / 0 
erklären lässt. Für die Berechnung der Abstände zu einer Oberfläche wird als Rezugsgröße der 
Abstand t[iZ zum Rotationszentrum verwendet . Die Lage der Fokuspunkte (fpX, fpY) auf der 
Scnsoroberftä. !i ■ v. ii i einem ■ l i • t : . ■ : i - n ■ t l - 1, ■ i j Oll's. 1 relal i\ /'im Mikrolinsenabstand d 
angegeben. 





Det 1 


Det. 2 


Det. 3 


Det 4 


Det. 0 


Tol. 


Ol 

2 

'S 


Projektion 
(BxH) 


[px] 


64x128 


128x128 




Dater 


Sinogramm 
(BxH) 


[px] 


64x153 


128x153 






5,0 


[px] 


30.8 


2.8 


60,5 


29,5 


63.0 


±0.5 


rameter 


S 100 


[px] 


33.1 




62,8 


32.0 


63,0 


±0.5 


9 


[rad] 


-1.1 


0.2 


07 


2,0 


0 3 


±0 1 


■ 

CL 
01 


Acp 


[rad] 


1.4 


0,1 


0.4 


1,7 


0 


±0.1 


o> 


Al 


[px] 


2 


0 


2 


0 




±0 5 




<RZ 


|mm] 


34 


39 


39 


34 




±1 


1! 


<P» 


0.05 


0.0 


0.0 


0.08 




±0,05 


U_ CL 




0,15 


0,95 


0 05 


0,0 




±0,05 
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(a) (b) 

Abb. A.4. Konvention für den Sc-annerobjektrantn (a): Transversalebene (TE): x-y-Ebene: 
Sagit.talebene (SE): y-z-Ebene: Koronarebene (KE): x-z-Ebene. Die Detektorprojektionen 
(b) werden (Mitsprechend der Skizze in s-z-Koordinaten angegeben und mit Winkel & um die 
Rotationsachse gedreht. 
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Abkürzungen 



Abk. 


Bedeutung 


a.u. 


willkürliche Einheit {engl, arbitrary unit) 


BLT 


Biolumineszenztomographie 


CCD 


engl, charge ooupled devioe 


CMOS 


engl, complemeutary metal oxide semiconduetor 


CT 


Computertomographie 


DKFZ 


Deutsches Krebsforschungszentrum 


DKS 


Detektorkoonliiiatensyslem 


DOT 


diffuse optische Tomographie 


FBP 


gefilterte Rückprojektion (engl, filtered backprojektion) 


FMT 


Fluoreszenztomographie (engl, fluorescence mediated tomography) 


FO\ : 


Gesichtsfeld (engl, field-of-view) 


FW HM 


Halbwertsbreite (engl füll width at half maximum) 


HKS 


Haupt koordinatensystem 


LSF 


Linienbildfunktion (engl, line spread funetion) 


LTI 


engl, linear tiine-invariant 


MLA 


Mikrolinsenarray 


MRT 


Magnetresonanztomographie 


MTF 


Modulationsübertragungsfunktion (engl, modulation transfer funetion) 


PET 


Positronen- Emissions- Tomographie 


PSF 


Punktbildfunktion (engl, point spread funktion) 


SXR 


Signal-Rausch-Verhältnis (engl signal-to-noise ratio) 


SPECT 


Einzelphotonen-Emissions- Tomographie (engl. Single photon emission 




oomputed tomography ) 
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Die vorliegende Arbeit zeigt, wie sich die optische Bildgebung 
kleiner Tiere platzeffizient mit flachen Mikrolinsen-Detektoren 
realisieren lässt. Als Grundlage der Bildberechnung dient ein 
hierfür entwickeltes Abbildungsmodell. Basierend auf dem 
Detektorprinzip ist ein kompakter Kleintierscanner entstanden, 
mit dem die optische Oberflächensignalmessung für Nackt- 
mäuse durchgeführt werden kann. Die Abbildungseigenschaf- 
ten der optischen Detektoren und der Scanneraufbau sind 
anhand experimenteller Phantomstudien charakterisiert und 
bewertet worden. 
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